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PREFACE 


This  study  report  contains  the  results;  obtained  to 
date  on  the  Angel  Clutter  Reduction  Techniques  Program  conducted  by 
the  Applied  Physics  Laboratory,  Johns  Hopkins  University,  for  the 
Federal  Aviation  Administration  under  Task  I£  of  Contract  No.  DOT- 
FA-72WA-2705,  issued  on  8  October  1971. 

The  report  Includes  the  results  published  during  the 
initial  phase  of  effort  in  three  interim  reports  (May,  July,  and  August 
1972)  and  the  results  of  an  extension  phase  which  investigated  the 
feasibility  of  applying  pattern  recognition  techniques  developed  by 
Bendix  Communications  Division  to  reduce  ASR  Angel  Clutter.  Mr.  0.  E. 
Mclntire  o£  ARD-200  was  the  FAA  Technical  Representative  for  this  effort. 

The  support  of  Mr.  K.  E.  Coonley  of  ARD-200,  Mr.  C.  Chapman  and  others 
of  the  National  Aviation  Facilities  Experimental  Center  (NAFEC) ,  and 
of  the  FAA  personnel  at  General  Mitchell  Field,  Milwaukee,  is  gratefully 
acknowledged . 

This  study  report  is  organized  into  two  volumes.  Volume 
I  contains  the  Study  Results.  Chapter  1,  Summary  and  Conclusions, 
contains  all  significant  results  and  a  recommended  course  of  action 
leading  to  realization  of  an  operational  angel  clutter  reduction  capability. 
Chapter  2  discusses  the  angel  clutter  problem,  its  sources,  and  its 
effect  on  air  traffic  control  operations.  Chapter  3  identifies  differences 
in  ASR  signal  return  characteristics  for  angels  and  aircraft  which  were 
measured  at  Milwaukee  airport  during  the  spring  bird  migration  period. 
Chapter  4  describes  angel  clutter  reduction  techniques  which  (an  exploit 
These  differences  in  a  manner  which  is  cost  effective  for  the  A5R-4,  5  and 
6  radars. 

Volume  II  contains  five  appendices  providing  a  summary  of 
field  rest  operations,  supporting  dsta,  supporting  analyses,  a  discussion 
of  pattern  recognition  techniques  and  hardware  design  data  for  the 
suggested  angel  clutter  reduction  techniques. 
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APPEND  LTA-1 

ANGEL  CLUTTER  TEST  INSTRUMENTATION 
AT  MILWAUKEE 

Instrumentation  Requirements 

The  overall  objectives  for  both  the  Radar/Beacon  Multisensor 
Tracking  study  and  the  Angel  Clutter  Characteristics  study  undertaken 
by  the  Applied  Physics  Laboratory  (APL)  for  the  PAA,  presented  similar 
instrumentation  requirements .  Each  required  mobile  radar  recording 
and  processing  equipments  which  could  be  located  at  designated  FAA 
field  test  sites.  The  field  installation  of  the  instrumentation 
required  interfacing  with  operational  ASR  systems  on  an  operationally 
safe  and  non-interfering  basis. 

Specific  data  requirements  for  the  Angel  Clutter  Characteristics 
study  necessitated  the  simultaneous  recording  of  normal  and  MTI  radar 
information  at  an  operational  ASR  site,  and  digital  radar  processing 
facilities  for  on  and  off-line  daca  collection,  formatted  and  recorded 
in  a  manner  which  was  compatible  with  the  available  laboratory  digital 
data  reduction  facilities.  Additionally,  a  source  of  precision  target 
position  measurement  was  required  to  permit  a  determination  of  angel 
clutter  range,  bearing,  altitude,  motion  and  velocity.  Clutter  position 
information  was  also  required  in  a  form  which  was  compatible  with  digital 
processing.  A  pictorial  diagram  of  the  angel  clutter  field  test 
Instrumentation  installation  is  presented  in  Figure  A-l. 

Instrumentation  Description 

Applied  Physics  Laboratory  experience  in  the  research  and 
development  phases  of  U.  S.  Navy  search  and  fire  control  radar  system 
programs ,  and  the  recent  laboratory  development  of  the  U .  S .  Navy 
Detector /Tracker  Set  AN/SYS-1  target  information  processing  system 
was  instrumental  in  the  rapid  design  and  construction  of  a  mobile 
Data  Acquisition  Module  for  use  in  FAA  field  testing.  The  major 
components  of  the  module  are  shown  in  Figure  A-2.  the  module  is 
designed  to  accept  analog  ASR  radar  video  with  separate  input  lines 
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for  normal,  Mil  and  beacon  video,  triggers,  aaiauth  XX  synchro, 
and  APG  (Aaiauth  PuXaa  Ganarator)  pulaas,  aacb  with  appropriate 
isolation.  Thasa  Inputs  may  ba  racordad  directly  on  an  analog 
Radar  Video  Recording  Syatsa  (RAVIR)  or  processed  by  a  radar 
digitiser  (referred  to  as  the  Adaptive  Video  Processor  in  the 
AN/ SYS- 1  system)  and  racordad  in  digital  fora  on  magnatic  tape. 

The  off-line  node  of  the  equipment  allows  the  replay  or  previously 
recorded  analog  video  tapes  into  the  digitiser  for  processing  and 
digital  recording.  A  Display  and  Contrr/1  Console  allows  an  operator 
to  view  normal,  MTI  or  beacon  video,  processed  video  in  -the  fora 
of  fixed,  mov.'ng  and  tentative  tracks  and/or  tracking  radar  positions. 
Ancillary  equipments  in  the  module  Include  a  VWV  receiver  and  Time 
Code  Generator,  VHF  communication  transceiver,  a  Teletypewriter 
and  electronic  test  equipment.  These  equipments  are  installed  in 
a  truck  mounted,  climate  controlled  8  x  16  foot  van.  The  Data 
Acquisition  Module,  AN/MPSr,l9  Tracking  Radar  and  Generator  van  were 
located  at  General  Mitchell  Field,  Milwaukee,  Wisconsin  during  the 
Angel  Clutter  Characteristics  testing. 

In  addition  to  the  instrumentation  provided,  a  Cessna  172, 
single  engine  aircraft  was  employed  for  the  collection  of  typical 
light  aircraft  radar  data.  These  data  were  collected  under  the 
same  conditions  as  the  angel  clutter  data. 

Test  Coordination  and  lnstrumentatxuu  uiuaw 

Data  collection  was  coordinated  by  a  Test  Conductor  positioned 
a t  the  Date  Acquisition  Module  Display  and  Control  Console.  An 
intercom  link  with  the  MPS-19  Tracking  Radar  permitted  voice  designation 
to  ASR  detected  targets  of  interest.  Television  monitors  in  the  MPS-19 
and  at  the  Display  and  Control  Console  provided  a  visual  display  of 
the  target  being  tracked  by  the  tracking  radar.  This  information  was 
obtained  from  a  television  camera  which  was  attached  to  the  MPS-19 
antenna  and  boresited  with  the  radar  beam.  (Aircraft  targets  were 
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generally  clearly  visible  on  the  television  system,  however,  the 
limited  power  of  the  system  precluded  visual  observation  of  angel 
clutter  being  tracked).  During  controlled  aircraft  data  collection, 
the  Test  Conductor  directed  the  aircraft  runs  via  the  VHP  communi¬ 
cation  facilities  provided  in  the  Data  Acquisition  Module. 

Angel  Clutter  Test  Instrumentation,  Detailed  Description 

Data  Acquisition  Module 

The  Data  Acquisition  Module  is  a  mobile  and  highly 
versatile  instrumentation  van  which  provides  the  facilities  for 
the  recording/reproduction  and  display  of  analog  radar  data;  the 
digital  processing  and  recording  of  radar  data  with  target  returns 
categorized  as  fixed  or  moving  ar.d  providing  specific  information 
such  as  video  amplitudes,  hit  counts,  run  lengths,  etc.  A  special 
mode  of  operation  allows  the  collection  of  data  in  a  finite  volume 
of  space  with  discrete  measurements  of  video  surrounding  a  designated 
point,  e.g.,  clutter  surrounding  an  aircraft  target,  the  distribution 
of  video  In  a  thunder  cell,  the  distribution  of  discrete  targets  in 
a  flock  of  birds,  etc.  These  unique  processing  and  recording  facili¬ 
ties  are  supported  by  subsystems  and  equipment  necessary  for  remote 
all-weather  operation. 

A  pictorial  diagram  of  the  Data  Acquisition  Module  was 
shown  in  Figure  A-2.  The  major  components  of  the  module  are  discussed 
in  the  following  paragraphs. 

RAVIR  is  a  wideband  radar  signal  recording  system  which 
was  designed  and  developed  by  the  Applied  Physics  Laboratory  for 
use  by  the  U.  S.  Navy  as  a  radar  training  device.  Seven  magnetic 
tape  recording  channels  are  available  with  interchangeable  input 
signal  conditioning  units.  This  flexibility  allows  the  simultaneous 
recording  of  normal  video,  MTI  video,  raw  defruited  beacon  video, 
radar  triggers,  APG  pulses,  IX  synchro  and  timing. 
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A  listing  of  RAVIR  recording  perimeters  is  provided  in 

Table  A-l . 

TABLE  A-l  RAVIR  RECORDING  PARAMETERS 


General  Tape  Recorder  Characteristics 

Parallel  Channels 

7 

Frequency  Response  (3  dB) 

400  B*  to  1.5  MBs 

Record  Time 

24  min. 

Tape  Speed 

120  inches/sec. 

Tape  Speed  Error  (Short  term) 

-  0.1Z  max. 

Dynamic  Skew  (adj .  channels) 

-0.3  micro  sec. 

Signal-Noise  Ratio 

23  dB  max. 

Crosstalk 

20  dB  min. 

RAVIR  Accuracy 

Bearing 

4. 

-  0*5  degrees 

Range 

1Z 

Output  Trigger  Delay 

0.25  microsec  max 

Frequency  Response 

•  DC  to  1.5  MHz 

During  radar  signal  processing,  the  normal  and  MTI  video 
signals  are  multiplexed  to  allow  MTI  video  utilization  out  to  a 
manually-selected  transition  range  and  normal  video  beyond  that 
range.  The  result  of  this  video  selection  is  directed  to  the 
Target  Information  Processing  System  (TIPS)  which  contains  the 
Adaptive  Video  Processor  (AVP) .  The  AVP  employed  in  the  module  was 
designed  for  a  high  level  of  effectiveness  in  a  land  environment  and 
in  rain  clutter .  The  processor  samples  the  video  amplitude  in  the  ten 
range  cells  (each  cell,  being  equal  to  the  pulsewidth  of  0.S33  usee) 
preceding  and  ten  cells  succeeding  the  ceil  being  tested  for  the  presence 
of  a  target.  These  twenty  cells  are  summed  to  determine  a 
mean  clutter  level.  The  cell  being  tested  is  compared  to  this  average 
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level  end  If  the  test  cell  exceeds  the  level  by  e  specified 
factor,  a  target  is  declared  present.  The  data  storage  is  advanced 
by  one  cell  and  the  process  is  repeated.  Thus,  the  output  of  this 
portion  of  the  digitizer  is  a  string  of  ones  and  zeroes,  where  the 
ones  indicate  the  existence  of  a  possible  target  and  the  occurrence 
in  time  indicates  the  range.  This  string  of  data  is  then  stored  in 
a  serial  delay  line  which  is  segmented  ao  that  each  segment  corresponds 
to  the  binary  data  obtained  from  one  azimuth  dwell.  Nineteen  such 
azimuth  strings  are  stored  and  at  each  range  cell  the  output  of 
all  nineteen  are  sampled.  If  at  least  m  out  of  the  nineteen  (n) 
bits  are  one,  then  a  target  "hit"  is  declared.  This  process  is 
referred  to  as  azimuth  correlation  or  "sliding  window  correlation". 

The  output  of  this  correlator  is  used  to  transfer  the  range  of  the 
detection  into  the  data  converter  which  formats  data  for  the  digital 
data  processor. 

In  addition  to  processing  all  targets  ns  Indicated  above, 
it  is  possible  to  collect  concentrated  data  about  a  single  target 
of  interest,  referred  to  as  a  "hooked"  target.  When  a  specific 
target  is  designated  (by  the  Test  Conductor) ,  the  digital  data 
processor  generates  a  grid  about  the  target  and  requests  the  ampli¬ 
tude  of  the  radar  returns  for  all  cells  within  the  grid  from  the 
digitizer.  The  grid  consists  of  twenty  cells  in  range  and  twenty 
azimuth  positions,  although  they  need  not  be  adjacent  azimuth  dwells. 
These  samples  are  converted  to  5  bits  of  amplitude  data  and  loaded 
into  the  digital  processor .  The  characteristics  of  the  digitizer 
are  listed  in  Table  A-2. 
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TABLE  A-2 

ADAPTIVE  l  jDEQ  PROCESSOR  (DIGITIZER) 
CHARACTERISTICS 


-  -  - 

- 

Delay  Line  Length 

20t  (t  *  0.8  microsec) 

Delay  Line  Spacing 

t/2  or  T 

Signal  Tap 

10t 

Average 

estimate  of  rms  noise 

Azimuth  Correlation 

Length  (n) 

variable,  1-19  (n  2  m) 

Minimum  Required  Azimuth 
Hits  (m) 

variable,  1-19 

In  addition  to  the  features  described  above,  a  constant 
false  alarm  loop  is  incorporated  around  the  signal  detection  threshold 
to  maintain  a  constant  false  alarm  rate  (CFAR) ,  as  averaged  over 
each  azimuth  position.  The  desired  false  alarm  can  be  manually  set 
at  0.5Z,  1.0Z,  5. OX,  10X,  or  20X.  A  vide  pulse  discriminator  (WPD) 
is  also  incorporated  following  the  threshold  detector.  Basically, 
the  assumption  is  that  an  aircraft  could  be  detected  in  two  adjacent 
range  cells,  but  due  to  the  length  of  aircraft  and  the  allowable 
minimum  spacing,  aircraft  detections  cannot  occur  in  three  or  more 
adjacent  range  cells.  In  strong  clutter,  however,  this  condition 
is  possible.  Consequently,  the  WPD  inhibits  detections  when  three 
or  more  consecutive  range  cells  experience  detections.  The  WPD 
can  be  manually  bypassed,  if  desired. 

The  Data  Converter  is  an  interface  unit  between  the  AVP 
and  the  digital  data  processor,  and  between  the  radar  and  the  processor. 
It  contains  the  range  counter  for  specifying  the  location  of  potential 
targets.  The  least  significant  bit  of  this  counter  corresponds  to 
approximately  1/16  of  a  nautical  mile.  (In  the  data  collect  mode 
(hsoked  target)  the  range  of  the  grid  points  are  specified  to  about 
1/64  of  a  nautical  mile) .  Angle  measurements  to  an  accuracy  of 
0.022  degrees  are  converted  to  the  appropriate  format  and  loaded  into 
the  digital  data  processor  when  Che  AVP  azimuth  correlator  declares  a  hi 
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The  Digital  Data  Processor  is  *  liuneyvell  DDF-516R 
computer.  The  inputs  to  the  processor  are  the  hits  as  declared 
by  the  AVP  azimuth  correlator,  MPS-19  tracking  radar  and  timing 
data.  Typically,  twenty  to  £orty  successive  hits  in  azimuth  would 
be  received  via  th«  AVP  £rom  a  large  aircraft  target.  The  processor 
shows  this  as  one  target  and  distinguishes  between  adjacent  targets 
by  the  absence  of  at  least  twelve  hits  in  azimuth.  For  the  hit 
group  assigned  to  a  single  target,  the  position  halfway  between 
the  first  and  last  hits  (centroid)  is  defined  as  the  target  center. 

A  detailed  description  of  the  process  involved  in  the  establishment 
of  target  files,  classification  of  targets,  modification  of  target 
classifications,  and  criteria  for  dropping  tracks,  etc.  is  beyond 
the  scope  of  this  description.  In  general,  however,  new  hits  which 
do  not  correlate  with  existing  tracks  are  established  as  net'  tentative 
tracks.  Depending  upon  the  number  of  successive  detections,  the 
number  of  misses  and  the  computed  velocities,  these  new  tentative 
targets  are  promoted  to  tentative,  firm,  or  fixed  tracks.  As  tracks 
develop  further,  the  classifications  may  change  and  some  tracks  may 
be  dropped  if  insufficient  detections  occur. 

Once  each  complete  ASR  antenna  rotation  the  processor 
dumps  the  contents  of  all  track  files  onto  magnetic  tape  along  with 
a  real-time  reading  which  has  been  received  from  the  module  time  code 
generator. 

Table  A-3  lists  the  significant  characteristics  of  the 
data  processor. 
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TABLE  A-3 

DDP-516R  DIGITAL  DATA  PROCESSOR 
CHARACTERISTICS 


Memory  Size 

8K  words 

Word  Length 

16  bits 

Memory  Cycle  Time 

1  microsec 

Direct  Memory  Access  Channels 

3 

Add/Sub tract  Time 

2  microsec 

Divide  Time 

11  microsec 

Multiply  Time 

6  microsec 

The  Control  end  Display  Console  serves  as  the  primary  means 
for  monitoring  the  performance  of  the  digitizer  and  digital  data 
processor.  Using  a  PP1  type  display.,,  the  operator  may  select  analog 
normal/MTI  video  from  the  ASR  in  real-time  or  KAVIR  playback  video, 
or  he  may  select  processed  video  with  the  appropriate  symbology  for 
tentative,  firm  and  fixed  tracks,  and/or  a  track  symbol  indicating 
the  tracking  position  of  the  MTS-19  radar.  An  azimuth  display  and 
range  rings  are  also  available.  The  track  ball  feature  of  the  console 
allows  the  operator  to  "hook"  a  target  of  interest  to  place  the  system 
in  the  data  collect  mode  previously  described.  The  "hooked  track" 
feature  at  the  console  is  primarily  employed  during  system  calibra¬ 
tion  and  test,  or  when  the  module  is  operated  at  remote  locations. 

The  "hooked  track"  is  initiated  by  positioning  video  cross  hairs 
over  the  target  of  Interest  with  a  track  ball,  thereby  positioning 
the  data  collect  aatrix.  Depressing  a  thumb  operated  switch 
activates  the  data  collect  mode. 

The  MPS-19  was  selected  as  the  instrumentation  tracking 
radar  due  to  the  near  identical  characteristics  to  the  ASR  series 
characteristics.  These  are  listed  in  Table  A-A. 
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TABLE  A- 4 

RADAR  CHARACTERISTICS 


ASR 

MPS-19 

Frequency 

2700-2900  MHz 

2700-2900  MHz 

Peak  Power 

400  kU 

137  kU  -  250  kW 

PRR 

1200  pps 

300-2000  pps,  adjustable 

Pulsevidth 

0.833  Us 

0.8  us 

Antenna  Beamvidth  Az. 

1.5° 

3.0° 

El. 

30° 

3.0° 

Antenna  Gain 

34  dB 

34  dB 

IF  Bandwidth 

2  MHz 

2  MHz 

An  analog  to  digital  converter  was  fabricated  and 
installed  in  the  MPS-19  to  allow  direct  tracking  data  input  to  the 
Data  Acquisition  Module  Data  Processor.  The  tracking  radar  was 
also  instruaented  for  the  recording  of  the  AGC  voltage  to  permit 
the  determination  of  target  aixe. 
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SUMMARY  OF  MILWAUKEE  DATA 
COLLECTION  OPERATIONS 

The  Milwaukee  field  operations  were  conducted  April  12-20 
and  May  10-12,  1972,  The  instrumentation  consisted  of  the  Data 
Acquisition  iiodule,  the  AN/MPS-19  Track  Radar  Module,  and  a  150  KW 
generator  van  previously  described  in  Appendix  A-l. 

In  addition  to  the  recorded  data  described  below,  observations 
were  made  and  local  reports  were  collected  cn  weather  and  bird  activity, 

ASR  Radar  Video  Data 

ASR  normal  and  MTI  video  was  recorded  on  :tfc?  Radar  Video 
Recording  System  (RAVIR)  installed  in:  the  Data  Acquisition  Module. 

It  was  also  used  in  the  off-line  mode  to  generate  digital  magnetic 
tapes  via  the  Target  Information  Processing  System  (TIPS)  portion  of 
the  Data  Acquisition  Module, 

The  TIPS  was  used  to  process  and  record  several  types  of  data 
which  would  be  required  in  subsequent  reconstruction  of  the  Milwaukee 
radar  environment.  These  included  angel  clutter,  land  and  lake  clutter, 
light  to  moderate  rain  clutter,  heavy  weather  cell  clutter,  specific 
aircraft  track  histories  and  the  Milwaukee  ASR  environment  under  various 
traffic  conditions 

A  high  speed  data  collect  mode  was  employed  primarily  in  the 
collection  of  angel  clutter  data.  In  this  mode,  system  tracking  of 
all  ASR  targets  is  suspended  and  replaced  by  an  expanded  data  collect 
mode  wherein  the  data  collect  matrix  center  is  positioned  by  the  target 
track  coordinates  of  the  Track  Radar  Module  This  enables  the  measure¬ 
ment  and  recording  of  the  position  and  amplitude  of  all  radar  returns 
within  the  matrix  area. 


A-13 


APPENDIX  A- 2 


A  low  speed  data  collect  node  was  utilised  primarily  in  the 
accumulation  of  ASH  track  histories.  These  histories  were  recorded 
during  various  terminal  traffic  condition*  and  in  various  weather  and 
AP  clutter  conditions. 

Track  Radar  Module  Data 

The  AN/MPS-19  Track  Radar  Module  contained  in  addition  to 
the  tracking  radar,  a  television  monitor  and  recording  system,  a 
strip  chart  recorder  used  to  record  radar  automatic  gain  control 
(ACC)  voltages,  and  an  analog  SONY  recorder  used  to  record  track 
radar  video. 

The  bore-sited  TV  system  was  Intended  as  a  means  of  Identifying 
bird  angels  in  track.  However,  this  was  Jrwf'"*-  « Jvo  due  to  a  combination 
of  darkness  and  excessive  range  for  visual  detection,  even  though  a 
high-power  lens  was  used. 

During  angel  or  aircraft  tracking,  the  A/D  converter  track  radar 
outputs  were  transmitted  to  the  Data  Acquisition  Module  in  the  high 
speed  data  collect  mode.  Track  range,  azimuth,  and  elevation  were 
recorded  every  0.4  seconds,  and  the  eight  bit  (255  level)  radar  video 
or  ACC  level  was  recorded  for  every  radar  pulse .  To  avoid  Interference 
with  the  ASR,  the  track  radar  transmissions  were  synchronized  with  the 
ASA  pulse  repetition  frequency  (PRF) .  Continual  use  of  the  ASR 
stagger  led  to  variations  in  track  radar  output  power  because  of 
uneven  modulator  charging  times;  hence  the  MPS-19  PRF  was  reduced  to 
approximately  575  Hz  ;»nd  was  maintained  in  synchronism  with  every 
other  ASR  transmit  pulse 

Data  Library 

Table  A-4  provides  a  summary  of  all  Milwaukee  digital  and 
analog  data  tapes  collected.  Table  A-5  provides  a  detailed  listing 
by  date  and  time  Of  the  data  collected. 


ASR  Environment 
Systen  Calibration 
AP  Track  ^rtory 
AP  MPS-19  Track* 


14 

4 

,5 

If) 


0 

2 

3 

9 


0 

2 

0 

4 


0 

4 

0 

9 


Total  No.  of  Tapes 


33 


14 


13 


Total  Minutes  (Approx) 


700 


330 


420 


450 


*  A  strip  chart  recording  cf  MPS-19  receiver  AGC  was  made 
during  all  runs  when  this  radar  was  used 
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SUMMARY  OF  PRELIMINARY  ANGEL  TRACKING 
TESTS  AT  NAFEC,  OCTOBER  1971 


Background 

Task  II  of  the  APL  contract  with  the  Federal  Aviation  Agency 
is  "Angel  Clutter  Reduction  Techniques  "  The  objective  of  this  effort 
is  to  evaluate  the  effects  of  angel  clutter  on  ASR  radars  and  to 
determine  cost-effective  means  of  reducing  these  effects.  Angel  echoes 
can  be  caused  by  meteorological  inhomogeneities ,  birds,  or  even  insects* 
Previous  experiments  have  indicated  that  the  majority  of  angel  echoes 
are  due  to  birds,  and  the  purpose  of  this  field  operation  was  to 
measure  the  detailed  characteristics  of  bird  angels  during  the  fall 
bird  migration  at  the  FAA’s  NAFEC  facility  near  Atlantic  City,  New 
Jersey.  Although  wideband  video  recorder  tapes  of  ASR  video  had  been 
made  previously  during  the  April  1971  migrations,  this  field  operation 
was  the  first  to  record  tracking  radar  data  as  well  as  ASR  data. 

Continuous  track  data  is  required  to  fully  assess  the  detailed  differences 
between  the  radar  target  characteristics  of  bird  angels  and  aircraf  t, 
and  to  develop  appropriate  methods  of  differentiating  between  them 
with  the  ASR. 

Instrumentation 

Figure  A-3  indicates  the  equipment  used  during  these  tests. 

A  wideband  video  recorder  (RAVIR  Mod  2A)  was  installed  in  the 
NAFEC  test  facility  to  simultaneously  record  MTI  and  normal.  ASR  video 
and  triggers.  A  tripod-mounted  Polaroid  camera  was  also  employed  to 
record  single  and  multiple-scan  PPI  video  presentations.  Twc  portable 
equipment  modules  were  parked  nearby:  a  150  KW  prime  power  generator 
and  the  Tracking  Radar  Module.  This  module  consists  of  an  S-Band 
(AN/KPS-19)  tracking  radar,  a  TV  system,  a  Sanborn  strip  chart  recorder 
to  record  ACC  data,  and  a  SONY  PV  120V  video  tape  recorder,  which 
was  used  to  record  'rack  radar  video  and  voice  annotations  of  the 
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type  of  target  end  its  identifying  characteristics  and  positional  data 
(range,  bearing,  and  elevation  as  a  function  of  tine).  The  track 
radar  transmitter  was  synchronised  with  the  A5R  transmit  pulse  to 
reduce  mutual  Interference. 

The  S-band  tracking  radar  had  been  received  on  loan  only 
three  weeks  earlier  from  the  Air  Force.  This  precluded  the  intended 
full  refurbishment  of  the  radar  prior  to  the  date  of  the  bird  migrations; 
however,  it  was  decided  to  conduct  the  operation  at  this  time  rather 
than  to  postpone  all  operations  until  the  spring  migrations.  The 
radar  on  the  whole  performed  well  during  the  test  period  (due  tc  a 
strong  effort  by  the  test  crew)  and  was  adequate  for  gathering  the 
desired  data. 

Chronology 

The  Track  Radar  Module  and  Prime  Power  Generating  Module 
arrived  at  NAFEC  late  in  the  evening  of  26  October  and  the  Track 
Radar  Module  was  in  operation  the  following  day.  The  wideband  video 
recording  system  was  mated  to  the  ASR-5  tc  record  both  normal  and 
MTI  video  simultaneously.  Little  or  no  migratory  bird  activity 
was  observed  on  the.  evening  of  26  October,  primarily  due  to  very  poor 
weather  conditions  (fog) .  A  combination  of  fog  and  night  precluded 
use  of  the  TV  system  to  gather  visual  data  on  bird  tracks  throughout 
most  of  the  operation. 

Early  in  the  evening  of  28  October,  fog  moving  in  from  the 
south  precipitated  very  heavy  bird  activity.  Bad  weather  normally 
precludes  migratory  activity  and  local  reports  had  indicated  considerable 
numbers  of  birds  in  the  area  (primarily  at  the  Brigantine  Wildlife 
Refuge)  waiting  for  clearing  conditions.  The  fog  on  the  evening  of 
28  October  was  generally  low-  lying  and  the  moon  could  be  clearly  seen 
above  it,  indicating  acceptable  migratory  weather.  Bird  angel  activity 
began,  shortly  after  dark  and  by  1930  EST  had  reached  lea*'"  fflcc 


A-21 


APPENDIX  A- 3 


A  series  of  five  15-minute  wideband  video  recordings  of  ASR  video 
were  made  between  the  hours  of  1541  and  2017  (Table  A-6) .  Almost  thirty 
Polaroid  photos  were  taken  between  1607  and  2150  hours  (Table  A-7) ; 

Figures  A-4  and  A- 5  show  representative  single  and  multiple-scan  PPI 
presentations  for  MTI  video  with  and  without  STC.  It  is  clear  from 
these  photographs  that  neither  MTI  nor  STC  is  a  complete  solution  to  the 
angel  problem  during  periods  of  hervy  bird  activity. 

The  Tracking  Radar  Module  PPI  display  presented  a  similar  picture. 
Aircraft  and  apparent  bird  angels  were  tracked  to  obtain  motional  character¬ 
istics  and  recordings  of  radar  video  and  ACC  voltage.  Figure  A-6  shows 
three  representative  ACC  recordings  and  the  A-scope  video  corresponding 
to  each. 

Since  all  required  data  were  recorded  on  the  evening  of  28  October, 
the  Tracking  Radar  Module,  generator  van,  and  wideband  video  recording 
equipment  were  returned  to  APL  on  29  October. 

Discussion  of  Preliminary  Results 

It  was  possible  to  make  several  statements  based  upon  the  data 
gathered  during  this  operation: 

(1).  ACC  recordings  of  tracking  radar  video  permitted  a  clear 
distinction  between  the  angel  and  aircraft  tracks.  (Aircraft 
were  identified  by  their  much  higher  velocity.)  This  difference 
is  clearly  visible  from  Figures  A-6 (A)  and  A~6(B).  It  should  be 
noted,  however,  that  the  AGC  output  is  smoothed  data  and  such 
distinctions  may  not  be  apparent  in  the  short  period  of  time 
during  which  the  ASR  scans  past  the  target . 
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MTI  VIDEO,  STC,  SINGLE  SCAN 


FIGURE  A— 4 

SINGLE  SCAN  ASR  PPI  PRESENTATION 
DURING  HEAVY  BIRD  ACTIVITY 
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TABLE  A-6 

WIDEBAND  VIDEO  IECOPDINGS  TAKEN  AT  NAFEC 
28  OCTOBER  1971 


Recorder;  BAVIR  Mod  2A 


Channel  1 

Normal  Video 

Track  2 

Chennel  2 

Mil  Video 

Track  4 

Chennel  3 

MPS-19  Triggers 

Track  5 

Timing 

ta  delayed  (-94  psee) 

Track  3 

T«p«  tei 

Number 

Start 

Stg£ 

28-1341 

1541  hre. 

1601  hre. 

28-1649 

1649 

1710 

28-1836 

1836 

1857 

28-1917* 

1917 

— 

28-2000* 

2000 

2017 

Broken  tepee.  #28-1917  not  rewound 
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TABU  A-7 

ASR-5  VIDEO  PHOTOGRAPH  LOG 
28  OCTOBER  1971 


PHOTO  APPROX.  RADAR  VIDEO 


MUtOlR 

=  /  /  - 

TUB 

NORM 

MTI 

STC 

DURATION  REMARKS 

1 

1607  hra. 

X 

1  acan  Caaara  Teat 

2 

1608 

X 

1  scan  Caaara  Teat 

3 

1609 

X 

1  scan 

4 

X 

1  scan 

3 

X 

60  sac 

6 

X 

120  aac 

7 

1634 

X 

300  sac 

8 

1810 

X 

1  scan 

9 

1818 

X 

240  sac 

10 

1821 

X 

90  aac 

11 

1830 

X 

X 

1  scan 

12 

1831 

X 

X 

90  aac 

13 

1834 

X 

1  scan  'V 

14 

X 

90  sac  \  In  Saquanca 

13 

X 

1  scan  J 

16 

X 

1  scan 

17 

1923 

X 

1  scan  ^ 

a 

11 

1925 

X 

90  sac  1 

19 

1926 

,  )  In  Sequence 

1  scan  f 

20 

1927 

X 

10  scans J 

21 

1930 

X 

1  scan 

22 

1931 

X 

X 

1  scan 

23 

X 

X 

1  scan 

24 

X 

X 

15  scans 

23 

X 

1  scan  "I  jn  ggqU#ncg 

26 

2140 

X 

60  sac  J 

27* 

2110 

X 

- 

1  scan  (20  Mila  Rang*) 

*20  ail*  radius;  remainder  lC-*ila  radius  with  2  alia  rang*  ring*. 
Oaa  aean  •  A  seconds 
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(2)  Visual  observation  of  the  A-scope*  video  (which 
contains  the  high  frequency  amplitude  fluctuation  components 
not  visible  on  the  AGC  recordings)  also  indicated  different 
fluctuation  characteristics  for  angel  and  aircraft  targets , 

Tho  observers  were  able  In  a  number  of  cases  to  recognize 
aircraft  and  bird  angels  froa  the  A-scope  video,  althougn 
the  AGC  recordings  appeared  to  be  a  more  reliable  indicator , 
There  was  a  general  feeling  that  the  bird  fluctuations  could 
be  associated  with  wing  beats,  although  analysis  of  the  video 
tape  recordings  are  required  to  determine  if  these  fluctuations 
corresponded  to  wing  beats  or  the  fluctuations  associated 
with  many  birds  in  the  same  radar  resolution  cell. 

(3)  TWo  distinct  types  of  angel  returns  were  identified. 

The  majority  were  of  the  type  shown  in  the  AGC  recording 
of  Figure  A-6(A).  The  type  shown  in  Figure  A-6(C)  were 
encountered  during  attempts  to  determine  the  maximum  altitude 
of  the  engil  activity.  These  angels  were  Isolated  echoes 
with  0.5  Hz  AGC  fluctuations  and  were  encountered  at  altitudes 
as  high  as  6600  feet. 

(4)  While  MTI  provided  much  improvement  ever  normal  video, 
and  MTI  with  STC  was  even,  better,  a  large  number  of  angel 
echoes  remained.  It  is  estimated  that  there  are  about  250 
angel  targets  on  the  MTI -STC  single  scan  photograph  in  F'gvr.- 
A-4  (excluding  the  automobile  t-affic  visible  on  a  nearbv 
road) . 


*  A  tracking  radar  A-scope  displays  signal  amplitude  (received  videc’ 
vertically  end  range  (time)  horizontally;  the  horizontal  sweep  moves 
to  the  right  with  increasing  range. 
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(5)  It  ./as  possible  to  distinguish  some  aircraft  returns 
from  angel  echoes  on  the  10-mile  PPI  display  on  the  basis 
of  angular  width  (run  length),  particularly  on  the  MTI-STC 
display  where  angel  clutter  was  somewhat  reduced.  This 
c;sn  be  seen  at  63°  in  Figure  A-5  (#24).  This  large  air¬ 
craft  produced  a  4°  bearing  width,  indicating  detections 
in  the  major  sldelobes  cf  the  AS&  radar  (1.5°  besmwldth) . 
This  would  not  be  the  case  for  smeller  aircraft. 

(6)  Preliminary  data  analysis  Indicated  that  the  aircraft 
targets  encountered  had  radar  cross-sections  of  from  25 

to  1000  square  meters,  while  the  angel  returns  were  in 
the  range  0.0001  to  0.01  square  meters.  With  somewhat 
limited  data  reduction,  it  was  not  possible  to  determine 
the  proportions  of  the  angel  echoes  due  to  single  birds 
and  groups  of  birds. 
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ANGEL  CLUTTER  TEST  SITE  SELECTION 

Fallowing  the  laitial  angel  clutter  teats  at  NAFEC  and  the 
refurbishment  and  optimization  of  field  test  Instrumentation,  It 
was  necessary  to  select  a  test  site  for  angel  clutter  data  collection 
during  the  spring  bird  migrations.  Criteria  for  this  selection  were: 

(a)  A  history  of  migratory  or  local  bird  activity  in  the 
area, 

(b)  a  history  of  angel  clutter  activity  with  some  level  of 
impact  on  Air  Traffic  Control  functions  and  coincident 
with  bird  activity  in  the  area,  and, 

(c)  a  physical  situation  at  the  radar  site  which  was 
appropriate  to  the  location  and  installation  of  the 
necessary  instrumentation. 

Subsequent  to  a  review  of  FAA  angel  clutter  field  reports,  bird  migration 
and  weather  forecasts,  and  considering  the  test  period  available,  four 
candidate  sites  were  selected.  These  were  Little  Rock,  Arkansas; 

Memphis,  Tennessee;  Minneapolis,  Minnesota  and  Milwaukee,  Wisconsin. 
Following  visits  to  each  of  the  airport  FAA  facilities  by  APL  personnel, 
General  Mitchell  Field,  Milwaukee  was  selected. 

Details  of  the  discussions  and  observations  at  the  four 
locations  are  provided  in  the  two  APL  trip  reports  Included  in  this 
Appendix. 
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1 •  REPORT  ON  VISITS  TO  LITTLE  ROCK  AND  MEMPHIS 
March  1972 


1,  Introduction 


APL  personnel  visited  Little  Rock  Municipal  Airport  (Adana 
Field) ,  Little  Rock,  Arkansas,  on  March  8-9  and  Memphis  International 
Airport,  Memphis,  Tennessee,  on  March  10,  The  purposes  of  the  visits 
vere  to  observe  and  discuss  local  angel  clutter  problems  and  to  evaluate 
the  potential  of  the  two  locations  for  Task  II  data  collection.  Little 
Rock  was  considered  a  likely  site  for  the  planned  field  measurement 
program;  therefore,  essential  preliminary  coordination  and  planning 
was  carried  out  with  FAA  personnel  there,  Robert  Stutsman,  SW-432, 

Task  II  Coordinator  for  the  Southwest  Region,  took  part  in  discussions 
of  Little  Rock  on  March  9-  This  report  summarizes  the  observations 
and  activities  of  the  trip  and  presents  recommendations  on  Little  Rock 
as  a  data  collection  site.  Key  contacts  in  Little  Rock  are  listed  in 
Enclosure  1.  A  brief  investigation  eliminated  Memphis  as  a  possible 
site  for  Task  II  measurements  (Section  8) • 

2,  Summary 

A  visit  to  Little  Rock  Municipal  Airport  on  March  8-9, 

1972  established  that  air  traffic  control  thei'e  is  significantly 
constrained  twice  daily  in  October-Ncvember  and  February-March  by  radar 
clutter  (Figures  l-4)*and  that  this  clutter  is  caused  by  birds  which 
roost  near  the  airport.  Frequently-observed  "clear-air"  clutter  (Figure  5) 
was  reported  to  be  as  serious  a  problem  as  the  bird  clutter.  The  nature 
of  the  local  problem,  the  physical  layout,  the  facilities,  the  air 
traffic  load,  and  the  local  FAA  personnel  make  Little  Reck  an  excellent 
site  for  a  field  measurement  program.  However,  Little  Rock  bird  experts 
(Enclosure  1)  predict  that  the  birds  which  cause  the  radar  clutter  will 
start  leaving  that  area  around  March  15  and  that  no  concentrations  are 
likely  to  remain  after  April  1, 

The  probable  bird  migration  schedule  together  with  the 
unlikelihood  cf  being  able  to  start  data  collection  m  Little  Rock  before 
the  last  week  of  March  prompts  the  following  tVccnanendations: 

Most  desirable  -  Locate  a  suitable  airport  site  in  the 

northern  Mississippi  Valley  where  birds  will 
be  common  in  April; 


*  AIL  figures  are  grouped  at  the  end  of  the  text. 


A—3 1 


APPENDIX  A-4 


Alternative  -  Begin  data  collection  In  Little  Rock 

as  soon  as  possible,  no  later  than 
the  last  week  in  March.  If  birds 
are  gone,  or  when  they  go,  collect 
data  on  clear-air  clutter  and  other 
local  angel  clutter. 


A  visit  to  Memphis  International  Airport  established  that 
terminal  air  traffic  control  there  is  not  affected  by  birds  nor.,  to 
any  significant  extent,  by  any  other  forms  of  "anomalous  propagation." 
Therefore,  Memphis  is  considered  to  have  no  value  as  a  data  collection 
site  for  FAA  Task  11.  Memphis  birds  migrate  at  the  same  time  as  Little 
Rock  birds. 


3.  Radar 


Little  Rock  ATC  has  an  ASR  4B  radar  with  parametric  amplifier 
and  PIN-dlode  attenuator  The  effect  of  the  latter  wss  described  as 
reducing  sensitivity  to  approximately  the  level  obtained  without  'the 
parametric  amplifier.  Antenna  height  was  estimated  to  be  45  feet. 

Antenna  tilt  was  reported  as  +2® . 

Standard  operating  procedure  is  to  use  MT1  video  with  SIC 
from  0  to  15  miles  and  NORMAL  video  from  15  to  40  miles.  NORMAL  video 
is  useless  inside  10  miles  due  to  ground  clutter  With  MTT,  no  ground 
clutter  appears.  Surface  traffic  is  observed  (Figure  4a,  240%  9  and 
14  miles)  and  recognized  by  location  and  speed,  and  then  is  ignored^ 
without  difficulty.  Precipitation  returns  are  reduced  by  use  of  circular 
polarization  (CP);  attenuation  is  satisfactory  in  light  rain  and  is 
useful  but  inadequate  in  heavy  rain.  When  AP  and  bird  clutter  appear, 
sensitivity  is  reduced  by  means  of  the  PIN  diode  attenuator,  known 
locally  as  "the  bird  switch. 'I  Slight  clutter  reduction  occurs;  the 
loss  of  sensitivity  on  aircraft  is  objectionable  FTC,  IAGC  and  CSS 
(Cross  Section  Sensitivity,  similar  to  STC)  are  available  but  are 
seldom  used  on  the  ground  that  they  offer  no  significant  performance 
improvement  in  clutter. 


4V  Radar  Clutter  Reports 

AE,  "anomalous  propagation",  was  found,  to  be  loosely  used 
to  refer  to  many  observed  radar  returns  not  identifiable  as  aircraft, 
surface  vehicles,  weather,  or  normal  ground  clutter.  It  may  include 
returns  from  birds,  insects,  atmospheric  inhomoganeities,  normally- 
undetected  ground  features,  etc.  In  information-gathering  discussions 
it  proved  important  to  establish  a  common  working  definition  of  the 
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term  and  then  to  distinguish  clearly  between  observations  of  AP  and 
interpretations.  An  effort  is  made  to  preserve  this  distinction  in 
the  following.  For  present  purposes,  birds  are  considered  separately 
and  AP  is  defined  as  follows: 

AP:  Any  observed  radar  returns  which  cannot  be 
associated  with  aircraft ,  surface  vehicles , 
birds,  weather  or  normally-detected  ground 
clutter. 

The  following  is  a  distillation  of  conversations  with 
approximately  15  FAA  personnel  at  Little  Rock  including  controllers, 
radar  technicians,  supervisors,  and  administrators.  A  high  degree  of 
consistency  was  noted  among  the  reports. 

Bird  clutter  at  Little  Rock  Imposes  operational  limitations 
on  traffic  flow.  Aircraft  not  equipped  with  beacons  cannot  be  tracked 
through  sectors  in  which  bird  clutter  appears.  Beacon-equipped  air¬ 
craft  must  be  skin-painted  to  initiate  tracking;  this  is  Impossible  in 
dense  bird  clutter.  When  bird  clutter  is  strong,  controllers  advise 
aircraft  of  the  lack  of  radar  coverage  in  certain  areas..  A  NOTAM  is 
normally  not  Issued. 

Bird  clutter  is  observed  throughout  October,  November, 
February  and  March.  It  occurs  twice  daily,  at  dawn  and  at  dusk,  for 
a  period  of  30  to  90  minutes.  Radar  clutter  is  correlated  in  location 
and  time  of  appearance  with  dense  streams  oL  birds  observed  visually 
from  the  tower. 

The  appearance  of  bird  clutter  on  a  PPI  is  shown  in 
Figures  1  through  4.  The  clutter  in  these  photographs,  taken  on 
March  8  and  9,  was  considered  typical  by  Little  Rock  personnel. 

The  characteristics  of  the  clutter  are  discussed  in  Section  5. 

AP  which  is  thought  to  result  from  causes  other  than  birds, 
is  usually  observed  to  the  northeast  and  southwest;  sometimes,  it  is 
seen  to  approach  from  the  west.  It  appears  on  a  PPI  as  a  moving 
swarm  of  discrete  point  echoes  resembling  aircraft  returns  (Figure  5) . 
The  AP  occurs  frequently;  we  were  assured  thatv  if  we  stayed  for  two 
days,  we  would  definitely  see  it.  It  occurs  more  often  at  night  than 
during  the  day.  The  AP  often  occurs  in  a  visually  clear  atmosphere 
and  cannot  be  associated  with  any  visible  clouds  or  precipitation. 

MT1  is  useful  in  reducing  the  clutter,  but  does  not  eliminate  it.  A 
figure-8  distribution,  as  seen  in  Figure  5,  is  typical,  although 
occasionally  the  clutter  appears  at  all  azimuths.  CSS  (Cross  Section 
Sensitivity)  is  of  some  use  in  attenuating  the  clutter;  it  attenuates 
aircraft  returns  as  well.  FTC  and  IAGC  are  of  no  use. 
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At  both  Little  Rock  end  Memphis,  returns  like  those  just 
described  vere  frequently  ettributed  to  ’'ducting,"  Other  character¬ 
istic  effects  ascribed  to  ducting  vere  a  reduction  in  aircraft  detec¬ 
tion  range  and  an  increased  intensity  and  extent  of  ground  clutter 
return.  Xnfonaatlqn  on  the  actual  presence  of  ducts  has  not  been 
sufficient  to  permit  a  correlation  to  be  established  between  ducts 
and  radar  effects.  The  extended  detection  ranges  corotonly  attained 
with  marine  radars  in  ducts  vere  never  noticed. 


5.  Observations 


Radar  clutter  was  observed  and  photographed  on  a  PP1  on 
the  evening  of  March  8  and  the  morning  of  March  9.  Table  1  is  a  log 
of  all  photographs  taken.  Figures  1  through  4  are  photographs 
selected  to  show  the  major  features  of  observed  clutter. 

Ground  clutter  (Figure  la)  extends  in  range  to  between  7 
and  15  miles  in  different  sectors.  It  renders  NORMAL  video  useless 
inside  10; miles.  Usually,  MTI  video  with  STC  is  used  inside  15  miles 
and  NORMAL  is  used  between  15  and  40  miles,  the  customary  maximum 
surveillance  range.  No  residual  ground  clutter  was  evident  with  MTI. 

In  the  evening,  March  8,  clutter  was  observed  to  appear 
rather  abruptly,  within  several  minutes,  on  a  scope  which  had  pre¬ 
viously  shown  a  practically  black  background.  It  lasted  without 
major  variation  from  15  minutes  before  sunset  until  about  45  minutes 
after  sunset,  1800  to  1900  CST.  The  occurrence  of  the  evening  air 
traffic  peak  in  this  Interval  increases  the  impact  of  the  clutter  on 
controllers.  The  clutter  was  most  dense  to  the  norths  northeast  and 
east,  where  bird  roosts  are  known.  This  specific  geographic  relation 
of  airport  and  roosts  aggravates  the  problem.  Because  the  two 
principal  runways  run  north-south  and  northeast-southwest,  either 
departure  or  approach  paths  pass  directly  through  the  bird  clutter 
region  on  controllers'  PPI's. 

Morning  clutter  began  in  the  east  about  40  minutes  before 
sunrise,  0628  CST.  It  faded  almost  to  the  point  of  disappearance 
after  about  30  minutes.  Clutter  then  reappeared,  most  strongly  in 
the  north  and  south  quadrants,  and  remained  for  over  an  hour.  It 
persisted  on  into  the  beginning  of  the  morning  air  traffic  peak.  At 
times  the  clutter  extends  to  ranges  of  approximately  12  miles 
(Figure  4). 
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Although  the  gross  clutter  patterns  differed  somewhat 
between  evening  and  morning  (Figures  2  and  3) ,  the  fine  detail 
visible  on  the  PPI  looked  identical  on  both  occasions.  The  video 
consisted  of  clouds  of  discrete,  sharply-defined  returns.  The  PPI 
blips  appeared  to  be  equal  to  or  smaller  in  size  than  those  from 
a  small  aircraft.  No  diffuse  or  extended  returns  were  noted,  (The 
evident  merging  of  retu:.*ns  in  the  photographs  is  a  photographic 
effect.)  Some  blips  could  be  singled  out  and  could  apparently  be 
tracked  by  eye  for  several  scans.  The  echoes  to  the  south  in 
Figure  4b  were  tracked  in  this  way  for  more  than  a  minute;  each 
appeared'  co  maintain  an  identity  and  a  fairly  steady  southerly 
course  for  that  long.  Eventually,  the  echoes  faded.  Other  echoes, 
strong  and  dear  on  one  scan,  would  fall  to  reappear  in  the  same 
neighborhood  during  the  next  three  or  four  scans. 

Visual  observations  were  made  at  about  the  time  Figure  4a 
was  taken.  A  serpentine  stream  of  blackbirds  was  seen  flying  south 
past  the  airport  from  the  vicinity  of  a  largi  roost  known  to  lie 
several  miles  to  the  north.  The  stream  was  estimated  as  about 
100  -  300  feet  wide  and  the  same  in  thickness.  It  extended  to  the 
limits  of  visibility  to  the  north  and  to  the  south  at  an  estimated 
altitude  of  500  -  1000  ft.  Similar  large  concentrations ‘could  be 
detected  at  greater  ranges.  A  few  flocks  of  20  -  50  birds  were 
visible,  also  heading  southward.  The  locations  and  extent  of  the 
bird  concentrations  correlated  plausibly  with  observed  PPI  clutter 
when  the  effect  of  Mil  was  included.  As  figures  show.  Mil  cancel¬ 
lation  was  more  effective  toward  the  west,  where  the  birds’  southerly 
track  was  roughly  tangential,  than  to  the  north  and  south,  where  the 
motion  of  the  birds  was  approximately  radial  with  respect  to  the 
radar. 


6.  Bird  Information 

Mr,  and  Mrs i  Henry- Halberg,  Thurman  Booth  and  Robert 
Pierce,  experienced  Little  Rock  bird  observers,  provided  much  useful 
consistent  information,  pertinent  parts  of  which  are  summarized 
below. 


Several  million  blackbirds ,  grackle*  and  starlings  roost 
(sleep)  in  railroad  yards  several  miles  north  of  the  airport,  and 
about  one  million  roost  near  the  river  to  the  east  of  the  airport, 
locations  corresponding  to  the  regions  of  dense  clutter  shown  in 
Figures  2-4.  The  birds  disperse  to  feed  after  first  light  and  then 
return  to  the  roost  in  a  more  concentrated,  briefer  flew  at  sunset, 
as  observed  on  radar.  In  previous  years  bird  populations  have  been 
as  much  as  tan  times  greater.  The  birds  roost  around  Little  Rock  in 
large  quantities  (2  -  50  million)  from  about  October  through  March. 
Usually,  Little  Rock  roosts  begin  to  break  up  in  mid-March  and  the 
large  concentrations  have  disappeared  by  about  April  1.  It  was  re¬ 
ported  as  very  probable  that  those  birds  causing  the  ASR  clutter 
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would  be  gone  before  April  1.  This  year,  migration  may  be  occurring 
about  a  week  earlier  than  usual.  Resident  birds,  in  the  area  through¬ 
out  the  sumner,  appear  in  small  numbers  and  small  concentrations. 

After  blackbirds  begin  migrating,  temporary  roosts  form 
tr  the  north,  lasting  from  several  days  to  several  weeks.  Day-to-day 
bird  motion  is  strongly  affected  by  local  weather;  on  a  weekly  time 
scale,  the  long-term  weather,  e.g.,  "a  warm  winter,"  is  the  major 
influence  on  when  bird  migrate. 


7.  Little  Rock  Facilities 


The  ASR  is  sited  on  airport  property  in  an  FAA-controlled , 
fenced  compound  of  adequate  size  to  contain  three  vans  within  approxi¬ 
mately  100'  of  each  other  and  within  100*  of  the  ASR  building.  The 
ASR  site  is  within  a  mile  of  the  Tower/IFR  room.  Clear  radar  surveil¬ 
lance  is  possible  on  all  azimuths  except  through  the  tci/er. 

Surrounding  terrain  is  very  flat  with  the  exception  of  low  hills 
about  10  miles  to  the  south.  The  ASR  site  is  accessible  from  public 
roads . 


Adequate  prime  power  (30  KVA  plus)  is  available  from 
main  feed  lines  at  the  ASR  site.  Local  loads  are  automatically  switched 
to  a  generator  if  the  normal  supply  fails.  In  the  event  of  a  power 
failure  when  APL  test  equipment  was  on  line,  the  APL  load  would  not 
be  switched  to  the  FAA  generator  but  would  be  dropped. 

We  indicated  our  Intention  to  provide  the  following 
protection  at  Interfaces  with  the  local  system: 

a)  Fused  prime  power  lines  to  APL  equipment 

b)  Master  cutoff  switch(es)  in  APL  prime 
power  lines 

c)  Isolation  amplifiers  m  signal  lines  (as 
at  Friendship) 

An  intercom  on  the  ASR  tower  can  be  relocated  to  provide 
communication  between  the  APL  van  and  the  controller/ technician  loca¬ 
tion  in  the  tower  building. 


8,  Memphis  Angel  Clutter  Situation 

For  the  last  three  years  Memphis  International  Airport 
has  had  no  operational  interference  from  birds .  Other  "anomalous 
propagation"  returns  cause  no  significant  problems  in  air  traffic 
control.  Informal  discussions  on  March  10  with  Memphis  Tower  Chief 
James  Arthur,  Deputy  Tower  Chief  Parsons  and  three  watch-standing 
traffic  control  supervisors  are  summarized  as  follows. 
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Until  three  years  ago  Memphis  ATC  observed  on  radar 
heavy  concentrations  of  birds  during  fall,  winter  and  spring. 
Returns  occasionally  extended  as  far  as  40  miles.  Traffic  control 
was  affected  because  primary  radar  coverage  was  impossible  in 
some  approach  and  departure  sectors.  Three  years  ago  the  city 
dump  in  which  the  birds  roosted,  located  near  the  edge  of  the 
airport,  was  eliminated  and  the  bird  problem  disappeared.  No 
modifications  to  the  radar  were  made  at  or  around  this  time. 

Since  then  birds  have  been  no  problem.  Occasional  groups  of 
migrating  birds  are  observed  in  transit  through  the  area;  these 
have  no  effect  on  operations .  Occasional  AP  in  the  form  of  moving 
clouds  of  discrete  echoes  is  observed,  but  it  is  not  considered 
an  operationally  significant  problem. 

The  Memphis  TRACON,  located  in  the  tower,  uses  an 
ASR-5  sited  on  the  opposite  side  of  the  airport.  Memphis  ARTOC 
is  located  near  the  ASR;  an  ARSR  radar  is  about  15  miles  sway. 

The  Memphis  TRACON  handles  about  800  to  1000  operations  on  an 
average  day.  As  a  rough  indication  of  radar  operation,  radar 
handoffs  between  terminal  and;  enroute  facilities  were  claimed  to 
be  normally  accomplished  without  difficulty  in  the  range  interval 
30  to  40  miles  from  the  ASR,  for  small  aircraft  as  well  as  for 
large. 


A.  0540CST 

MTI  AND  STC  WITH  NORMAL  BACKGROUND.  TWO  SCANS, 
10  MILE  SCALE,  2  MILE  RANGE  RINGS.  MAPCH  9. 

BIRD  CLUTTER  STARTING  TO  GROW  AT  085°,  4.5  Ml. 


B.  0604CST 

MTI  AND  STC.  ONE  SCAN,  2  MILE  RANGE  RINGS. 
MARCH  9.  NO  AIRCRAFT  KNOWN  TO  BE  IN  AREA. 

FIGURE  1 
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A.  1810 —SCALE  10  MILES, 
2  MILE  RANGE  RINGS 


B;  1815-SCALE  10  MILES, 
2  MILE  RANGE  RINGS 


C.  1820— SCALE  6  MILES  , 
2  MILE  RANGE  RINGS 


D.  1825— SCALE  20  MILES, 
5  MILE  RANGE  RINGS 


FIGURE  2 

DUSK  CLUTTER,  MARCH  8.  MTI  &  STC.  TWO  SCANS 
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A.  0600 


B.  0603 


C.  0609 


D.  0613 


FIGURE  3 

DAWN  CLUTTER,  MARCH  9,  Mil  &  STC 
TWO  SCANS,  10  MILE  SCALE 
2  MILE  RANGE  RINGS 


magnetic  month  towand  top  of  PAGE 
TIME* _ CST 


A.  0635 _ MTI  &  STC.  TWO  SCANS,  20  MILE  SCALE, 

5  MILE  RANGE  RINGS.  MAR  9 


B.  0702 _ MTI  ft  STC.  TWO  SCANS,  6  MILE  SCALE, 

2  MILE  RANGE  RINGS.  MAR  9 


FIGURE  4 
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B— 2320CST 

FIGURE  5 

ANGEL  CLUTTER  ON  LITTLE  ROCK  ASR 


(30  milk  ranw.  courtesy  or 
C.  ORAHAM  .  FAA,  LITTLE  ROCK  ) 
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TABLE  1 

Little  Rock  ASR  4B  PPI  Photographs 
March  8-9.  1972 


Unices  otherwise  Indicated,  photos  show  MTI  video,  STC  on,  parametric 
amplifier  on,  10  mile  scale,  2  mile  range  rings  (RR),  two  sweeps /photo 

Polaroid  camera  with  portrait  lens,  f/5.6,  exposure  8  seconds  (2 
sweeps)  unless  otherwise  indicated,  ASA  3000  film,  21.5"  from  portrait 
lens  to  PPI  plastic  overlay  (approx.  22,5"  to  CRT  face). 


No. 

Time  (CST) 

Remarks 

2 

1805 

3 

10 

- 

11 

Sunset 

4 

12 

5 

15 

Parametric  amplifier  off 

6 

15+ 

11 

20 

6  ml  scale,  2  mi  RR 

12 

22 

6  mi  scale,  2  mi  RR,  1  sweep 

13 

25 

20  mi  scale,  5  ml  RR 

21 

0540 

MTI  with  NORMAL  background  (low 
NORMAL  video  gain) 

22 

48 

MTI  with  NORMAL  background,  air¬ 
craft,  birds  090/10 

23 

50 

24 

52 

(wind  reported  250/4  KTS) 

25 

53 

26 

56 

Some  birds  beginning  to  appear 
at  'v  010 

27 

0600 

28 

03 

Birds  to  east  fading 

29 

04 

1  sweep 

30 

08 

Birds  appearing  to  north 

31 

08+ 

20  mi  scale,  5  mi  RR 

32 

09 

33 

13 
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\  Table  1  Continued 


No. 

Time  (CST) 

Remarks 

34 

15 

40  ml  scale,  10  ml  RR,  4  sweeps, 
surface  vehicles  240,  9  ml, 
aircraft  NE  35  mi 

35 

22 

MTI  and  NORMAL 

36 

23 

MTI  and  NORMAL 

37 

0624 

6  ml  scale,  2  mi  RR,  MTI  and  NORMAL 

38 

26 

6  ml  scale,  2  mi  RR, 

- 

28 

Sunrise 

39 

32 

20  ml  scale,  5  mi  RR 

40 

35 

20  mi  scale,  5  mi  RR;  200-600 
discrete  blips  to  north 

37-45 

Visual  observations  of  bird  streams 
(see  text) 

41 

48 

20  mi  scale,  5  ml  RR,  5  sweeps 

42 

0702 

6  mi  scale,  2  ml  RR;  clumps  to  south 
are  moving  south;  runway  MTI 

beacons  NE  and  SW  approx.  1  ml 
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Enclosure  (1) 

Principal  Contacts  In  Little  Rock 


FAA.  Little  Rock 

FAA  Sector  Chief  Noteboom 
Asst.  Sector  Chief  Daniel  Gardner 
Radar  Unit  Chief  Everett  Graham 
Tower  Chief  Uunsekker 


Airway  Facilities  Sector 
Building  C-8,  Adams1  Field 
Little  Rock,  Arkansas  72202 
Telephone  No.  501-374-9038 


Little  Rock  Avifuana  Contacts 
Thurman  Booth 

Federal  Blackbird  Control  Officer 
U.S.  Fish  and  Wildlife  Service 
Room  315,  Post  Office  Building 
Little  Rock,  Ark. 

Telephone  No.  501-378-5382 


Robert  Pierce 

Cooperative  Extension  Service 

Unlv.  of  Arkansas  Graduate  Institute  of  Technology 
MacArtbur  Park,  Little  Rock,  Ark. 

Telephone  No,  501-376-6301 
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II.  REPORT  OF  TRIP  TO  MINNEAPOLIS /MILWAUKEE  AIRPORTS 

March  1972 


The  airports  and  FAA  facilities  at  Minneapolis ,  Minnesota 
and  Milwaukee,  Wisconsin,  were  visited  during  the  period  15  -  17  March 
1972  by  APL  personnel.  The  primary  purpose  of  these  visits  was  to 
determine  the  feasibility  of  conducting  FAA  Task  II  test  operations 
at  either  of  these  locations. 

Upon  arrival  at  Minneapolis  a  tour  of  the  Tower,  IFF. 
room  and  maintenance  areas  was  provided  by  the  Assistant  Radar  Chief. 
Discussions  with  personnel  in  these  areas  revealed  a  history  of  un¬ 
defined  angel  clutter,  primarily  in  the  spring  and  fall.  It  was 
expected  that  the  spring  clutter  problems  would  begin  about  1  April 
and  continue  for  approximately  two  months.  Controllers  interviewed 
could  not  recall  any  occasions  when  clutter  had  become  an  operational 
problem,  although  it  was  an  annoyance.  Maintenance  personnel  were 
cordial  and  receptive  to  the  possibility  of  locating  Task  II  equipment 
at  Minneapolis. 

Conversations  with  several  local  "bird  watchers"  were 
consistent,  with  the  arrival  of  waterfowl  anticipated  about  1  April 
followed  closely  by  other  types ,  with  peak  activity  in  mid  to  late 
April. 


Upon  arrival  at  Milwaukee,  the  Task  II  requirements  were 
discussed  with  the  Assistant  Sector  Chief  and  a  tour  of  the  Tower  and 
IFR  room  was  made  in  company  with  the  evening  maintenance  technician. 
All  personnel  contacted  agreed  that  heavy  angel  clutter  could  be  anti¬ 
cipated  during  the  spring  and  fall  of  each  year  with  the  spring  clutter 
historically  beginning  about  1  April.  In  contrast  with  Minneapolis, 
however,  the  Milwaukee  air  controllers  interviewed  stated  that  clutter 
has  been  an  operational  problem.  The  source  of  the  heavy  clutter  at 
Milwaukee  is  undefined  after  several  attempts  to  correlate  with  atmos¬ 
pheric  conditions.  Photographs  of  previous  clutter  observations  dis¬ 
closed  figure  eight  patterns  centered  on  the  radar  and  extending,  ip 
some  cases,  to  twenty  miles.  The  one  point  of  correlation  noted  was 
the  tendency  of  the  clutter  pattern  to  align  with  the  wind.  Peak 
clutter  periods  tend  to  occur  around  sunrise  and  sunset,  however, 
clutter  has  been  observed  at  all  hours. 
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All  facilities  at  Milwaukee  are  centrally  located 
providing  eaay  access  to  the  ASR  (located  on  the  field  In  front 
of  the  control  tower),  the  tower  cab,  IFR  room,  FAA  offices. 
Weather  Bureau  and  Flight  Service  Station.  The  Tower  Chief 
and  Sector  Chief  each  expressed  high  interest  in  the  project 
and  their  willingness  to  assist. 

The  Milwaukee  "bird  watchers"  contacted  reported 
soma  migratory  waterfowl  activity  beginning  now,  with  overall 
bird  activity  expected  to  reach  a  peak  in  April. 

In  view  of  /the  apparently  heavier  clutter  conditions 
and  resultant  operational  problems,  anticipated  bird  activity 
and  reduced  travel  distance,  it  is  recommended  that  Milwaukee  Air 
port  be  selected  as  the  FAA  Task  11  test  site  with  an  equipment 
arrival  date  planned  for  29  March  1972. 
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STATISTICAL  DATA  FOR  BIRDS  AND  AIRCRAFT 

Spatial  Autocorrelation  Function!  for  ASR  Video  Matrix 

The  characteristic  distributions  in  azimuth  and  range 
bf  ASR  video  returns  are  of  importance  for  the  optimization  of 
a  video  processing  scheme.  They  determine,  respectively,  the 
azimuthal  extent  over  which  pulse-to-pulse  correlation  is  profitable 
and  the  sampling  interval  in  range  needed  to  yield  independent 
video  samples.  A  measure  of  these  characteristics  is  the  one-dimen¬ 
sional  autocorrelation  function  presented  below.  Initial  results 
include  range  and  azimuth  functions  for  a  known,  controlled  Cessna  172 
aircraft  and  for  an  angel  clutter  target.  (Figure  B-l) 

The  plotted  functions  ate  autocorrelation  functions  as 
defined  by  Bendat  and  Piersol  (Reference  *  ),  normalized  by  the 
value  of  the  function  at  zero  lag.  The  input  data  from  each  of 
several  scans  consisted  of  a  string  of  ASR  amplitude  values 
(37  in  azimuth  or  40  in  range)  from  a  selected  row  or  column  in 
the  Data  Collect  Matrix.  The  autocorrelation  function  of  each 
string  was  found  and  normalized.  Then  the  average  autocorrelation 
function  over  the  several  scans  was  obtained.  Figure  B-2  shows 
Ear;’*-  autocorrelation  functiors  for  a  Cessna  172  aircraft  and  for 
an  angel  track  (MKE  16-1,  Track  I)..  Figure  B-3  shows  azimuth 
autocorrelation  functions  for  the  same  two  targets. 

The  range  autocorrelation  results  show  that  both  the 
aircraft  and  angel  returns  usually/  fall  into  a  single  range  cell 
but  occasionally  appear  in  two  range  cells,  as  is  expected  from 
straddling  of  two  range  cells.  The  azimuthal  autocorrelation 
functions  show  a  slight  difference  for  the  aircraft  and  the  angel, 
with  the  angel  track  having  slightly  higher  autocorrelation  for 
lags  less  than  0,45°  and  slightly  lower  thereafter.  Both  targets 


Bendat,  J.  S.  and  A.  G.  Piersol,  Random  Data:  Analysis  and 
Measurement  Procedures,  Wlley-Interscience,  New7 York,  1971 
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Tracks  Uaad  For  Statistical  Analysis 


CESSNA  172  (MKE  21-1)  194  SCANS 
BIRD  ANSEL  CLUTTER  (MKE  16-1 )  90  SCANS 
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Figure  b-2  ASR  Mean  Range  Autocorrelation  Coefficient  For  A  Lag  Interval1  Of  One 

Pulaevidth 


Figure  B-3  ASR  Mean  Azinuth  Autocorrelation  Coefficient  for  a  Lag  Interval 
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reach  a  value  of  0,5  at  2  lags  (0.3°),  which  is  close  to  the 
value  of  2-3  lags  expected  due  to  the  antenna  pattern  itself 
on  a  point  target.  These  results  confirm  single-scan  observations, 
i.e.,  if  targets  are  smaller  than  the  radar  resolution  cell  (one 
pulsewidth  by  one  beamwidth)  then  the  run-length  similarity  for 
the  Cessna  and  the  angel  track  persists  over  many  scans. 

The  results  are  similar  for  both  the  aircraft  and  the 
angel.  From  this,  one  may  conclude  that  angel  clutter  (presumably 
birds)  resembles  general  aviation  aircraft  in  their  mean  spatial 
distribution.  Analysis  of  additional  aircraft  and  angel  data  is 
required  to  determine  whether  the  angel-aircraft  differences  in 
Figures  B-2  and  B-3  are  significant. 

Amplitude  Distributions  for  ASR  Video  Matrix 

When  the  MPS-19  radar  tracked  an  aircraft  or  angel,  the 
ASR  Data  Collect  Matrix  was  automatically  positioned  in  the  area 
containing  the  tracked  target.  ASR  video  amplitude  in  each  of 
the  40  x  37  matrix  cells  was  sampled  each  scan,  while  both  radars 
observed  the  same  target  and  its  surroundings.  Approximate 
probability  densities  in  the  form  of  relative  amplitude  frequency- 
of-occurrence  historgrams  and  the  corresponding  cumulative 
distributions  have  been  plotted:. 

Frequency-of-occurrence  as  a  function  of  relative  ASR 
video  amplitude  is  shown  in  Figure  tM  for  a  Cessna  172  and  in 
Figure  B-5  for  two  angel  tracks..  Abscissas  in  the  plots  are  the 
logarithm  of  the  range-normalized  video  amplitude,  discussed  below; 
ordinates  are  the  number  of  occurrences  in  a  target  track. 

The  amplitude  distributions  were  obtained  by  assuming 

that  the  distribution  type  did  not  change  with  target  range.  Then, 

_2 

the  decrease  in  amplitude  with  increasing  range,  e.g.  as  R  for 
a  point  target,  was  compensated  for  by  multiplying  target  amplitudes 
in  the  matrix  from  each  scan  by  the  square  of  the  range  measured 
at  the  time  the  matrix  data  were  collected.  Accurate  range  was 
obtained  from  the  continuous  target  track  maintained  by  the  MPS-19. 
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Figure  B-4  AS*  Amplitude  Frequency-Of-Occurrence  for  Cessna  172 
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In  Initial  data  analysis,  range  factors  of  R  ,  R  ,  and  R 

_2 

were  used,  and  in  all  cases,  the  density  obtained  with  the  R 

factor  had  the  smallest  half-width.  This  was  interpreted  as 

an  indication  that,  amo;'<*  the  three  range  dependences  tried, 

—2 

R  most  closely  described  the  target  amplitude  variation* 

since  compensation  with  an  inapplicable  range  factor  broadens 

_2 

the  density.  The  R  amplitude  dependence  (which  corresponds  to 
-4 

an  R  dependence  of  signal  power)  implies  that  the  angel  targets 
behaved  approximately  as  point  targets,  like  aircraft,  which  is 
in  agreement  with  the  conclusion  from  single  scan  analysis  above, 
that  the  angels  were  not  comparable  to  or  larger  in  size  than  an 
ASR  resolution  cell. 

In  addition  to  amplitude  variations  caused  by  target 
cross-section  fluctuation  and  by  target  range  change,  the  video 
amplitude  varies  due  to  the  gain  change  in  the  direction  of  the 
target  as  the  antenna  rotates,.  The  result  is  a  smearing  of  the 
target  amplitude  density  toward  lower  amplitudes.  To  deduce  target 
characteristics  from  the  observed  distribution,  the  known  antenna 
contribution  must  be  removed.  In  the  present  problem,  it  is  the 
ASR  video  characteristics,  including  the  antenna  effect,  which  are 
of  interest.  These  are  shown  in  the  figures. 

At  low  amplitudes  (0,  1,  2,  3,  on  a  0-31  scale),  the 
distributions  of  target  amplitude  and  of  Data  Collect  Matrix 
background,  which  consists  of  receiver  noise  and  MTI  residue, 
overlap  and  initially  cannot  be  separated.  Therefore,  the 
histograms  shown  were  formed  from  all  740  independent  samples  in 
a  matrix.  In  these  histograms,  over  90%  of  the  total  number  of 
samples  lie  in  the  first  two  intervals.  These  are  predominantly 
matrix  background  amplitudes;  the  typical  mean  value  of  background 
was  between  0.7  and  1.0.  Occasional  returns  from  targets  other 
than  the  one  being  tracked  may  be  received;  these  targets  pass 
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through  the  Data  Collect  Matrix  within  a  few  scans  and 
their  contributions  to  the  resulting  histograms,  which  are 
accumulated  over  many  scans,  are  negligible.  Amplitudes  greater 
than  approximately  three  are  therefore  assumed  to  be  due  solely 
to  the  tracked  target ,  with  high  probability . 

From  the  amplitude  frequency-of-occurrence  data, 
cumulative  probability  distributions  were  formed  for  the  Cessna 
172  aircraft  and  for  angels.  To  reduce  the  strong  effect  of 
background,  the  densities  ab've  were  truncated  or,  equivalently, 
an  amplitude  threshold  was  established.  Thresholds  of  0,  1,  2,  3 
and  4  on  a  scale  of  0-31  were  used.  The  resulting  distributions 
are  shown  on  Weibull  distribution  plots  in  Figures  B-6  through  B-8 
for  angel  targets  and  Figures  S-9  and  B-10  for  the  Cessna. 

Similar  normal  probability  plots  are  shown  in  Figures  B-ll  for 
the  Cessna  and  Figure  B-12  for  an  angel  target. 

The  Weibull  distribution  is  a  two-parameter  distribution 
with  density  given  by 

p(x)  -  lx)*"1  e~(£ 

X  1 

and  cumulative  distribution 


The  Weibull  family  includes  the  exponential  (g  *  1)  and  the 
Rayleigh  distribution  (B  “  2).  The  distribution  is  useful  in 
searching  for  an  analytical  description  of  a  measured  distribution 
or  in  comparing  measured  to  theoretical  distributions.  A  straight 
line  on  the  Weibull  plot  indicates  exact  agreement  with  the 
Weibull  distribution;  the  slope  of  the  line  can  be  used  to  determine 
the  value  of  the  Weibull  parameter,  £. 
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THRESHOLD  *0 


Figure  B-6  WEIBULL  PLOT  OF  MATRIX  AMPLITUDES  -  ANGEL  TRACK  15-2 


WE I BULL  PLOT  OF  MATRIX  AMPLITUDES  -  ANGEL  TRACK  16-3 


Figure  B-ll  LOG  NORMAL  PLOT  OF  MATRIX  AMPLITUDES  -  CESSNA  TRACK  22-1 
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Cumulative  Probability 
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The  Weibull  plots  of  aircraft  and  angel  distributions 
(Figures  B-6  through  B-10)  show  similar  features.  With  threshold 
values  of  0  and  1,  the  compound  nature  of  the  distribution  is 
evident.  Background  is  included  and  dominates  the  distribution. 

When  the  threshold  is  raised  to  three  or  four,  on  the  other  hand, 
several  differences  appear.  For  a  given  target,  the  slopes  of 
the  curves  are  quite  similar  for  threshold  amplitude  aT  *  3  and  ■  4. 
Furthermore,  the  values  of  the  slopes,  giving  the  parameter  g  in 
the  Weibull  distribution,  lie  between  0.9  and  1.4  for  the  aircraft 
and  between  1.1  and  1.7  for  the  angels;  actual  values  are  given  in 
the  figures .  These  results  indicate  that  the  observed1  distributions 
lie  between  the  exponential  and  the  Rayleigh  distributions  in  type. 
Small  amplitudes  are  probable  but  smallest  amplitudes  (0,  1,  and  2) 
less  likely,  as  shown  by  the  small  change  when  the  threshold  is 
raised  from  3  to  4.  (For  the  exponential  distribution,  smallest 
amplitudes  afc  most  probable;  for  the  Rayleigh  type,  smallest 
amplitudes  are  least  probable  of  all  amplitudes  below  the  mean) . 

The  plot 8  of  cumulative  distribution  on  Log  Normal 
probability  paper.  Figures  B-ll  and  B-12,  are  similar  to  the  above 
in  variation  as  a  function  of  threshold.  Again,  curves  for  aT  *  3 
and  4  lie  close  with  similar  slopes.  The  estimated  goodness-of-fit 
to  a  straight  line  is  not  markedly  better  for  either  Weibull  or 
normal  plots.  For  a  threshold  of  3,  the  Log  Normal  parameter  p 
(mean-to-median  ratio)  was  about  10  for  the  aircraft  and  3  to  6.2 
for  the  angel  track;  for  a  threshold  of  4,  the  values  were  6.2  and 
2.3  respectively. 

Figures  B-13  and  B-14  show  a  typical  set  of  angel  track 
video  amplitudes  recorded  by  the  Track  Radar  Module.  Each  track 
radar  pulse  is  converted  to  digital  format  and  transmitted  to  the 
Data  Acquisition  Module,  where  it  is  written  on  digital  magnetic 
tape  along  with  range,  bearing,  and  elevation  measurements.* 

The  ordinate  in  the  figures  is  the  A/D  converter  level  recorded 
for  each  pulse.  The  sample  rate  was  the  track  radar  PRF  of  half 
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the  ASR  PRF.  Figures  B-13  and  B-14  thus  cover  0.243  and  0.208 
seconds,  respectively.  Both  figures  were  made  of  the  same  angel 
track  and  were  separated  in  time  by  four  ASR  scans  or  16  seconds. 

The  track  radar  was  in  the  automatic  range  track  mode  but  manual 
angle  track  was  used  to  avoid  the  30  Hz  conical-scan  modulations 
of  the  MPS-19  angle  tracking  function.  Automatic  Gain  Control 
(AGC)  was  used,  maintaining  the  signal  level  approximately  constant. 

Fourier  Analysis  of  LdPS-19  Video  Fluctuations 

To  observe  the  frequency  spectra  of  the  target 
fluctuations,  the  video  from  the  MPS-19  tracking  radar,  which  was 
used  to  place  a  video  gate  in  the  ASR  video  during  the  generation 
of  Data  Collect  Matrices,  was  simultaneously  sampled  and  recorded 
for  several  angel  and  aircraft  targets.  These  samples,  collected 
every  radar  pulse  while  the  radar  was  locked  on  the  target,  were 
used  to  generate  Fourier  Spectra  of  the  target  fluctuations  as 
seen  by  the  MPS-19.  Since  the  characteristics  of  the  MPS-19  are 
almost  identical  to  those  of  the  ASR,  and  as  the  two  radars  were 
colocated  in  Milwaukee,  the  results  should  be  representative  of 
the  ASR  video  for  the  targetc.  For  each  spectra,  210  simples 
were  used  giving  a  resolution  of  approximately  0.6  Hertz-  for  a 
600  Hertz  sampling  rate.  The  frequency  distributions  were  averaged 
for  20  consecutive  sets  of  data  oh  each  target,  representing 
approximately  34  seconds  of  video  (or  -  20,000  samples) .  Figure  B-15 
shows  two  sample  frequency  distributions  for  a  bird  angel  and  the 
Cessna  172  which  result  from  this  procedure.  Note  that  the  ordinate 

for  the  aircraft  plot  covers  approximately  twice  the  range  of  the 

\ 

ordinate  for  the  angel  plot.  In  both  plott?  the  conical  scan 
frequency  (30.5  Hertz)  of  the  tracking  radat.  is  clearly  visible. 
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The  general  characteristics  of  the  plots  indicate  the  angels 
appear  to  be  represented  by  low  frequency  fluctuations,  with  a 
prominent  peak  about  10  Hertz  corresponding  to  the  "wing  beat" 
frequencies  found  by  other  investigators.  The  aircraft,  on  the  ether 
hand,  had  higher  frequency  modulations,  generally  characterised  by 
multiple  frequency  series  similar  to  those  found  in  the  normal 
Fourier  decomposition  of  periodic  nonsinusoidal  functions.  For 
example,  there  are  peaks  at  82  Hertz,  164  Hertz,  248  Hertz.  Eighty- 
two  Hertz  appears  to  be  the  fundamental  frequency  of  this  series, 
with  164,  248,  etc.,  being  the  first,  the  second,  etc.  overtones 
normally  found  in  the  Fourier  decomposition.  Frequency  multiples 
above  the  Nyqulst  frequency,  in  this  case  300  Hertz,  are  found  to 
be  "folded  over"  and  observed  at  lower  frequencies.  For  example, 
the  series  with  a  fundamental  frequency  of  80  Hertz  has  a  multiple 
equal  to  320  Hertz  which  is  higher  than  Nyqulst  frequency.  As  a 
result,  this  multiple  is  folded  over  and  observed  at  280  Hertz. 

Similar  effects  occur  for  higher  multiples.  Table  I  lists  the 
observed  and  folded  frequencies.  It  is  interesting  that  the  aircraft 
engine  used  for  this  run  nominally  operated  at  2400  rpm,  which 
with  the  two-bladed  propeller  corresponds  to  4800  Hertz  per  minute, 
or  80  Hertz  per  second,  the  observed  fundamental  frequency. 

The  apparent  dependence  of  the  major  structure  of  the 
video  frequency  spectrum  on  the  Cessna  propeller  rate  points  out  the 
fact  that  great  care  must  be  used  in  developing  angeT/aircraft 
discrimination  techniques,  in  order  to  avoid  developing  discriminators 
which  operate  on  characteristics  peculiar  to  certain  targets.  The 
most  practical  way  of  avoiding  this  is  to  collect  a  large  data  base 
of  different  types  of  aircraft  and  angels.  While  the  80  Hertz  lines 
seen  in  the  track  radar  video  spectrum  are  not  necessarily  so  prominent 
in  the  ASR  MTI  video,  they  should  be  closely  investigated  in  future 
data  collection  efforts. 


B-23 


i 


r 

y 

t 


APPENDIX  B-2 

PROCESSOR  SIMULATION  RESULTS 

Detailed  results  of  each  Azimuth  Pattern  Processor  are 
listed  in  Tables  I  -  X-  The  performance  curves  in  the  test  were 
derived  from  the  scan-weighted  averages  of  aircraft  and  angel 
detection  percentages . 

The  track  numbers  of  each  angel  and  aircraft  do  not 
explicitly  appear  in  the  tables,  However,  they  can  be  inferred 
from  the  number  of  scans  of  data  for  each  track;  this  number  is 
listed  in  the  left-most  column,  separated  by  a  vertical  line,  in 
each  table  The  correlation  between  track  numbers  and  number  of 
scans  is  as  follows  for  the  Milwaukee  data: 

Aircraft  Track  Number  of  Scans 


6-1 

58-59 

6-2 

97 

3-2 

73 

3-3 

82 

21-1 

256 

22-1 

160 

Antjel  Track 

Number  of 

10-1 

24 

10-6 

69-70 

11-3 

19 

y\~u 

5 

M_5 

20 

1  ’-6 

35 

11-7 

7 

16-1 

62 

16-3 

30 
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For  the  Baltimore  data,  the  correlation  is 


Aircraft  Track 
2F 

3-BC 

A/C 

4C 

4D 

2A 

2B 


Number  of  Scans 

120 

34-35 

308 

129-132 

163-164 

155-180 

211-218 


A  short  description  of  each  processor  simulation  is  given  below  for 
each  table. 


1.  Simple  Amplitude  Detector  (Milwaukee  data) 

Table  I  lists  the  results  of  processing  the  Milwaukee 
data  with  a  fixed  threshold;  an  aircraft  is  declared  present  if  at 
least  one  amplitude  exceeded  the  threshold.  The  value  of  each 
threshold  tested  is  listed  after  "THR,"  in  the  top  line  of  the  table. 

The  maximum  video  amplitude  was  31,  corresponding  to  5-bit  quantization. 

All  of  the  Milwaukee  data  consisted  of  37  amplitude  samples 
taken  on  every  other  radar  sweep  (pulse  period) .  The  range  cell 
processed  was  chosen  by  summing  all  amplitudes  in  each  range  cell 
for  all  37  azimuth  samples  and  choosing  the  range  cell  with  the 
largest  amplitude  sum. 
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2.  H/H  Azimuth  Correlator  (Milwaukee  data) 

Tables  II  and  III  list  the  results  of  subjecting  the 
data  to  a  fixed  threshold  and  requiring  that  M  consecutive  samples 
exceed  the  threshold  for  aircraft  detection.  Table  II  applies  for 
M/M  *  2/2  and  Tabic  III  applies  for  M/M  *  3/3  as  the  thresholds 
were  varied  from  4  to  21  out  of  the  31  possible  levels. 
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3.  Dual  Threshold  Processor  (Milwaukee  data) 

This  processor  uses  two  thresholds ,  one  low  (given  by 
"THR."  in  the  tables)  and  one  high  (set  "MR."  levels  above  the  first). 
The  number  of  ACI  (Azimuth  Correlation  Intervals)  is  counted  at  each 
threshold.  A  third  threshold  is  set  on  the  ACI  difference.  An 
aircraft  is  declared  if  the  ACI  difference  is  less  than  this  third 
value,  given  by  "MN."  in  the  tables. 

Table  IV  lists  the  results  for  combinations  of  the  three 
threshold  values  when  all  37  azimuth  samples  were  processed.  Table  V 
applies  when  only  10  of  the  37  samples  are  processed  at  any  one  time; 
at  least  one  aircraft  detection  was  required  as  the  10-sample  window 
was  swept  across  the  full  37-sample  azimuth  coverage. 


B-31 


APPENDIX  B-2 


X  SO  TO 


X  X  04 


X  (N  TO 


»\C» 


3N!N 


(N  ,  o  to 


M 

<n  x  to  a 

•€ 

H 

X 

(N  X  (N  M 

O 

X 

« 

nr(N  a  e* 

a  m 

*i  c 

O  ON'n  H 

X  H 

to  O 

U  M 

«  O  (N  fM  fH 

X  W 

H  a . 

»••  •  •  •  *4 

«  x  «  .  w 

a  xxx  tfl 

a  f-me  r  x 


r*P»TO©o«—OOX 

TO  TO  P»  TO  O  O'  O  O  O' 


ff  <S  ®  O  (ft  9  r  vO  o 
p*  to  vo  ®  to  r*  r»  ^  to 


tfi  to  to  ®  to  m  tfi  top* 


mxooooooo 

TOTOOOOOOOO 


*“  P*  TO  O  0:0  O  O  TO 
TOTOTOOOOOOTO 


S'^OOOtflOCCO 
TO  TO  TO  O  TO  O  TO  TO 


TO»-ffloinrsooTO 

r**TOTOTOTO®OoTO 


isaiNoifloif  bo 
TOr^tftTOTOTOP^TO® 


SflXOO'-MO'S 
in  X  TO  X  X  TO  TO  O'  TO 


OXTOOOOOOO 

OTOTOOOOOOO 


ri-OiOQrOOO 

TOTOTOOOTOOOO 


NOflOONrif,  O 
XP-XXTOXP'TOX 


or^'nointfitNon 

f''TOXXr*>®TOOX 


<N  tfl  TO  OONO^n 
TO  TO  TO  TO  TO  TO  TOTO 


TO  TO  TO  ©  TO  r?  O  CO  TO 
(N  TO  r-  (N  JT  TO  TO  TO 


O  ©  ©  ©  ©  f*  ©  ®  © 
O  O  O  TO  TO  TO  O  TO  O 


TO  ©:©  O  ©  ©  f  TO  © 
TO  O  ©  O  TO  O  P»  TO  TO 


r*  p»  TO  O  ©  ©  •—  TO  TO 
TO  TO  TO  TO  TO  CO  P*  p 


TO  TO  TO  TO  ©  TO  TO  <N  O 
TO  X*-  TOTO  TOTO 


TO 

O  TO  ©  O  TO  TO 

TO 

TO 

O  TO  O  O  TO  TO 

TO 

^  r»  r~ 

TO 

O  TO  X  CO  X  X 

O 

X 

O  X  TO  TO  TO  TO 

TO 

*— 

X 

X  TO  tf>  TO  TO  M 

TO 

TO 

1 

TO  X  TO  TO  TO  TO 

X 

TO 

» 

©  X  ©  o  o  o 

TO 

O  9*  O  ©  O  O 
r-  *»*»•*•*• 

TO 

TO 

o  p*  o  o  o  o 

TO 

TO 

I 

©  O'  O  ©  O  o 

^  r»  r-  ♦— 

TO 

TO 

X  X  O  C  X  X 

1- 

| 

90 

TO  X  O  O  TO  TO 

TO 

C/3 

*-  r— 

u 

X 

o 

1 

X  TO  O  O  CD  X 

P* 

O' 

^oooooo 

TO 

Q 

r- 

o 

TO  *“  ©  P»  VO  X 

as 

r* 

St 

C/3 

TO 

c*  x  ©  <yv  O'  O' 

TO 

£ 

X 

TO 

r-  p*  x  to  to  o 

X 

p 

TO 

» 

X  SO  TO  TO  TO  TO 

X 

a 

TO 

X  X  O  O  o  O 

9v 

© 

o 

TO 

TO  TO  O  O  O  O 

TO 

M 

*" 

a 

TO 

X  TO  O  O  X  TO 

X 

o» 

O'  O'  ©  0-9*  O' 

TO 

C/3 

t 

•— 

X 

1 

TO  tfl  O  P~  O  SC 

it 

ro 

TO 

SO  X  O  TO  TO  TO 

TO 

1 

to 

r* 

> 

M 

\C 

ui 

P*  TO  O  O  X  « 

M 

X 

1 

© 

P*  X  O  OTO  TO 

TO 

» 

*. 

r- 

►4 

t- 

5 

X 

to  m  x  x  x  x 

TO 

h 

TO 

1 

M 

X  tfl  TO  TO  TO  TO 

X 

•  • 

U 

to 

ee 

X 

w 

TO 

1 

M 

X  *”  X  tfi  r»  TO 

u 

r* 

© 

tfi 

O. 

TO  TO  TO  TO  TO  TO 

o 

X 

■< 

<*. 

03 

1 

X 

as 

w 

TO 

o 

p*  p.  o  o  to  o 

M 

it 

TO 

1 

M 

P-  o  o  to  o 

► 

TO 

*< 

H 

r*  r* 

< 

U 

a 

r* 

1 

W 

TO  O  ®  *“  © 

a 

TO 

N 

TO 

1 

X  X  vo  p-TO  O 

M 

X 

1 

M 

f-t 

X 

a 

u 

TO 

X  X  (N  OTO  To 

o 

X 

H 

M- 

X 

X 

H 

-  ik,  - 

■TO 

SO  tf)  TO  P*  tC  TO 

M 

r- 

w 

31 

-f- 

{£ 

x  p*  to  fM>x  o 

X 

X 

TO 

u 

tf»  TO  P*  X;tf>  X 

X 

(N 

* 

TO 

X 

<N  r- 

<< 

p* 

u 

M 

u 

w 

X 

w 

APPENDIX  B-2 


* 


®  ®  <*t  CP  <N  vO  ©  m  O  P-  IT  I  O 

*-  (N  W(£|<»9^>0 


®  SO  fN  CO  »■  O  O'C  9-N-fl  \C 

*“  «-  *-  a  a  a 


cr.  \£  T- 

OOOOOOOa© 

r- 

x  rs  r. 

CT-  U1  v£i  O  in  O  P-  S3  O 

r" 

<N  m  »C  m  af  irs  oc  ® 

<x  rvi  rvj 

ST<  N  O  O  O  ®  (N  O  vr 

*— 

•*  a  03  iC 

x  rvj  *- 

O  «-  O  C  O  o  o  -a  c 

cl  <n 

vcvon 

inmir.  ooshnci 

r» 

n  *-  a  in  in  m  oo  r*. 

'C  vC  OJ 

Ofooin»co>£n 

*■» 

fM  M'On 

vT  \0  *- 

oocoooo 

m 

vi  r\| 

croinooxr-no 

*■» 

<s  *-  a  ni  <n  m  s?\  m 

vc  04  o: 

3  nt  o  o  e  X  0C  <N  c 

r- 

fM  a  r-  N<tnO 

O  04  r- 

o  in  o  o  in  iNiO  r-  rr 

r* 

ic  « 

\0  ^ 

oaeoc'-acxvc 

r- 

Cu 

i-  N  N  r-  (M  X  iTi 

\o  rj 

o 

© ©  o  t,n  cni©  in  «s 

r* 

■c 

r*  oj 

as 

^  vO 

CC 

o  o  O  O  O  O  O  lC  o 

r- 

U 

IT,  (N 

Cf. 

6J 

=t  oi  n 

0. 

inOtCOO  p-  <N  r-  fn 

** 

(N  a  <N  3  0"  P- 

o 

Pc 

.J 

O 

o 

c  oj  o* 

H 

c  <c  m  c  c  o  ^  o 

* 

r- 

t- 

IN  iC  <N  *“  vC  a 

m 

u 

w. 

(C 

a. 

ONr 

t- 

lC  O  IT.  O  O  N  O  X  C 

«— 

x 

*“  IT.  «-  a  «- 

H 

►- 

•  ♦  • 

<A  ~ 

3  CO'  l'  C  Iflp  NC 

e.  *  » 

(C 

in  vc  »-  «Nr-  xr 

© 

nr  c.  sk 

o 

c 

p-rr. 

X 

< 

r“ 

3  O  n  ■£  X  N 

m 

O'  ?  JUPPX 

r> 

O' 

PPOOOlfl 

o 

r* 

x  in  ooNifin 

in 

r* 

O  X  m  r,  cx 

vr 

i 

i 

IT.  *-  IT.  m  r- 

04 

vi 

i 

i 

0:0  lINn* 

04 

z* 

i 

t 

a  a  a  ©  a  x 

cv 

O 

i 

i 

3  N  r-  ITX'3 

T~ 

m 

i 

i 

ff  lePiCX® 

X 

IT. 

i 

i 

aniocnc 

c 

*•» 

i 

t 

X  (N  X.  r«  li 

Vi 

X 

t 

i 

1C  Pvl  p»  X  C  3 

o 

ro 

i 

i 

a  ffl  a  <o  x  a 

X 

r 

i 

i 

a-o  vO  o  o  r~ 

04 

(N 

i 

i 

®  a  B3M(I 

sC 

e 

i 

i 

IT3  3  0SN 

a* 

i 

i 

OrXNlT.  n 

c* 

co 

i 

i 

CCJ  «N  X  P»  O  VO 

5f 

i 

i 

O'  CC  O'  O'  O'  O' 

C** 

ZT 

i 

i 

anxhp x 

o* 

i 

i 

a  if  a  x  c  if. 

X 

o 

i 

i 

a:ir  it  r»  r-<  <*■> 

CM 

i 

i 

p  i\  p  ill  p  ifi 

vi 

1 

i 

cc 

rvi 

i 

w 

C1  -■»  X  r-  i-  O' 

en 

i 

(5 

cs  p*  ^  ®  e '  « 

X 

i 

< 

t 

H 

c 

i 

PC 

a  n-  o  cr  in  c 

vC 

04 

hi 

x  a  o-  vc  x  x 

e 

i» 

i 

U 

•i 

w 

i 

X 

X 

*n 

i 

a 

n  ifl  o  r  •“  O 

CC 

O 

i 

a- 

o  r-  x  in.  n-  in 

o 

»r# 

«< 

i 

■e 

IX 

i 

X 

Cti 

X 

i 

c 

a  a  opns 

(c 

04 

> 

i 

M 

oi  r»  c  e*  a  cr- 

> 

c 

< 

t 

1 

1- 
l  i 

v— 

l« 

o 

o 

i 

w 

X 

o  r*  x  x  x  m 

o 

in 

U 

z* 

i 

1- 

in  m  vc  c  co 

w 

si 

H 

t 

X 

f- 

X 

c 

X 

15 

n*  »-  CN  r-  N  X 

15 

H 

<N 

t> 

n  m.  c-  r* 

H 

in 

W 

X 

W 

x“ 

tC 

*- 

r~ 

a. 

X  r*  riNOif 

«c 

X 

e 

U 

m  <r  r>  w  io  m 

X 

oc 

m\ 

04 

ICC 

*-  fN 

-l« 

o* 

C! 

H 

o 

« 

< 

in 

B-34 


TABLE  IV  37  SAMPLE  DUAL  THRESHOLD  RESULTS  (Cont’d) 


DUAI.  THRESHOLD 


APPENDIX  B-2 


X  X  IN 


VC  X  IN 


XX*- 


X  <N  IN 


W 

O'  X  IN  C3 

*■*  «c 

as 

if  VC  *»  K 

*“  U 

as 

K 

»(Nn  x 

*— 

c 

OINN  H 

*"  H 

U 

X 

ONr  H 

*“  W 

c  « 

•  *  •  *4 

as  *  •  W 

S  K  *  O 

her  e 


O'  IN  X  ©  IT.  C  P"X  O 
IN  m  VC  **k  »  IT  f>*  eo 


00  f“  O  O  O  if  IN  VO 
r»  *“  O'  it  it 


O  ©  ©  ©  O  O  OX  c 


O'  X.  XO  X  C  r-  ®  o 
in  <"  x  er  a  m  ec  O' 


<CNOCi  O  .*■  IN  C  O 
r-  «-  a  CO  X 


o*-oocccvcc 
in  in 


ifliNinooccnnn 
<n  r  o  N  n  in  to  n 


o  r*  o  o  x.  x  xx  m 

IN  (Nvflri 


O  O  O  O  O  O  O  rVX 


O'  m  m  o  o  x  r»  w  © 
<Nm*-aiNiNX0'0' 


esccoKajNO 

IN  a  f“  (N  CD  13 


o  ir  o  ©  iniNO  n  vc 

VC  ON 


Olf  COOi"  ec  CD  VO 
<N  NfflBKl 


Of  ooinnioime 
r-  in 

o  o  o  ©  c  o  c  X© 

IT.  (N 


IT  O  f  O  O  i»  it'- 0-0 
IN  a  in  :in  f"i  r*  *“  O'  vc 


O'  O'  0  .  0  ©  O'  ©  O'  X 
IN  VO  f  i-)  r-  \C  a 


vc  r»  x.  ©  ©  in  ©  oe  © 
x.  r-  m  m 


a  o-.  m  in  c  in  n  n  o 

(NvCr-  IN  rr>  v a  n>, 


f— 

a 

O' 

*— 

m 

1 

X 

1 

1 

1 

it 

1 

| 

o 

1 

1 

m 

1 

1 

in 

( 

r- 

1 

1 

X 

1 

1 

m 

1 

1 

r« 

1 

1 

IN 

1 

1 

X 

1 

1 

1 

O 

1 

I 

1 

X. 

1 

1 

in 

1 

1 

i* 

1 

1 

o 

1 

1 

IN 

I 

1 

1 

V) 

IN 

1 

M 

n 

1 

a 

<< 

1 

o 

1 

X 

m 

1 

« 

II 

1 

u 

to 

1 

C£ 

w 

m 

1 

« 

Id 

1 

a, 

< 

1 

1 

X 

w 

X 

1 

o 

► 

O’ 

1 

M 

< 

1 

h 

1 

o 

Q 

X 

1 

X 

W 

-a 

1 

l-l 

H 

» 

A! 

X 

Q 

O 

a 

K 

IB  . 

m 

H 

_  »v. 

s 

•< 

K 

as 

r« 

U 

< 

IN 

X 

u 

M 

v> 

X 

a  O  ®  X  IN  X 
O'  Ov  O'  X,  X  p~ 


O'  n  M  O  If.  OD 
x  x  x  in  en  m 


VC  VC  X -1*1  O'  x 
If  rif  «—  re 


vD  C  If  M  if  n 
O'  O'  CT1  ®  O'  O' 


a  in  n  if  o  x 
O'  vo  o'  vc  ro  ® 


sntoaoc 
x  in  f-  m  io  n 


VO  INI'  ev  o  vf 
O'  ®  O'  vc  n*  X 


o  o  n  on  o 
O'  a  co  a  m  n 


mam  otNo 
if  *•  if  in  n  m 


ec  in  ®  vc  x 
O'  X  O'  O'  O'  O' 


a  n  r,  r-  x  c* 
c  in  x  or  «  c 


e*  in  x  f~  n  o- 
r  in  n  *  x  r- 


X  a  TO  r-  0“  f 

x  n*  o>  cc  ®  o' 


a  rs  r-  O'  O  IT 
X  9  O'  X  X  X 


m  n  iN  r.  O'  I— 
x  f  <rx  m  r» 


f-  r  ffl  r-  VJ  f. 
®  in  ov  o>  x  x 


IN  f  f  f  X  C 

n  f  ®  c 


m  in  af  o  O'  « 
*-  x  x  x  x 


X  r*  l*1;<N  X  C 

x  a  r-  «  in  x 

IN  *— 


X 

*  IN 

«e  r- 
o 

M 


TABLE  IV  37  SAMPLE  DUAL  THRESHOLD  RESULTS  (Cont’d) 


\  I  ©  IN  Jt 


i 


O  /N  <~4 


O  IN  IN 


0004 


O  ®  M 


o  ac  .*m 


00  IN  Jt 


OO  Ji  OlflfOr  O 
f^<N®OXO'00'0 


4*5  «-  x  o  o  inx  in  fl 

<*5*-r,'®jtx®®o' 


MlNiOCOlTN  »  ("I 
*-  <N  <N  •-  ISI  O  C*  X 


WMOlOOfOOC) 

h«ffl0009'000 


oioaoinxiToo 
m  Jt  t^  o  x  x  x  o  o 


®  4  a  O  i/l  1(1  IN  «*>  4*5 
IN  P'1  IN  IN  rsi  st  o  as 


n  n  4  O  IP  fNI  O  o  o 
OMCCjiOOOO 


s 

e 


X  04  PI 

O'  X  X  O  X  O  O'  ®  o 

r— 

IN,  nansiflco  ji 

X  04  04 

N40O04IN00 

V 

r» 

*-  *-  *-  4t  «  f*> 

’${S  ®  it 

O  X  Jt  O  X  O  o  o  o 

\ 

O'  ®  ®  ®  O'  ®  o  o  o 

**  ^  r- 

\ 

oo  sovn 

NBOOOSOXa 

\ 

Xt-sfXXXXOO* 

®  ac  04 

\ 

O'  X  OCM  04  O'  O 

\ 

\ 

/NJtlNCNfN/NJtW® 

X  04  St 

'\’/5  0'0'000*“®0 

«■» 

U4 

\ 

X  04  X 

19 

O'  <n\X  O  O  W  O'  05  0 

« 

H 

IN  4*5  ;f  PM  04  IT  05  O' 

X 

\ 

X  04  04 

M 

04  x  O  O'-'O  ®  «  04  m 

L» 

*"  04  sf  **  04  »  X 

as 

M 

X  00  st 

tu 

x  s»  n  o  x  oi'  o  x  o 

0-  x  x  »  ®  ®  e.  co  ® 

a 

X 

r- 

-4 

o 

o 

O  ®  4*5 

04 

®  O  04  O  O  X  «t  I-*  >*5 

X 

n 

X  O  Jt  Jt  X  X  f  ®  X 

M 

a 

« 

M 

OS 

O  ®  IN 

»- 

f»XXOXXJtJtX 

X 

w 

05  ®  r»  04  <t  »•  t*  f*  St 

H 

a- 

-1 

*  «  • 

Jt  O'  O'  X  O  X  O'  04  O 

< 

BS  »  * 

u 

04  X  *“  04  05  X  4*) 

» 

X  OS  * 

t s 

a 

H  C  C 

SB 

« 

APPENDIX  B-2 


•« 

w 

w 

o 

< 

as 

w 

► 

< 

o 

w 

(-> 

s 

19 

H 

w, 

a 


U 
v ) 


Ol 

O  O'  o  o>  ®  ® 

05 

r> 

O  O'  O  ®  ®  O' 

O 

r»  r" 

X 

o  o  6  O'  x  o- 

05 

X 

O  O'  O  X  O'  ® 

T—  4— 

® 

*- 

Jt  X  X  X  O'  O' 

r- 

n 

4 

1 

O'  0*  X  4*5  X  O' 

X 

( 

4 

oo-o®oo 

O' 

® 

1 

O  O'  O  O'  o  o 

o 

! 

1 

t—  T— 

O' 

O  •—  O  ®  O'  » 

t" 

X 

4 

O  O'  O  O'  O'  O' 

o 

4 

4 

4 

r' 

X 

®  O  O  X  os  IN 

r- 

& 

O'  ®  0»  O'  Oi  O' 

O 

X 

1 

4 

O  «  ®  4*5  O'  Ol 

® 

X 

JO  O'  05'0.  O-  O' 

Ol 

1 

4 

r- 

X 

X  O  X  04  X  1*5 

IN 

Jt 

1 

®  O'  O'  ®  O'  O' 

O' 

jt  cn  n  x  Jt  ® 

r- 

Ol 

O'  X  O'  X  X  ® 

X 

Ol 

O  ®  X  O  O'  O' 

X 

® 

1 

O  O  Ol  o  Ol  Oi 

o 

r"  r- 

Jt 

1 

O  «—  X  O'  O'  I"" 

o 

O' 

O  O'  O'  O'  O  O' 

05 

o- 

tOw*l(')(*l«- 

o 

X 

M 

On  X  O'  C  O’  O' 

« 

X 

tu 

O  O  X  ®  X  ® 

,*• 

X 

1 

u 

O  O  O'  o  O  O 

O' 

*< 

T™ 

H 

O 

* 

®  IN  ®  P>  X  X 

& 

X 

w 

Ol  ®  O'  O'  o  o 

O' 

’O 

I, 

X 

w 

O' 

M 

o  o'  m -r»  ®  r» 

M 

4*5 

m 

04 

o  X  o  X  X  X 

» 

X 

■< 

1 

re 

« 

X 

O 

oxo  ®.®  X 

w 

X 

® 

04 

0X0  O'  O'  O' 

> 

O' 

H 

*—  4— 

« 

O 

(*> 

M 

X  X  X  O =r*  O 

Q 

r- 

O' 

H 

X  05  x  ®  O  o 

U 

r- 

« 

1— 

H 

43 

X 

® 

st  r-  05  st  x  o 

o 

o 

X 

H 

05  *»  9  r*  ®  © 

H 

r» 

(At  . 

•« 

—  “ri 

—  W 

» 

X 

®  r>  n  in  x  o 

r* 

U 

x  oi  r-  ®  x  x 

X 

IN 

« 

Mr- 

■< 

H 

u 

< 

w 

B-36 


| 

f 


TABLE  IV  37  SAMPLE  DUAL  THRESHOLD  RESULTS  (Cont’d) 


DUAL  THRESHOLD 


APPENDIX  B-2 


9X9 

r-r-ooinmocio 

n 

o  o  o  m  cm  ® 

in 

r- 

®  in  O  O  ©  X  O  O'  © 

® 

o  o  ©  ®  ©  © 

© 

♦—  f—  T— 

r-  r- 

9  CO  «*» 

9  ©9  Olfl  POOO 

o  n  n  ®  r»  n 

m 

r- 

iTl  (N  ®  O  ©  r*  O  O'  O 

© 

©  ©  ©  ©  c-  © 

X 

r»  r- 

9  00  M 

m  O  r»  O  U'  O  *-  (N  O 

m 

©  cm  m  ©  r-  m 

X 

f- 

fM<-90fM®r'r'r> 

9 

i 

i 

©  ®  ©  m  tn  t" 

© 

1 

i 

<-4 

9  9  9 

t-  fM  o  o  o  ®  o  o  o 

fM 

i 

o  ©  o  o  O  o 

0  X 

r* 

xxoooxooo 

© 

i 

o  o  o  ©  o  o 

0  ll 

r-  f»  r  4—  4—  f— 

i 

4--  r*  4—  44“  44-  4“ 

-  G 

9  9  Cl 

XCI9009000 

O 

i 

i 

0  ®  ®  0  ©  0 

©  O 

r- 

m  ©  ®  ©  ®  r*  o  o  o 

X 

i 

0  ©  ©  0  ©  0 

©  ^ 

4-  r“  r*  r- 

i 

4—  4—  r- 

CO 

9  9  CM 

0'9r^omrMf-iri® 

i 

i 

0  n-  r-  r»  © 

in  | 

*— 

CM  O'  9  0  9  ©  f-  O'  O' 

© 

i 

©  ®  ©  ©  ©  © 

©  5 

T- 

i 

r— 

w 

u 

(M  00  9 

inr'90o»~o®o 

o 

i 

i 

X  00  O  ©  ©  O 

X  ® 

c- 

r-  9  o»  o  or  o'  o  ©  o 

® 

i 

©  ©  O  ©  ©  O 

©  Q 

r-  r* 

i 

r-  «4_ 

Id 

i 

0 

CM  X  C| 

4-  «  ®  O  O  r-  o  n  o 

n 

i 

x  r-  0  r-  r-  r~ 

©  * 

r* 

9  ^*  r>  o  in  r-  o  o'  o 

© 

i 

©  ©  0  ®  ©  © 

©  is 

*”  r—  4— 

i 

4— 

S 

CM  X  fM 

©  X  CM  O  U1  m  r-  U1  © 

9 

i 

i 

c>  9  9  O  X  CM 

c-  ^ 

4— 

4-  9  9  9  c*  w  w 

9 

i 

i 

©  X  ©  X  © 

"  ^ 

fM  9  9 

inr*»9oor*ooo 

r-> 

i 

i 

OXOOOO 

©  § 

O'  O'  O'  O  O'  O'  O  0-0 

O' 

i 

O  ©  O  O  O  O 

O'. 

*—  *—  r*  »- 

i 

r—  — -  r—  r-  r- 

j 

►J 

(N  9  n 

OffiJOOUlOtDO 

© 

i 

x  r-  0  0  0  ® 

X  ^ 

1~ 

r»x®o®®oo'o 

® 

i 

©  ©  0  0  0  © 

O'  3 

4*™  4—  4— 

i 

4“  r*  r- 

to 

fM  9  fM 

o*  cm  x  o.tr  fs  mm® 

i 

i 

©  cm  p*  p'  ©  © 

m  p; 

n  (0  ©  ©  m  ©  ®  O'  O' 

c* 

i 

©©©©©© 

i 

i 

V) 

> 

0*9 

©  ci  ©  o  o>»-  o  o  o 

in 

i 

u 

0  r-  0  x  0  0 

O'  ci 

4-* 

r*  ©  x  ©  x  ©  ©  o  o 

® 

i 

o 

0  ©  0  ©  0  0 

O'  rj 

^  ^  C  *“ 

i 

< 

r~  1—  v- 

3 

W 

i 

H 

X 

o  ®  m 

U 

n  ©  «  ©  in  ®  in  o  o 

© 

i 

X 

O  O  X  ©  X 

r-  g 

r* 

<* 

M9noin«cDOo 

© 

i 

M 

O  ©  O  ©  ©  © 

©  H 

H 

r-  r- -t* 

44 

i 

U 

r- 

44 

X 

w 

i 

03 

M 

O  X  CM 

M 

©  ©  ©  o  o  m  fM  m-vo 

M 

X 

i 

W 

X  O  O  ©  cn  CM 

U 

r- 

r* 

u 

<N  n  fM  fM  ff>  CM  9  O'  X 

U 

9 

i 

04 

©  X  ©  ©  ©  © 

U 

© 

m 

<m 

i 

< 

M 

X 

9 

i 

X 

« 

0  9  9 

04 

o  f*  ©  o  in  ©  o  o  o 

w 

® 

i 

o 

0  in  ©  000 

© 

O  O'  ©  O  O'  ©  o  ©  o 

► 

© 

i 

M 

©  ©  0  0  0  0 

> 

© 

IB 

4—  ^  *■*  4f>  r«  «-• 

« 

i 

h 

4T»  T”  f*  4"  r" 

< 

o 

i 

u 

O  9  C> 

H 

r*  ©  ®  o  ©  in  fM  cm  © 

Q 

9 

i 

w 

©  ©  t~4  CM  CM  O 

a 

© 

H 

®  o  r*  ®  c*  «  9  x  r* 

M 

X 

i 

H 

in  n  m  ®  ©  © 

M 

r» 

U 

r- 

H 

i 

M 

r— 

H 

M 

X 

a 

X 

O  9  fM 

H 

©  m  n  o  in  r  x  <m  © 

u 

fM 

n  i-  n  0  ©  © 

© 

r- 

« 

Q 

©  O'  ©  9  m  f-  <m  r*  © 

M 

M  . 

P* 

H 

_  CL.  . 

9  fM  9  (*•  X  © 

M 

—>  W 

r~ 

■ 

< 

X 

«  •  • 

9  O'  O'  m  o  m  c-  n  o 

r— 

Oi 

X  C-  cl  CM  ©  O 

© 

BS  »  • 

W 

CM  ©  4~  cm  n  ©  n 

X 

U 

m  ©  r-  ®  m  x 

X 

CM 

X  X  X- 

ca 

< 

04 

X 

fM  *- 

■< 

© 

H  *  X 

X 

u 

H 

u 

v> 

4« e 

w 

APPENDIX  B-2 


y£>  CD  si 

*■* 

*"  m  ft  oor-rmo 
a  r-  a  oo  a  id,  ft  ®  ft 

rr\ 

in 

O  01  x  r-  ft  ft 

C  O'  9  CO  CD  (Tv 

r- 

88 

O  CO  n 
r- 

iftamoocoNNo 

<N  **aNNlflNN 

in 

vc  on  a  <c  on  r» 

(T  ffl  Jv  ID  VC  f~ 

-n 

vO  OC  (N 

r- 

0»-OOlDOCOOvOO 
<N  Nan 

vO 

r-  1 

| 

VD  CD  ON  n  vC  oo 
xaecNDi/v 

49 

^  ^  3 

*“  o  ft  o  o  n  •-  oo  o 
aNaocainn  a* 

1 

r*  i 

in  | 

i 

o  a  oo  «v«  -x 

C  T  (TIT  a  T 

ft 

97 

vc 

*■» 

ID  ID  ID  O  O  00  ft  ID  O 

ON  «“  a  CM  ON  ID  00  00 

i 

C  | 

1 

vC  ON  00  i*  pj  n 

ft  ®  ft  x  ft  ft 

o 

or 

C  2t  fM 
r- 

O  Os!  o  O  ID  00  OC  O  ft. 

ON  N  OC  a 

25 

«-  ©  ft  ID  ■—>  ID 
ft  ID  ft  X  X  X 

x 

D 

it  CO  it 

O  »-  VO  o  ID  ft  ft  00  vO 
Nf(N0Oni-ID®N 

(N  ! 

it  1 

I 

co  vc  o  ft  n  n 
ft  X  O  00  *  ft 

c» 

50 

CJ  GO 

ONOCO  r-  00  ID  VO 
<N  <N  f.<N  x  a 

cr 

M  | 

1 

n  a  mo  o  n 
oo  ft  ft  vc  .ft  x 

ft 

ft 

it  oo  rg 

COCOiDONOONvO 

VC  ON 

t***  1 

r-  | 

aoNOftan- 

ac  a  x  m  id  ft 

C 

vr 

it  Sf  « 

movcooonovo 

DaNoouiNinoio 

ft 

vC  | 

in  i 

1 

oo  vo  c.  o  ft  id 

ft  CD  O  O  ft  ft 

*—  rr 

ID 

ft 

^  ^  ^ 

r* 

CN000>-It-0 

'n^aNr-NffN 

r*  } 

t  j 

| 

i"  a  a  i:  n  n 

ft  M  ft  ft  ft  ft 

ft 

5?  it  (N 

ClTCOONOnC 
*-  or  n 

1 

X  1 

rsj  j 

> 

vc  m  n  ii  n  f. 

oo  a  ft  .t  ao  oo 

c 

MCD9 


W 

nosm  O 


(N  QD  ON  W 

*“  u 

06 
OB 

N99  ft, 

♦» 

at 

o  a  h  S 

"  I- 

U 

u 

©  a  ON  I* 

r-  W 

•  •  •  «J 

ft  •  i  W 

X  ft:  *  o 

(hex;  se 


oo*_oooa©ioi"'> 

(N  ON  *“  r-  ft  vc 


o  x  o  c  m  n  o  oo  vc 
oo  a 


O  a  ©  O  ©  ©  O  h  vC 

.a  n 


ft  r*  m  o  ID  *-  a  ft)  .4) 

(N  V£,  (N  m  ft> -r-  os  se 


o  onid  o  m  a  a  on o 

in  00=»-  ON  OO  ("I  f  ©  ID 


•*t  (N  o  c  in  h  a  a  vo 
M  N  *■  v  r  f*  |s  a 


a  o  ft  u-  o  in  n  n  o 
on  ft  *-  <n  on  vc  n 


1 

10 

r** 

1 

w 

ID  ft  O  ON  a  ID 

*  A 

n 

1 

o 

ft  a  o  c  ft  ft 

£ 

<€ 

1 

ft 

SO 

ON  X  X  r-  ft  O 

X 

ON 

1 

u 

X  tN  ft  X  X  ft 

r* 

•  • 

u 

•  t 

«/> 

fic 

V) 

w 

© 

1 

u 

a  n  r-  oco  n 

w 

cc 

c? 

ON 

1 

fr 

a  *-  ft  x  x  « 

\C 

< 

< 

05 

at 

w 

,  a 

O 

in  ft  c  e  a  rn 

w 

in 

> 

ID 

1 

H 

ft  n  o  ft  ft  ft 

-> 

X 

*c 

1 

*■< 

r- 

< 

1 

u 

o 

ft 

1 

w 

X  <N  ft  •»  *“  C 

o 

vc 

w 

ID 

1 

f"< 

in  ON  ID  X  ft  O 

t- 

1 

M 

•• 

h 

X 

C 

X 

u 

ft 

t  cr  (N  n  o 

o 

<r 

H 

a 

1- 

r>  ^  r>  oo  c 

H 

in 

w 

4<  M 

w  ^ 

m 

* 

X 

ON 

a 

x  ft  on  Dj  vo  © 

'C 

SB 

ft 

u 

in  a.  t-.  x  id  <c 

as 

(N 

«e 

ON 

K 

ON  0* 

U 

H 

u 

c/J 

« 

CO 

B-38 


TABLE  V  10  SAMPLE  DUAL  THRESHOLD  RESULTS 


APPENDIX  B-2 


v m c  ooc  c  no 

co 

CM 

a»-«®ainr*®a 

m 

vo  cxco 

in  a  uo  o  o  ®  p»  a  © 

CO 

ly 

CM  f-  a  CM  CM  l/o  VO  p. 

co 

VO  ©  CM 

O  *“  O  O  UO  ®  ®  CO  Cl 

in 

09 

CM  CM  5t  CO 

*" 

1 

vo  >o  a 

r-  p*  P-  ©  0  >  r-  o  O 

in 

1 

1 

*“ 

»f-S«9inr*BV9v 

m 

1 

1 

vOvOrt 

Ifltf  IflOOttMOfl 

« 

1 

1 

r— 

CM  c-  a  CM  CM  U0  ffl  p- 

CO 

1 

1 

VO  MS  CM 

o*~oom«®voo 

1 

1 

CM  CM  vo  CO 

CM 

1 

1 

9  0  9 

OPVCOlf.ee  C-  CO 

r» 

1 

CVI 

CM*-CM®Cir-UO»p- 

a 

a  o  ro 

O  CM  O  O  O  *-  ®  CM  CO 

vo 

» 

1 

CM 

r-  CM  CM  CM  P-  a 

CM 

1 

1 

59  O  N 

O  O  O  O  MO  CM  O  CO  re 

CM 

1 

1 

CM 

9  «“ 

r* 

9  <X>  9 

omooitccoo 

1 

r- 

CM<c-CM®CO*-mo« 

a 

1 

C* 

a  vC  CO 

OlTOOOFJJOJii 

CM 

1 

1 

_r- 

*-  cj  CM  r*  CM  ®  UO 

co 

1 

5»  Vfi  CM 

O  t-  O  O  in  CM  O  9  CO 

o 

1 

1 

1 

.  ^ 

r»  cm 

CM 

1 

- 

1 

1 

in 

I 

NO» 

®  ®  O  O  0  9  O  CO  O 

re 

1 

u 

CM 

f—  CM  c  i*  9  vo 

CO 

1 

15 

l 

«c 

U 

1 

H 

CM  O  CO 

15 

o  in  o  o  m  cm  o  m  re 

CM 

1 

X 

CM 

« 

f>*  CO 

CM 

1 

u 

H 

•1 

1 

u 

X 

1 

« 

CM  O  CM 

H 

o  o  o  o  o-o  o  a-  c 

M 

9 

1 

X 

CM 

a 

m  cm 

y» 

r» 

1 

CL, 

X 

< 

1 

K 

K 

1 

X 

CM  vo  9 

CL 

cmoo  oc  o  co 

w 

CO 

1 

o 

V* 

r>  a  CM  CM  T”  CP  VO 

► 

a 

1 

H 

o 

m 

< 

« 

H 

i-9 

o 

1 

o 

O 

©  re  co 

M 

re  a  o  o  in  ®  ©  co  vo 

a 

® 

1 

w 

X 

r- 

n 

r-  r-  CM  CM  ®  UO 

X 

a 

1 

n 

« 

u 

H 

1 

w 

M 

u 

X 

o 

« 

O  v£  CM 

CM  ©  O  O  O  CM  O  a  VO 

15 

a 

X 

V- 

w 

c*  VO  *-  p.  a 

H 

co 

o 

(U 

X 

►J 

«  •  t 

►J 

a  o  a  m  o  uo  c  cm  © 

—  X  ' 

CM 

X 

«r 

X  »  t 

M 

CM  C  C  CM  co  10  CO 

X 

f- 

u 

*•* 

*c» 

X  X  IB 

U 

«c 

CM 

X 

o 

H  «  B5 

r. 

u 

H 

< 

in 

< 

O  Cl  O  r  05  C 

vO 

©  a  o  ®  p>  ® 

r-  *“ 

« 

co  Cm  co  vo  o  ® 

vO 

a  ®  a  in  m  re 

VO 

a  cm  co  a  cm 

P» 

P-  a  ®  CM  V-  ro 

a. 

o  a  ro  a  p»  uo 

»- 

o  a  a  ®  ®  a 
•" 

a 

a  cm  p*  n  are 

o 

a  ®  a  vo  r-  ® 

® 

o  a  r-  o  p*  o 

CM 

a  a  ®  a  m  r~ 

rem  o  a  ®  ® 

in 

a  ®  o  ®  r-  ® 

® 

re  a  in  ®  a  re 

o 

ocamiop 

p' 

vo  r~  a  a  a  o 

o 

p-  a  ®  cm  co  re 

uo 

vo  vc  ©  m  re  in 

a 

a  ®  c o  a  a*  a 

a 

r- 

®  a  ®  •-  a  *- 

p* 

®  p-  a  ®  «  a 

® 

r-  o  *—  a  o  in 

uo 

®  a  a  \c  co  co 

p* 

re  c*  ®  re  a  ro 

CM 

p»  a  a  «  ®  a 

re 

ar  prcn 

® 

m  cm  a  p-  p-  ® 

*• 

re 

in 

cm  c  a  co  ®  co 

M 

uo 

a  *-  ®  a  m  p- 

15 

uo 

< 

X 

f-COIfi  CM  CO 

u 

a 

p  c,  o  a  a  a 

> 

® 

r- 

•« 

P-  f-  VO  *”  O 

o 

m. 

in  i-  m  «  a  o 

X 

r* 

r— 

H 

X 

co  a  a  co  a  c 

15 

a 

co  co  p-  ®  e 

H 

re 

_  IU_ 

X 

a  r-  co  cm  re  © 

p* 

m  e  p*  «  uo  re 

•V 

CM 

CM  *- 

P* 

u 

to 

B-39 


FT 


I 

j 

I 


i 

i 

r 

t 


i 

i 

r 

\ 

t 

f 

i 

[ 

i 

\ 

| 

l 


t 

1 

i 

r 


APPENDIX  B-2 


cm  cc  a 

in  in 

3  0  0 

*-  o  in 

O 

p* 

X  X  OX® 

L'' 

*“  T— 

IV  St 

®  O  X 

mo  x 

O 

p* 

m  m  o  p»  p* 

'r 

T3 

r- 

r- 

*“ 

f— 

cm  oc  x 

r»  in 

«  c  m 

*"  O  SO 

SC 

in 

x  m  m.  x 

*• 

in 

V—  r— 

it  *~ 

p*  ©  3 

p-  o  ic 

m 

m 

m  m  m  3  3 

X/ 

r™ 

r- 

N  X  PI 

*0  *- 

in  o  in 

in  *-  p> 

sO 

X 

p*  m  3  r  3 

*v 

sC 

s’-  »— 

f— 

3  3*“ 

3  r>  CM 

m 

CM 

1 

| 

P*  P*  -»“  f 

r 

PI  it  Of 

in  m 

3  0  0 

*  OP 

o 

P* 

1 

f 

x  x  •;  «.  x 

r 

r“  f*- 

v  it 

m  O  X 

m  o  so 

o 

P* 

t 

P  SIC  (CP 

cr 

£ 

r- 

T~ 

r* 

1 

r~ 

(N  if  CP 

•—  h 

«  o  m 

«“  c  o 

X 

X 

1 

x  P»  p*  3  a. 

s/j 

r- 

r*  r- 

5T  r» 

f“  O  3 

p*  o  r* 

m 

m 

1 

os  m  m  in  3 

/• 

r" 

r* 

1 

01,3  (N 

^  (N 

cm  o  in 

m  *T'X 

X 

X 

I 

NOSPi-C 

r*" 

»->  *- 

3  3*“ 

3  pp  x 

ii. 

CM 

1 

CC  P  05  IN  ** 

/■*• 

• 

f 

0X3 

o  r- 

moo 

•“OX 

o 

o 

1 

OPO'-IS 

■  a 

^  r* 

p*  cm 

X  ©  X 

mo® 

o 

p* 

1 

o  m  o  ®  p* 

X 

T 

«“ 

r™ 

r- 

r~  r- 

O  CC  1*1 

1*1  o 

«  o  o 

IN  Id  P 

O 

o 

i 

x  o  m  x  r~ 

»/, 

r-*“ 

.x  *- 

f“  X  3 

SO  X  SO 

m 

m 

1 

m  m  m  m  3 

O  OC  IN 

CM  r* 

SC  O  O 

in  x  m 

X 

r» 

1 

r*  x  in  0  m • 

-if 

r*  *•* 

** 

CM  IN  *“ 

X  3  X 

3 

CM 

1 

| 

PP«INr 

©3  if 

O  X 

mom 

p»  O  r- 

O 

<N 

1 

1 

0  r*  0  *-  pi 

c* 

r-  r~ 

r-  cm 

CO  O  S£ 

m  p?m 

o 

P> 

1 

0  m  0  x  sc 

r 

<5^ 

T- 

r* 

1 

r-  f— 

Of1*. 

X  *“ 

X  o  o 

x  mm 

O 

CM 

1 

0  0  x  ®  r 

• 

it 

r»  ir* 

X  ** 

r*  »  3 

so  X  X 

m 

IT 

1 

( 

0  m  m  u.  0 

f* 

0  3  CM 

X  <N 

sT  O  O 

m  pi  m 

X 

X 

1 

1 

X  X  P1  C  p 

T—  r- 

PM  CN  *- 

X  3  3 

■f 

CM 

I 

1 

xpx  iNn 

n 

un 

X  OP  it 

CM  *“ 

3  o  m 

cm  o  x 

© 

CM 

l 

w 

u 

X  X  X  P!  T 

iT 

r* 

X  CM 

X  X  3 

Sf)  O  X 

o 

X 

1 

o 

m  m  m  x  p- 

T 

X 

r— 

T” 

1 

■< 

XXX 

1 

H 

e 

01  CN 

SO  O  IT 

o  P  m 

o 

o 

* 

3  0  X  C 

r- 

r* 

< 

PM 

x  x  X 

3  in  sc 

X 

3 

w 

m  m  cp  m  m 

H 

•  • 

u 

•  * 

«  X  PI 

X 

in 

1 

OS 

O 

cu 

X  *“ 

©  O  O 

3  CM  P* 

X 

w 

X 

1 

w 

p*  x  c  0  0 

r" 

ft 

a 

*r 

U 

rr 

*-  3  X 

3 

o 

x 

cu 

®  in  x  cm  *” 

/n 

r 

© 

-< 

i 

w 

01 

i 

m 

a 

ea  f 

c. 

in  r- 

p*  o  m 

3  OX 

o 

w 

X 

i 

o 

®  •»  r-  a  - 

j 

it 

r- 

3-  r* 

mm  3 

in  o  x 

m 

-> 

m 

I 

H 

m  m  m  p*  p* 

> 

05 

c 

ss 

< 

H 

rf 

c 

i 

U 

o 

C  if  CP 

M 

O  CM 

mo  m 

3  3  X 

X 

o 

m 

i 

ft) 

X  CM  P*  «“  *- 

/*« 

{Tft 

s 

TT 

1* 

in  « 

«—  IN  x 

X  *“-06 

in 

ft) 

m 

i 

P 

m  pi  m  x  m  c. 

»*- 

w 

u 

M 

ft) 

t- 

w 

w 

X 

o 

(K 

o  ar  cm 

fl 

p»  in 

O  O  IT. 

p*  3  CM 

X 

o 

p* 

3  m  pi  x  o 

1, 

a  ► 

X 

r- 

u 

x  r* 

*-  T" 

r  »“  p* 

3 

M 

3 

H 

X  3  p*=x 

h 

\C 

H 

o 

M 

ft. 

•— 

MM 

—  »— 

-  <3 

:c 

w) 

*  •  » 

►J 

3  ©  m  m  o  m  r*  cn  c 

CM 

X 

X  P-  X  CN  X 

w 

w 

< 

CK  •  • 

hi 

nj  r* 

*-  CM  X  XX 

as 

r* 

o 

m  o>  r»  x  m. 

W 

VL 

r. 

I  K  S’! 

C 

3 

CM 

K 

n 

r* 

o 

k  »: 

7f. 

u 

M 

* 

< 

< 

* 

B-40 


TABLE  V  10  SAMPLE  DUAL  THRESHOLD  RESULTS  (Cont'd) 
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4.  M/M  Azimuth  Correlator  (Baltimore  data) 

The  additional  aircraft  data  gathered  at  Baltimore's 
Friendship  airport  sampled  40  consecutive  radar  sweeps,  rather 
than  37  alternate  sweeps  as  done  at  Milwaukee.  Results  for  a 
4/4  Azimuth  Correlator  are  given  in  Table  VI;  Table  VII  applies  to 
5/5. 

To  evaluate  the  effects  of  using  only  alternate  sweeps, 
the  processor  was  also  run  for  2/2  on  alternate  sweeps.  Table  VIII 
applies  to  the  even  samples  and  Table  IX  applies  to  the  odd  samples 
from  the  Baltimore  data. 
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5.  Dual  Threshold  Processor  (Baltimore  data) 

This  processor  was  Identical  to  that  used  on  the  Milwaukee 
data  except  that  40  consecutive  radar  sweeps ,  as  opposed  to  37 
alternate  radar  sweeps  were  processed.  Table  X  lists  the  results. 
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AZIMUTH  PATTERN  VIDEO  AMPLITUDE  DATA 

Azimuth  Pattern  Analysis 

To  help  visualize  the  type  of  returns  used  in  the 
calculations,  a  computer  was  programmed  to  plot  data.  The  amplitudes 
for  the  video  Tin  the  row  with  Che  maximum  sum  of  amplitudes)  from 
the  Data  Collect  Matrix  for  each  scan  were  used.  The  enclosed  figures 
illustrate  three  examples  of  two  scans  for  each  of  these  plots.  Each 
plot  has  target  amplitude  plotted  on  the  ordinate  and  target  azimuth 
on  the  abscissa.  The  amplitude  is  scaled  0  to  31  and  the  azimuth 
integers  correspond  to  alternate  scans.  (The  plus  signs  along  the 
edges  represent  5  amplitude  divisions).  Figures  3-16  and  B-17  are 
angel  returns  illustrating  Che  large  variety  of  returns  found  for 
angels.  These  examples  are  typical  of  the  plots  obtained  for  angel 
tracks.  The  most  striking  characteristic  of  these  plots  is  their 
width.  Whereas  the  aircraft  plots,  of  which  Figure  B-l°.  is  typical, 
appear  to  be  generally  confined  to  about  40  dwells  in  azimuth. 

The  angels  appear  to  fill  the  full  74  (twice  37)  dwells  observed, 
indicating  that  perhaps  an  extended  bird  flock  was  tracked,  or 
perhaps  several  flocks.  In  addition,  scan-to-scan,  the  angel  plots 
do  not  vary  in  a  consistent  manner,  while  the  aircraft  track  almost 
always  consists  of  one  principle  peak  with  few,  if  any,  side  peaks. 
This  peak  appears  to  be  one  of  Che  principle  reasons  for  the  success 
of  the  5/5  detector.  The  width  difference  also  explains  the  reason 
tor  degradation  in  the  discrimination  capabilities  of  other  detectors 
.'ployed,  as  the  width  of  the  sampled  region  was  decreased.  Another 
general  characteristic  was,  that  except  for  small  fluctuations,  the 
larger  fluctuations  for  aircraft  almost  tlways  vary  from  zero  to  a 
maximum  or  vice  versa.  Angels,  on  the  other  hand,  appear  to  have 
more  intermediate  types  of  fluctuation,  which  accounts  for  'he 
capabilities  of  the  dual  threshold  detector. 
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Figure  B-16  Angel  Azimuth  Patterns 
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Tigure  B-17  Angel  Azimuth  Patterns 
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Figure  B-ll  Aircraft  Azimuth  Patterns 
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ASR  ANGEL  CLUTTER  DENSITY  OBSERVATIONS 


Introduction 

Total  numbers  and  densities  of  angels  as  observed,  for 
example,  on  a  PPI  are  of  interest  because  they  are  numerical  measures 
of  what  is  meant  by  "an  angel  clutter  problem"  in  the  present  air 
traffic  control  system.  Furthermore,  they  are  indicative  of  the  angel 
clutter  which  must  be  handled  in  some  way  by  the  developing  automated 
traffic  control  systems.  This  Appendix  presents  some  quantitative 
results  on  observed  numbers  and  densities  of  angels  based  on  data 
cell acted,  for  FAA  Task  1!-  Observed  densities  depend  strongly  on  signal 
processing  and  antenna  tilt  as  well  as  on  angel  cross  sections  and 
spatial  distributions.  The  results  presented  below  are  from  ASR's 
which  were  set  up  primarily  for  a  satisfactory  tradeoff  between  aircraft 
and  ground  clutter  detections.  The  data  represent  typical  conditions 
of  moderate-to-heavy  angel  clutter  in  Milwaukee,  Little  Pock  and 
Atlantic  City  (NAFEC). 

Observed  Angel  Numbers  and  Densities 

Several  PPI  photographs  of  ASR  Mil  video  we*-e  examined; 
the  photographs,  obtained  in  preparation  for  and  during  FAA  Task  II, 
are  Identified  in  Table  1.  Tne  number  of  angels  detectable  by  eye 
were  counted  in  two-mile  range  intervals.  The  result^  are  shown  in 
f  B-19.  In  areas  of  dense°t  -lucter,  clustered  angels  cculd  not 
b«  ^s. -.‘-rarely  counted  by  eye.  In  chess  areas,  the.  count  g’ven  is  an 
underestimate  The  actual  number  of  detectable  targets  In  these  range 
intervals  could  be  as  much  as  fifty  per  cent  higher  than  the  number 
gi<>er.  The  total  counts  of  detectable  angels  are  listed  in  Table  1. 

The  identity  of  the  source  of  these  angels  is  of  some  interest  in  the 
search  for  angel  clutte»  reduction  techniques;  it  is  considered  later. 
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TABLE  1 

Angels  in  ASR  MTI  Video 


PPI  Photograph 

Date 

Time 

CSS/STC 

Total  Angels 
Counted 

Reference 

Milwaukee  Q972} 

April  15 

2122  CST 

CSS-2 

225 

a. 

April  i7 

0620 

CSS -2 

310 

a. 

April  18 

0602 

None 

284 

b. 

Little  Rock  (1972) 

March  8 

1805  CST 

STC 

230 

c. 

March  9 

0553 

STC 

167 

c. 

NAFEC  ( 1971 ) 

October  28 

1931  EDT 

STC 

225 

d. 

References : 

a.  Figure  5. 2. 1-1,  MSO-F-133,  Task  II  Interim  Report,  22  May  1972 

b.  Figure  5. 2. 1-3,  MSO-F-133,  Task  II  Interim  Report,  22  May  1972 

c.  On  file;  aimilar  to  Figures  2  and  3,  MFS-0-197,  23  March  1972 

d.  Figure  2,  MRD-3-394,  10  November  1971 
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MILWAUKEE  »  a  *  #  j 


0  2  4  6  8  10 


RANGE 

(MILES) 

Figure  B-19  Number  of  angels  counted  per  range  interval  in  PPI  photographs 

of  ASR  MTI,  video.  Probable  underestimates  shown  by  *  (see  text). 


APPENDIX  B-4 


Mean  angel  densities  were  determined  from  the  data  of 
FigureB-lS.by  dividing  the  number  of  angels  counted  in  a  range  interval 
by  the  area  of  the  interval.  Table  2  gives  the  results.  Included 
in  the  Table  are  results,  obtained  in  the  same  way,  for  several 
additional  photographs  of  Little  Rock  video.  These  photographs  were 
taken  within  minutes  of  ones  listed  in  Table  1. 

Mean  densities  in  Table  2  are,  for  the  most  part,  on  the 
order  of  1  per  square  mile  with  a  range  of  about  two  orders  of  magnitude. 
The  minimum  densities  shown  represent  a  few  angels  in  88  square  miles. 
Inside  six  miles,  these  densities  must  he  considered  underestimates 
for  two  reasons.  First,  as  indicated  above,  where  density  is  high, 
some  echoes  appear  to  be  merged.  An  effective  automatic  detection 
system  or  an  optimized  combination  of  PPI  and  camera  would  probably 
often  show  two  distinct  angels  where  the  visual  counting  yielded  one. 
Second,  angels  are  usually  not  distributed  uniformly  in  azimuth 
Large  sectors  of  a  particular  range  interval  may  contain  few  or  no 
angels  while,  in  other  sectors,  the  mean  density  is  several  times  the 
tabulated  value. 

To  obtain  more  detailed  information  on  angel  densities,  two 
photographs  of  Little  Rock  MTI  video  were  further  analyzed.  A  square 
grid  was  superimposed  on  each  photograph  and  the  number  of  angels 
detectable  in  each  cell  was  counted.  The  grid  and  results  are  shown 
in  Figure  20.  (Note  that  cell  size  in  miles  differs  in  the  two  cases.) 
The  azimuthal  nonun.ifor.tity  which;  was  mentioned  above  is  evident. 
Although  detection  density  generally  decreases  with  range,  the  variation 
is  not  smooth.  In  Figure  20a,  the  distribution  appears  to  peak  about 
four  to  six  miles  from  the  radar.  These  distributions  are  fully 
consistent  with  the  known  customary  paths  of  the  birds  believed  to  be 
the  source  of  these  angels. 
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* 

TABLE  2 

-2 

Observed  Mean  Angel  Densities  (mile  ) 
in  Range  Intervals 


Interval  '(miles! 


Ml lwaukee 

0-2 

2-4 

4-6 

6-8 

April  15 

>16 

3,3 

1.0 

0.4 

17 

11.7 

2.7 

0.5 

0.3 

18 

9.9 

2.6 

0.9 

0.07 

Little  Rock 


March  8 

1805 

6.8 

2.7 

0.7 

0.03 

1812 

8.3 

2.4 

0.5 

0.03 

March  9 

0552 

3.1 

1.5 

0.7 

0.1 

0553 

3.9 

1,4 

0.8 

0.2 

0556 

3.3 

2.4 

0.6 

C.4 

0648 

5.0 

1.8 

0.7 

0.1 

NAFEC 

October  28  3.6  1.7  0.7  0.4 
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Figure  B-20  Angels  counted  in  PPI  photographs  of  Little 

Rock  ASR  MTI  video.  Underestimates  are  shown 
by  *. 

a.  Cell  size  1.9  x  1,9  miles.  0556CST,  9  March  1972 

b.  Cell  size  1.56  x  1.56  miles.  1805C3T,  8  March  1972 
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Cumulative  distributions  of  the  vie an  densities  calculated 
from  the  data  of  Figure  20  are  given  in  Figure  21.  Maximum  angel 
densities  observed  were  near  10  mi  2.  Note  that  this  value  depends 
on  the  cell  size  chosen.  If  smaller  cells  had  been  used,  higher 
maximum  densities  would  probably  have  been  found.  Mean  densities 
taken  over  all  non-empty  cells  were  2.4  mi"2  over  105  square  miles 
and  4.4  mi”2  over  58  square  miles. 

Angel  Identification 

Far  the  data  given,  a  strong  circumstantial  argument  can 
be  made  that  the  angels  were  due  to  birds.  At  Little  Rock,  the 
appearance,  disappearance  and  approximate  spatial  distribution  of  angels 
was  very  closely  correlated  with  easily-observed  mass  movements  of 
blackbirds  near  the  airport.  Local  air  controllers  confirmed  that 
this  correlation  was  normal.  At  Milwaukee,  migratory  activity  in  the 
area  was  common  and  widespread  during  the  times  of  observation.  At 
NAFEC ,  the  area  of  angel  activity  was  a  bird  sanctuary  known  to  be 
active,  and  early-evening  bird  activity  was  a  regular  occurrance. 

Ir.  all  cases,  angels  which  could  be  followed  on  the  PPI  behaved 
like  birds  and  not  like  ground  clutter,  air  or  ground  vehicles,  or 
recognizable  atmospheric  phenomena.  Weather  was  clear  and  stable 
enough  in  each  case  to  be  eliminated  as  a  plausible  clutter  source. 
Satisfactory  MT1  cancellation  of  stationary  targets  and  ground  clutter 
was  verified  before  and  after  angels  became  numerous  and  was  verified 
In  clear  sectors  and  at  longer  ranges  while  angels  were  present. 

These  latter  observations  also  confirmed  that  equipment  noise  was  not 
a  factor.  The  conclusion  dsawn  on  the  identity  of  the  observed  angels 
is  that,  with  high  probability,  almost  all  were  due  to  birds.  Although 
some  may  have  been  due  to  individual  birds,  the  characteristics  of 
typical  bird  activity  at  each  site  during  the  period  when  observations 
were  made  suggests  that  most  angels  were  probably  due  to  groups  of 
several  birds  or  more. 
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120 


30 


1 
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Figure  B-21  Cumulative  distribution  of  mesn  eneel  densities  in 

Little  Rock  ASR  MTI  video,  @«nd(t>)  refer  to  FigureB-20, 
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Summary 

In  PPI  photographs  of  MT1  video  at  three  ASR  sites,  typical 
angel  densities  on  the  order  of  1  mi“2  were  found  in  moderate-to-heavy 
clutter.  Maximum  mean  densities  of  approximately  10  mi“2  were  observed 
in  averages  over  areas  of  six  to  eight  square  miles.  The  total  number 
of  detectable  aagels  ranged  between  approximately  150  and  300.  With 
hig,i  probability,  the  angels  were  due  to  birds  in  groups. 
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ANALYTIC  MODEL  OF  PPI-OBSERVER  INTERACTION 


Abstract 


An  analytic  model  of  PPI-observer  interaction  is 
investigated  with  a  view  toward  evaluating  the  impact  of  Clutter 
Reducing  Devices  (CRD)  on  the  detectability  of  selected  target  typer. 
Optimum  conditions  are  found  within  this  model  for  the  reduction 
of  clutter  masking  effects  and  several  properties  of  this  system 
are  illustrated  using  particular  examples  of  the  model. 

The  results  indicate  that  as  the  clutter  component 
increases  the  strength  of  the  clutter  reducing  device,  as  measured 
by  the  percentage  of  valid  targets  rejected,  must  be  increased  to 
produce  an  Increase  in  target  detectability:  l.e.,  although  a 
greater  percentage  of  valid  targets  are  rejected  by  the  CRD,  more 
targets  become  detectable  due  to  a  reduction  in  clutter  masking 
effects.  Also,  as  the  clutter  Increases,  the  gain  in  target 
detectability  for  this  optimum  CRD  Increases.  This  result  can  be 
misleading,  however,  as  an  improvement  of  from  one  to  ten  targets 
detectable  out  of  a  thousand  possible  targets  amounts  to  a  ten  to 
one  increase  in  detectability  but  is  not  sufficient  improvement 
for  the  conditions.  Thus,  eventually  this  "real  life"  device  also 
"saturates"  as  far  as  usefulness  is  concerned.  Finally,  the 
receiver  will  be  shown,  in  the  models  considered,  to  actually 
Increase  target  detectability  slightly. 


Introduction:  Clutter  Elimination  and  Target  Detectability 

In  designing  a  radar  system,  the  objective  is  always  to 
obtain  the  maximum  detectability  of  some  given  target  type. 

Increasing  detectability,  however,  implies  increasing  the  sensiti/ity 
of  the  radar.  Increasing  clutter  which  can  reduce  detectability  by 
obscuring  or  masking  valid  targets.  Thus,  tradeoffs  must  be  made 
between  the  cost  of  n**t  detecting  a  tir^jet  and  the  amount  of  clutter 
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masking  that  can  be  tolerated  in  order  to  maximize  the  probability 
of  detecting  a  target.  It  is,  however,  possible  to  use  the 
characteristics  of  the  expected  targets  and  clutter  to  build 
devices  which  will  selectively  reject  clutter  while  retaining  valid 
targets,  increasing  detection  probability.  Ideally  such  a  device 
would  reject  only  clutter  and  pass  only  targets,  so  that  the 
sensitivity  of  the  radar  could  be  made  as  high  as  desired  and 
target  detectability  could  be  chosen  at  will.  Unfortunately,  real 
life  devices  will  always  pass  some  clutter  and  reject  some  targets 
so  that  sone  reduction  in  actual  number  of  targets  on  the  PPI 
occurs,  and  some  targets  are  still  masked  by  clutter  when  the 
device  is  incorporated  into  the  system.  If  this  Clutter  Reducing 
Device  (CRD)  is  designed  properly,  however,  large  numbers  of 
targets  previously  obscured  by  clutter  will  now  be  easily  detected 
so  that  an  overall  increase  in  detectable  targets  occurs.  In 
addition,  as  a  large  amount  of  clutter  has  been  eliminated,  the 
sensitivity  of  the  system  can  be  increased  until  the  clutter 
passed  by  the  CRD  reaches  maximum  tolerable  levels .  With  proper 
design  this  can  compensate  for  some  of  the  targets  rejected  by  the 
CRDs  (perhaps  even  increase  the  efficiency  of  the  CRD)  and 
contribute  to  the  overall  increase  in  detectability  of  the  desired 
targets.  The  problem  comes  in  finding  the  parameters  for  the  CRD 
which  will  maximize  this  detectability,  particularly  in  view  of  the 
probability  of  rejecting  a  valid  target  during  the  processing. 

If  one  is  too  worried  by  the  possibility  of  the  CRD  rejecting  a 
valid  target,  the  parameters  may  be  chosen  to  minimize  target 
rejections  so  all  targets  appear  on  the  PPI.  Unfortunately, 
this  situation  also  implies  that  a  large  amount  of  clutter  will  be 
passed  by  the  CRD  and  obscure  many  of  the  targets  on  the  PPI, 
producing  the  same  effect  as  a  rejection.  Thus,  little  or  no 
improvement  is  obtained.  The  opposite  extreme  also  yields  no 
improvement.  However,  if  the  designer  is  willing  to  experience 
the  loss  of  a  few  targets  due  to  CRD  rejection,  many  other  targets 
may  become  detectable  yielding  net  increase  of  targets  detected. 

The  choice  of  parameters  determines  which  case  occurs  and,  by  a 
proper  choice  of  the  operating  parameters,  the  maximum  improvement 
for  a  given  device  can  be  obtained.  (Figure  1  shows  how  a  CRD 
might  be  incorporated  into  a  radar  system). 

The  purpose  of  this  memo  is  to  provide  an  analytic  approach 
for  optimizing  the  parameters  of  a  CRD  in  the  above  situation  and  to 
illustrate  the  results  for  a  particular  model  of  aircraft  detectability 
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for  a  man-radar  system.  In  particular,  the  characteristics  of  a 
device  designed  to  detect  aircraft  amidst  angel  returns  will  be 
used  to  characterize  the  decision  making  scheme.  The  concept  of 
this  device  was  derived  from  the  results  of  angel  clutter  study 
conducted  for  the  FAA  and  so  represents  a  possible  system. 

General  Considerations 

Basically  the  problem  is  formulated  as  follows.  At  any 
given  time  there  is  generally  some  finite  number  of  aircraft 
targets  (N^c°)on  a  PPI  display.  Due  to  the  interaction  of  clutter 

with  these  targets  and  to  basic  limitations  in  the  skill  of  the 
operator  observing  the  PPI,  some  targets  go  undetected  by  this 
observer.  In  principle,  however,  all  targets  could  be  found 
given  sufficient  time  and  information.  The  actual  number  of  targets 
detected  (N^J  is,  therefore,  some  fraction  of  the  actual  number 

of  targets  and  can  be  related  to  this  actual  number  by 

"Aij-»ac°  f(1W’  (l> 

where  F(x)  represents  the  observer's  skill  in  detecting  targets  as 
a  function  of  a  clutter  parameter  (NPT).  Here  F(x)  is  a  function 

with  range  £o,  1 

appearing  on  the  PPI  prior  to  application  of  the  CRD.  The  nature 
of  F(x)  depends  on  the  capabilities  of  the  PPI  observer  but  can 
be  generally  expected  to  be  a  monotonically  decreasing  function  of 
x.  This  follows  from  the  assumption  that  an  increase  in  the 
clutter  will  lead  to  a  reduction  in  the  number  of  targets  detected 
by  the  observer.  Thu3  in  general,  for 


X  -  Y 

(2a) 

F(x)  _  F(y) 

(2b) 

with 

0  J.  F(x)  <  F(y)  ^  1 

(2c) 

and  Ncl  is  a  measure  of  the  clutter  density 
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If  now  a  CRD  is  incorporated  into  the  system,  the  number 
of  targets  actually  on  the  PPI  is  reduced  to 

<V>'  -  V  £r  (3) 

while,  at  the  same  time,  the  clutter  has  been  reduced  to 


(!W’ 


CL 


W 


Here  1  -fr  and  1  ~fr'  are  the  fraction  of  targets  and  clutter, 
respectively,  removed  by  the  CRD.  Both  of  these  fractions  are  a 
function  of  the  design  parameters  at  the  CRD  and  may  vary  between 
0  and  1  as  determined  by  the  operating  characteristics  of  the 
devi ce . 


Presumably  the  detecting  capabilities  of  the  PPI  observer 
are  unchanged  by  the  CRD,  so  the  number  of  targets  actually  detected 
by  the  observer  should  be  related  to  the  targets  actually  displayed 
in  the  same  manner  as  prior  to  Installation  of  the  device,  i.e., 

(NAC>’  ■  <NAC°>'  F«“cL>,>  ■  fr  »AC°  F(£r'  KCL>  <5)  . 

The  object  of  building  such  a  device  is  '-hat 


AC 


(6) 


i.e.,  more  targets  are  detected  with  the  device  than  without  it. 

In  other  words, 

(NACr 

Improvement  factor  *  — - -  ••  1  (7) 

AC 

which  in  terms  or  characteristics  of  the  observer  and  a  CRD,  becomes, 


fr  F(fr*  Ncl) 
F  <HCI> 


(8) 


If  this  inequality  is  not  satisfied  for  a  CRD,  then  it  is  actually 
decreasing  target  detectability  and  should  no  longer  be  considered. 
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Note  that  the  dependence  on  the  desired  target  numbers  has  cancelled 
out,  so  that  only  the  operating  characteristics  of  the  CRD  and  the 
observer  as  functions  of  the  initial  clutter  density  remain.  The 
optimization  occurs  by  choosing  the  values  of  fr  and  fr’  which 
maximize  thi3  ratio.  As  fr  and  fr’,  which  describe  the  operating 
point  of  the  CRD,  are  dependent  parameters,  differentiation  with 
respect  to  one  of  these  parameters  will  lead  to  a  maximizing  condition. 


i.e. , 

d 

***  1 

l-l 

l _ 

[£t’  scl1  ' 

d  fr 

L 

This 

leads  simply 

to 

fr  X(fr*  Ncl) 

ncl 

d  fr 
d  fr 

where 

1 

r 

d_ 

dx 

F(x) 

=  0 


x*a 


(9) 


(10) 

(11) 


As  F(x)  is  monotomic  decreasing, ^(x)  is  positive,  although 
the  exact  form  depends  on  the  model  assumed  for  the  observer.  In 
view  of  this  result,  mrximization  of  the  detectability  require.')  that 
the  operating  parameters  of  the  CRD  must  not  be  selected  to  place  the 
operating  point  on  a  portion  of  the  CRD’s  characteristic  operating 
curve  with  negative  slope  for  all  other  terms  in  equation  10  are 
positive.  The  exact  point  of  operation,  of  course,  depends  on  the 
exact  characteristics  of  the  device  used.  Once  both  functions  are 
known,  however,  the  exact  optimum  operation  point  can  be  found. 


The  second  criteria  for  a  maximum,  (the  third  criteria  for 
an  operational  CRD),  namely 


d2 


d  fr' 


fr 


F  1,CL) 

F  «o> 


<  0 


(12) 


leads  to  the  expression 


y  (fr'  NCL}  <  f T (Ncl  g)2 


2  1  d_ 

fr  g  d  fr 


S 

/ 


(13) 
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where 


y(x)  =  F(k)  dx7  F(x) 


If  this  condition  is  not  satisfied  no  real  optimum  operating  point 
exists  for  the  clutter  conditions  and  the  device  provides  no  benefit 
to  the  observer.  In  this  case,  the  optimum  possible  situation  is  to 
switch  the  device  from  the  circuit,  and  allow  the  observer  to  select 
targets  on  his  own. 

Exponential  Observer 

To  illustrate  the  technique  just  described,  it  will  be 
useful  to  assume  an  exponential  model  for  the  operating  characteristics 
of  the  observer.  The  exact  relation  assumed  is 

F  (NCI)  -  e  u*  (14) 

where  r  is  a  parameter  describing  the  operator’s  clutter  handling 
capability.  As  one  operator's  "heavy  clutter"  is  another  operator's 

"light  cluttec",  only  the  ratio  N^/ x  allows  some  measure  of  the 

"clutter  situation"  for  each  observer.  Thus,  when  this  ratio  is 
large,  the  clutter  level  is  sufficient  to  cause  the  operator  to  lose 
many  targets,  although  another  operator,  viewing  wit'  the  same 
clutter  level,  may  have  enough  experience  to  result  in  a  smaller  ratio 
of  N  /t  and  therefore,  a  smaller  fraction  of  lost  aircraft.  This 

ratio  then  determines  which  is  heavy,  and  which  is  light  clutter  for 
each  observer  and  is  the  characteristic  parameter  describing  the  man- 
clutter  situation.  The  lower  curve  in  Figure  2  illustrates  this 
exponential  function  where  regions  of  light,  moderate,  and  heavy 
clutter  have  been  more  or  less  arbitrarily  assigned. 


In  principle  a  CRD  can  be  constructed  to  separate  any  twe 
(or  more)  target  types.  By  way  of  example,  however,  it  will  be 


™  A'sr  *  ^  r-'T  ^T^J* "' 


t  -g^g ^yrjgtgng^pyy^^^y^g^fK^^y  "wfl'1^ 
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assumed  that  the  desired  targets  are  aircraft,  and  the  clutter  Is 
produced  by  angel  returns  generally  from  birds  and  anomonalous 
propagation  effects.  A  device  capable  of  discriminating  between 
these  two  types  of  returns  based  on  the  azimuth  run-length 
characteristics  of  the  video  has  been  devised  in  response  to  a 
study  program  conducted  by  APL  for  the  FAA.  Figure  3  is  a  plot 
of  the  operating  curve  for  this  device  as  derived  from  actual 
data  collected  from  angel  and  aircraft  video  returns.  Also 
plotted  is  a  hyperpollc  fit  to  this  data  given  by 


fr,  .  .0636 
11  1.06-fr 


-0.6 


(15) 


which  is  the  desirec  relationship  between  fr  and  fr'  for  this 
device. 


Using  these  relationships,  the  improvement  in  target 
detections  can  be  calculated.  Figures  4  and  5  show  improvement 
curves  for  both  moderate  and  heavy  clutter  situations  with  this 
approximation  of  the  CRD  operating  curve.  The  improvement  factor, 
(N.-'/N.p)  can  be  seen  to  decrease  below  one  for  some  values  of 

fr  indicating  that  this  region  of  operation  is  to  be  avoided  as 
more  aircraft  are  lost  by  the  device  than  are  gained  by  the  reduced 
clutter  level.  The  optimum  operating  point  is  that  which  maximizes 
the  improvement  factor,  and  is  indicated  on  each  plot.  Both  plots 
also  indicate  that  both  the  maximum  improvement  factor  and  the 
optimum  operating  point  change  with  changing  clutter  conditions. 
Figure  6  shows  a  plot  of  the  maximum  improvement  factor  as  a  function 
of  the  clutter  parameter  N^/t.  Interestingly  k.hi',  device  yields 

little  improvement  in  light  or  moderate  clutter  but  results  in 
considerable  improvement  in  severe  or  heavy  clutter,  although  this 
is  model  dependent.  Figure  7  shows  how  the  optimum  operating  point 
for  this  CRD  varies  with  clutter  conditions.  To  implement  this 
device,  some  means  of  adjusting  the  operating  point  must  be 
incorporated  30  as  to  allow  adjustment  to  the  optimum  parameters 
for  the  clutter  situation  as  required  by  this  curve. 
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Finally  returning  to  Figure  2,  the  upper  curve  represents 
the  new  operating  curve  for  the  radar-CRD-man  system  and  indicates 
the  maximum  improvement  which  can  be  achieved  within  this  model  for 
this  device,  This  improved  curve  is  very  much  model  dependent  as 
Figure  8  illustrates.  Here  the  same  calculations  have  been 

-x2 

performed  with  an  e  model.  In  both  esses..,  the  horizontal  dashed 
line  represents  the  ideal  case  of  no  target  losses  regardless  of 
clutter  conditions  and  the  closer  the  final  '.-urve  approaches  this 
line,  the  mote  effective  is  the  CRD. 

It  is  perhaps  of  interest  to  consider  these  two  sets  of 
curves  in  more  detail  in  order  to  gain  some  appreciation  of  the 
effect  of  this  process.  Table  1  illustrates  the  pertinent  points, 

-x2 

For  the  e  model  in  moderate  clutter  and  before  application  of 
the  CRD,  the  radar-man  system  is  finding  78%  of  targets  actually 
there  Afterwards,  the  radar-CRD-man  system  finds  92%  of  the 
original  number  of  targets.  As  the  clutter  level  increases,  the 
improvement  also  increases,  although  as  can  be  seen  for  the  simple 
exponential  racdel,  an  improvement  of  a  factor  of  10  to  1  is  not 
enough  as  only  about  51%  of  the  targets  are  detected  after  application 
of  the  CRD. 

There  exists  an  intriguing  special  case  of  a  CRD  for  which 
fc  =  fr'  describes  the  operating  curve.  This  device  is  utilized  almost 
daily  in  every  radar  system  in  existence.  It  is  the  familar  gain 
control  Interestingly,  if  the  gain  is  reduced,  both  targets  and 
clutter  are  reduced  equally,  however,  with  either  of  these  models, 
due  to  the  non-linear  characteristics  of  observer,  there  is  an 
actual  increase  in  the  detection  capabilities  of  the  system.  (This 
cf  course  assumes  no  difference  in  the  amplitudes  found  for  the  clutter 
and  the  targets.)  Figures  9  and  10  are  the  corresponding  curves  for 
the  two  models  considered  with  the  gain  control  used  as  a  CRD.  It  is 
clear  that  some  improvement  can  be  achieved  by  decreasing  the  gain 
but  only  in  very  heavy  clutter  situations  and  even  then  the  fraction 
of  targets  detected  is  still  miniscule.  Although  if  it.  is  absolutely 
necessary  co  operate  with  such  high  clutter  densities,  then  by  all 
means  turn  the  gain  down,  as  the  improvement  is  significant. 

If,  in  fact,  amplitude  is  a  useful  discriminate  between 
clutter  and  valid  targets,  the  gain  control  will  act  as  a  CRD  and 
even  greater  improvement  factors  can  be  achieved  for  such  targets. 
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TABLE  1 

FRACTION  OF  OBSERVABLE  AIRCRAFT 


Model 

Clutter  Condition 

-V 


-x 


Before 

CRD 


After 

CRD 


-x 
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Before 

CRD 


Af  te 
CRD 


Moderate  Clutter 
.5 


.61 


.78 


.78 


-92 


Heavy  Clutter 
3 


.05 


.51 


1.2x10 
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As  It  turns  out,  this  is  indeed  the  esse  for  angels  and  aircraft,  as 
the  great  preponderance  of  angels  are  ^ery  weak  targets.  Thus,  reducing 
the  gain  can  selectively  reduce  clutter  returns,  increasing  the 
detectability  of  aircraft  as  described  above.  Figures  9  and  10 
therefore,  represent  the  minimum  possible  improvements  within  these 
models  for  the  gain  control  CRD,  with  more  improvement  possible  if 
clutter  and  targets  can  be  reasonably  separated  by  amplitude. 

Other  models  of  the  operator-receiver  system  and  other  CRD's  will 
yield  similar  conclusions  and  although  specific  results  will  be  model 
dependent,  it  is  clear  that  with  a  mono tonic  decreasing  non-linear 
operator  detection  curve  similar  maximization  effects  can  be  found 
for  realistic  CRD's.  Such  devices  will,  in  fact,  increase  the  number 
of  targets  which  can  be  detected.  The  price  is  the  loss  of  a  few 
targets  rejected  by  the  device,  but  this  loss  can  be  made  smaller 
than  the  increase  in  targets  newly  detectable,  resulting  in  a  net 
increase  in  detectable  targets. 

Adaptive  Systems 

The  concern  so  far  has  been  with  the  effect  of  clutter 
masking  as  viewed  over  the  whole  PPI.  Normally  occurring  clutter, 
appears  only  in  subregions  of  the  PPI,  with  other  subregions  remaining 
relatively  clear.  This  means  that  subregions  of  the  PPI  display 
should  be  individually  processed  by  a  CRD  operation  with  the  optimum 
parameters  for  that  subregion.  This  effect  can  be  achiev’d  by  use 
of  various  adaptive  feedback  systems.  What  is  needed  is  some  measure 
of  t.he  clutter  intensity  as  indicated  by  the  ratio  Ncl/t. 

It  is  possible,  by  monitoring  both  the  input  and  output 
signals  of  the  CRD  to  calculate  the  clutter  parameter  N^.  Figure  A 

below  illustrates  the  operation  of  the  CRD  in  terms  of  number  of 
input  hits  (N°)  and  number  of  output  hits  (N') 


N 


N  °  +  N  ° 
AC  CL 


Input 


N' 


CRD 


V  +  ncl' 


fr  N 


AC 


Output 


+  fr' 


FIGURE  A 
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These  exact'  relationships  are  undoubtedly  too  complicated  for  quick 
evaluation,  thus  some  approximation  will  most  likely  be  made  in 
practice.  In  principle,  however,  the  feedback  can  be  accomplished 
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exactly  within  these  models  using  this  procedure.  The  feedback  delay 
between  the  actual  calculation  of  N°  and  N*  and  the  actual  modification 
of  the  operating  parameters  of  the  CRD  will  depend  on  the  type  of 
clutter  considered,  although  a  delay  of  one  or  two  sweeps  would  appear 
to  be  appropriate. 


Target  Identification 


For  many  purposes  it  is  not  desirable  to  merely  eliminate 
clutter  from  the  PFI-observer  system.  The  elimination  of  the  clutter 
from  the  system,  as  shown,  can  improve  target  detectability,  but  can 
then  result  in  the  observer  vectoring  aircraft  into  the  clutter 
regions  on  the  PPI  which  have  been  cleared  by  the  CRD.  If  the  clutter 
consists  of  birds  as  in  the  example  considered,  the  results  can  be 
particularly  disastrous.  Thus  it  is  particularly  desirable  to  display 
the  clutter  on  the  PPI.  This,  in  turn,  would  reduce  target  detectability, 
at  least  with  present  display  systems. 


The  obvious  conclusion  is  to  hav  .  some  unique  method  of 
identifying  targets  and  clutter  on  the  PPI  display  (such  as  color, 
or  symbology).  Target  detectability  is,  after  all,  really  “target 
identification".  The  analysis  so  far  has  considered  either  the 
presence  or  absence  of  clutter  and  targets  and  their  interaction  as 
the  primary  effect.  Actually,  the  same  effects  and  conclusions  can 
be  drawn  if  the  problem  is  viewed  as  a  target  identification 
problem,  with  all  returns  continuously  on  display,  but  separated  now 
by  various  display  techniques  into  target  types.  Now  the  display 
of  the  clutter  targets  can  prevent  the  vectoring  of  aircraft  into 
clutter  targets,  measurably  improving  system  performance  and  safety. 


Generalization 

Although  the  presentation  has  been  covered  in  terms  of  the 
interactions  between  a  radar  system  and  its  ovserver,  the  conclusions 
actually  are  more  general,  applying  to  any  information  handling  system 
where  some  desired  information  must  be  sorted  from  undeoired  clutter 
information.  In  general,  this  unwanted  clutter  information  will 
mask  the  actual  Information,  reducing  information  transmittal  rate  of 
the  system  and  increase  the  need  for  redundancy  to  assure  the  complete¬ 
ness  of  information  actually  transferred.  Thus,  this  more  general 
system  may  be  described  by  a  monotonically  decreasing  function 
describing  the  fraction  of  the  original  information  received  from  the 
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system  by  the  observer  in  the  presence  of  clutter.  This  function  will 
satisfy  a  relation  equivalent  to  equation  1  where  the  symbols  now 
represent  information  content  instead  of  target  counts.  As  a  result, 
the  derivations  of  the  previous  sections  can  be  developed  for  this 
system  and  similar  conclusions  will  be  found. 

For  example,  if  there  exists  some  sort  of  clutter  reducing 
device,  which  of  course  will  also  eliminate  some  of  the  desired 
information,  and  if  it  is  possible  to  satisfy  equation  8  with  this 
device,  then  the  information  transmitted  can  be  increased  by  reducing 
the  effects  of  clutter  masking.  This  indicates  how  the  generalization 
to  arbitrary  systems  can  be  developed,  although  the  existence  of  the 
device  with  the  required  parameters  is  not  always  assured. 
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ANALYSES  OF  AREA  MTI  PERFORMANCE  FOR  ANGEL  CLUTTER  REDUCTION 


Summary 


An  analytical  model  of  an  Area  MTI  is  derived.  The  model 
reflects  the  interactions  of  random  bird  velocity,  single-scan  thres¬ 
holding,  different  area  MTI  designs,  dense  angel  clutter  effects, 
and  operator  masking  effects.  A  practical  baseline  design  is 
evaluated  with  the  model,  resulting  in  the  conclusion  that  the 
Area  MTI  improves  on  single  scan  processing  only  if  1)  angel 
densities  are  light,  and  2)  aircraft  detection  by  an  operator 
would  be  severely  masked  by  angel  clutter  breakthrough. 


2.  Introduction 

Scan-to-scan  or  Area  MTI  is  a  technique  which  discriminates 
against  slow  moving  and  stationary  targets  by  means  of  the  change 
in  their  position  from  scan-to-scan.  The  measured  position  of  all 
declared  targets  for  one  scan  is  stored  in  a  digital  memory.  On  the 
next  scan,  the  measured  positions  of  new  targets  are  compared  with 
the  stored  position;  if  there  are  any  scan-to-scan  correlations, 
corresponding  to  low  velocities,  the  new  data  is  not  displayed. 

For  the  ASR  scan  rate,  a  60  knot  target  moves  410  feet  (one  pulsewidth) 
from  scan  to  scan. 

While  Doppler  or  pulse-to-pulse  MTI  discriminates  only  on 
the  basis  of  radial  velocity  (range  rate),  Area  MTI  discrimination 
involves  total  velocity  in  any  direction.  Since  bird  velocities  are, 
on  the  whole,  significantly  smaller  than  aircraft  velocities.  Area  MTI 
seems  to  be  of  potential  value  in  reducing  angel  clutter.  Area  MTI 
can  be  interfaced  with  single-scan  processing  (such  as  5/5  correlation) 
and  an  operator  as  shown  in  Figure  1. 

There  are  two  main  types  of  Area  MTI,  with  several  variations 
on  each  type:  Sector  Map  and  Target  Store. 

The  Sector  Map  Area  MTI  divides  the  surveillance  area  into  a 
grid  on  non-overlapping  sectors;  each  range-bearing  sector  (or  cell)  is 
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represented  by  one  bit  in  a  sector  map.  If  a  target  is  detected 
in  a  certain  sector,  the  storage  bit  for  that  sector  is  set  to  1; 
otherwise,  the  bit  is  reset  to  0.  When  the  data  for  the  new  scan 
comes  in,  the  cell  is  displayed  only  if  the  old  scan-new  scan  bit 
pattern  is  0,  1. 

The  sector  map  can  be  implemented  fairly  easily.  If  cell 
quantization  of  1/16  nmi  by  0.75°  is  used  in  a  map  out  to  15  nmi, 
chon  digital  storage  of  115,200  bits  is  required.  In  addition, 
fairly  simple  logic  can  be  constructed  to  report  the  presence  of  a 
hit  within  a  cell. 

The  Target  Store  Area  MTI  stores  the  range  and  bearing  of 
the  centroid  of  each  target  detected  during  a  scan.  When  a  target  is 
detected  on  the  new  scan,  data  from  the  old  scan  are  searched  to  find 
a  target  within  a  range-bearing  correlation  window  corresponding 
to  low  velocity.  If  a  correlating  target  is  found,  the  new  target 
is  not  displayed. 

The  Target  Store  Area  MTI  is  most  easily  implemented  with 
a  small  digital  computer;  the  target  store,  search,  and  correlation 
operations  make  hardware  designs  fairly  complicated.  To  avoid  a 
delay  of  the  displayed  data  relative  to  radar  input,  two  target 
stores  are  required:  Input  (to  store  new  targets  as  they  come  in) 
and  Process  (to  search  for  correlating  targets).  If  500  targets  are 
to  be  processed,  about  2K  to  4K  of  sixteen-bit  memory  is  required. 

In  the  following  sections,  simple  analytical  models  will 
be  derived  and  used  to  determine  the  performance  of  the  Area  MTI 
under  various  conditions.  In  Section  3,  the  basic  mechanism  for 
bird-aircraft  discrimination  is  explained  and  a  first-cut  design  at 
an  Area  MTI  is  suggested  Details  of  this  design  are  attended  to 
in  Section  4,  5,  and  6.  The  model  is  extended  in  Section  7  to 
describe  the  interaction  effects  of  multiple  angels  in  a  dense  clutter 
environment.  In  Section  8,  the  effect  of  operator  interaction  is 
examined. 


3.  Display  Probability 


The  performance  of  any  Area  MTI  in  angel  clutter  can  be 
modeled  very  simply  by  two  numbers.  One  represents  the  probability 
that  an  aircraft  leaves  the  correlation  bin  it  was  in  during  the  last 
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scan;  the  other  is  the  probability  an  angel  leaves  the  bin.  The 

bin  corresponds  to  a  sector  cell  in  the  Sector  Map  Area  MTI  or  to 

the  correlation  window  in  the  Target  Store  Area  MTI.  An  ideal  Area 

MTI  would  be  configured  so  that  all  aircraft  move  so  fast  that  they 

always  leave  their  original  bins,  while  birds  never  leave  their 

original  bins.  Denoting  the  aircraft  and  bird  "jump"  probabilities 

by  F,  and  P..  respectively,  the  ideal  Area  MTI  has  P,  =1,  P,,  *  0. 
ja  jb  ja  jb 

A  practical  design  objective  might  be  P  *.9,  P.,  3  .2;  90%  of  the 

ja  jb 

aircraft  and  20%  of  the  birds,  on  the  average,  leave  the  bins  from 
scan  to  scan. 

The  probability  that  a  bird  or  angel  target  is  displayed 
can  be  found  easily  by  enumerating  the  two  ways  display  can  come  about. 
First,  if  a  given  target  is  detected  by  the  radar  on  the  present 
scan  but  not  on  the  previous  scan,  the  Area  MTI  detects  no  correlation 
and  thus  displays  the  target  (It  is  assumed  here  that  there  is  no 
interaction  from  other  targets;  angel  densities  are  frequently  low 
enough  so  that  this  assumption  is  approximately  valid  The  more 
general  case  is  treated  in  Section  7.)  The  other  way  a  target  can  be 
displayed  is  if  it  is  detected  on  both  scans  and  also  jumps  out  of 
its  correlation  bin  between  the  scans.  If  aircraft  and  bird  detection 
probabilities  are  denoted  by  p  and  p,  and  detection  is  assumed 

cl  D 

independent  from  scan-to-scan,  then  the  aircraft  and  bird  display 
probabilities,  (pdig)a  and  (Pjjg^  respectively,  can  be  found: 


(!Wa  "  Pa 

(1  -  >a>  +  Pa*  PJa 

(t) 

(Pdis>b  ‘  Pb 

<l  -  V  +  V  pjb 

(2) 

Figure  2  shows  how  the  display  probability  depends  on 
detection  and  jump  probabilities,  which  are  a  function  of  bin  size 
relative  to  angel  and  aircraft  motion  between  scans  and  the  target 

position  within  the  initial  bin.  For  a  properly  selected  bin  size, 

aircraft  would  tend  to  fall  on  one  of  the  curves  with  high  P,  (0.8 

to  1-0);  birds  on  a  curve  with  low  Pj  (0.0  to  0.2).  It  is  3 

interesting  to  note  that  an  Area  MTI  which  is  effective  against  birds 
(Pj  *  0.0)  discriminates  better  with  a  of  1.0  than  with  a  of 

0.5.  This  corresponds  to  the  fact  that  stationary  targets  are 
displayed  only  when  a  miss-hit  pattern  occurs;  this  pattern  is  most 
likely  at  =  0.5. 
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Figure  3  illustrates  the  effect  of  an  Area  MTI  on  aircraft 
and  bird  display  probabilities.  Both  curves  show  display  probability 
with  Area  MTI  versus  display  probability  without  Area  MTI.  Note 
that  in  all  cases,  the  Area  MTI  reduces  display  probability  (all 
points  are  below  the  45°  line).  The  significant  point  about 
Figure  3  is  that  birds  with  a  high  display  probability  before  Area 
MTI  turn  out  with  much  lower  display  probability  after  Area  MTI;  for 
example,  80%  display  probability  is  reduced  to  29%  by  Area  MTI. 

This  means  that  birds  which  single-scan  processing  fails  to  remove 
are  removed  by  Area  MTI.  When  single-scan  processing  is  effective, 
i.e  ,  =  .2  or  less,  Area  MTI  does  little  to  reduce  this  probability 

further.  Thus,  Area  MTI  might  be  viewed  as  insurance  on  single  scan 
processing  which  becomes  effective  precisely  when  the  latter  fails. 


The  single-scan  detection  probabilities,  p^  and  p^  are 
related  to  each  other  by  means  of  the  single-scan  thresholding 
technique  used  prior  to  the  Area  MTI.  Desirable  single-scan  processing 
techniques  maximize  p  while  minimizing  p,  ,  Figure  2,5-6  of  Reference 

cl  D 

2  in  effect  shows  p^  versus  1  -  p^  for  three  detectors:  an  amplitude- 

only  threshold,  a  maximum  run-length  detector,  and  an  M/N  detector. 

The  values  of  p  and  p  are  average  values  for  the  bird  and  aircraft 

cl  D 

tracks  analyzed.  To  illustrate  the  relationship  between  single-scan 

processing  effectiveness  and  Area  MTI  effectiveness,  we  will  consider 

P  *  (p,  )**  as  a  reasonable  approximation  to  the  measured  points  of  the 
a  b 

first  of  these  three  curves.  This  function  is  shown  in  Figure  4  as  a 
dashed  line.  (Actually,  what  is  shown  as  a  dashed  line  is  the  display 
probabilities  for  an  Area  MTI  which  passes  all  aircraft  and  birds. 

From  equations  (1)  and  (2),  if  P^a  =  P  ^  a  1,  then  (pdis)a  *  Pfl  an<* 

(P,,  ),  *  p.  .  Thus,  the  meaning  of  the  dashed  curve  Is  chat  with  no 
dis  D  D 

Area  MTI,  both  display  probabilities  are  low  for  very  effective  single¬ 
scan  processing  which  detects  few  birds  and  few  aircraft;  as  single  scan 
processing  is  weakened,  thus  passing  more  birds  and  more  aircraft  on 

to  the  Area  MTI,  the  locus  of  (P,  )  uad  (P,.  ),  follows  the  dashed 

,  dis  a  dis  b 

path. 


If  an  Area  MTI  is  used,  (P,,  )  and  (P,,  ),  are  related  to 

dis  a  dis  b 

p  and  p  by  equations  (l)  and  (2).  If  a  particular  set  of  Area  MTI 

2  D 

parameters  is  chosen,  P,  and  P,,  ,  the  jump  probabilities,  are  fixed 

Ja 

so  that  a  particular  curve  (P^  )  vs  (P^  *s  traced  out  as  the 


.2  .4  .6  .8 

(5*0**).  (WITHOUT  AREA  MTl) 


MT1  ON  AIRCRi 
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first  threshold  is  varied  parametrically.  Thus,  each  pair  (1* 

ja , 

P^)  defines  a  different  performance  curve  which  can  be  used  to 

judge  various  Area  MT1  parameters  using  display  probabilities  as 
criteria.  Several  of  these  curves  are  shown  in  Figure  4. 


For  example,  it  can  be  seen  that  any  Area  MTI  that  treats 
aircraft  and  birds  equally  (i-e.,  P^a  =  will  actually  degrade 

performance  relative  to  no  Area  MTI  at  all.  This  is  proved  in 

Enclosure  II. As  an  example,  the  case  of  P  =  P.  =  .5  is  shown. 

ja  Jo 

Every  point  on  this  curve  lies  downward  and  to  the  right  of  the 
dashed  curve;  in  other  words,  a  (.5,  .5)  Area  MTI  will  decrease  the 
aircraft  display  probability  and  increase  the  angel  display 
probability. 

An  ideal  Area  MTI  (P  =  1,  P^  =  0)  always  does  better 

than  no  Area  MTI  at  all;  its  points  lie  up  and  to  the  left  of  the 
dashed  line. 


A  practical  Area  MTI  design  might  achieve  P 


0.8, 


jb 


0.2.  Figure  4  indicates  that  this  Area  MTI  improves  on  single 


scan  processing  for  some  values  of  the  single-scan  threshold,  and 
does  worse  for  other  threshold  values.  In  order  for  the  (.8,  .2) 
curve  to  lie  to  the  upper  left  of  the  dashed  line,  it  turns  out  that 


must  be  larger  than  66%. 


Thus  for  this  Area  MTI,  single  scan 


processing  must  pass  at_  least  66%  of  the  bird  reports  in  order  to  make 
Area  MTI  worthwhile. 


A  more  desirable  Area  MTI  would  correspond  to  the  (.9,  .2) 
curve  In  Figure  4.  In  this  case,  p,  must  exceed  24%  in  order  to 

D 

improve  on  single  scan  processing  Points  A,  B,  and  C  on  this  curve 
which  correspond  respectively  to  p^  «  .5,  .75,  and  1.0,  illustrate 


that  performance  is  improved  as  is  raised. 

point  might  be  point  B,  for  which: 

A  reasonable  design 

P  ■  .9 
ja 

(3) 

n  i=  .2 

-jb 

(4) 

p  »  .93 
*a 

(5) 

j 

\ 
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pb 

*  .75 

(6) 

^dis^a 

=  .875 

(7) 

^dis^b 

-  .27 

00 

4.  Area  MTI  Design 

Design  parameters  such  as  type  of  Area  MTI  (Sector  Map  or 
Target  Store),  correlation  box  size,  and  data  quantization  will  now 
be  selected  to  satisfy  equations  (3)  and  (4). 

The  objective  of  this  design  is  to  pick  the  smallest 
possible  bin  size  and  largest  possible  range-  angle  quantization  LSB 
size  satisfying  the  20%  jump  probability  for  birds.  The  smaller  the 
Lin,  the  more  sensitive  the  Area  MTI  will  be  to  slow-moving  aircraft. 
The  larger  the  quantization  LSB  the  fewer  bits  that  will  be  needed  to 
store  Area  MTI  data. 

A  reasonable  bird  model  would  include  random  heading  of 
birds,  a  Rayleigh  distribution  of  air  speed,  and  a  ground  velocity 
which  is  the  vector  sum  of  air  velocity  (i.e.,  bird  velocity  with 
respect  to  the  wind)  and  wind  velocity.  It  is  shown  in  Enclosure  I 
that  this  model  corresponds  to  a  joint  Gaussian  velocity  vector  with 
a  30  knot  standard  deviation  and  a  mean  velocity  in  the  direction  of 
the  wind.  Wind  speed  will  be  modeled  as  0  knots  on  a  quiet  day  and 
30  knots  on  a  windy  day. 

With  a  jointly  Gaussian  velocity  vector,  it  turns  out  that 
leaving  the  bin  in  the  angle  direction  (probability  P.)  and  leaving 

the  bin  in  the  range  direction  (probability  P^)  are  independent 

events,  so  that  the  total  probability  of  leaving  a  bin  is: 

PJb  "  PR  +  P9  -  PR  P6  <9> 

The  Area  MTI  will  be  designed  so  that  PR  and  PQ  are 
approximat  ly  equal;  in  order  chat  ■  0.2,  PR  and  P^  must  be 
about  0. 1. 
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5.  Range  Parameters 

The  correlation  bin  is  first  considered  in  the  range 
dimension  only.  Let  the  bin  length  be  A.  A  stationary  bird  has  no 
chance  of  leaving  the  bin  in  one  scan,  while  a  bird  traveling 
with  a  radial  speed  of  A/T  (where  T  is  the  scan  period,  4  seconds) 
will  always  leave  the  bin  in  one  scan.  If  the  bird  is  randomly 
located  in  the  bin,  then  velocities  between  0  and  A/T  result  in 
intermediate  jump  probabilities,  as  shown  in  Figure  5(a).  This  is 
the  case  for  a  target  with  deterministic  velocity  and  the  Sector 
Map  Area  MTI. 

If  a  Target  Store  Area  MTI  is  used  with  infinitely  fine 
quantization  of  the  range  measurement,  the  target  is  certain  to 
leave  the  bin  if  its  velocity  exceeds  A/2T,  since  it  always  starts 
out  exactly  in  the  center  of  the  bin.  Any  velocity  less  than  A/2T 
ensures  that  the  target  never  leaves  the  bin  from  scan  to  scan. 

These  conclusions  are  reflected  in  Figure  5(b). 

Finally,  if  quantization  q  is  used  with  the  Target  Store 
Area  MTI,  the  resultant  uncertainty  in  target  position  spreads  out 
the  slope  of  the  curve,  resulting  in  Figure  5(c).  This  curve  can 
be  viewed  as  intermediate  between  Figures  5(b)  and  (a). 

If  the  above  jump  probabilities  are  integrated  over  the 
Gaussian  bird  velocity  distribution  derived  above  (30  knot  standard 
deviation,  0-30  knot  mean)  then  the  average  jump  probabilities  of 
Figure  6  are  obtained.  The  solid  and  dashed  lines  represent  the 
limiting  cases  of  coarse  and  fine  quantization  (Figures  5(a)  and 
(b))  respectively.  See  Enclosure  III  for  the  mathematical  derivation 

In  order  to  achieve  the  design  goal  of  P  =  0,1,  it  is 

K 

evident  from  Figure  6  that  A/T  must  be  at  least  120  knots  and  that 
some  finite  quantization  (intermediate  between  the  solid  and  dashed 
lines)  will  be  acceptable.  For  example,  a  range  quantization  cell 
of  q  =>  410  ft.  with  a  range  bin  size  of  A  a  1230  ft,  gives  PD  =  0.031 

for  zero  mean  wind  and  Pp  *»  0.041  for  30  knot  mean  radial  wind. 

Although  these  values  (q  =  410  ft.,  A  =  1230  ft.)  give 
conservatively  low  Pp,  they  will  be  chosen  for  practical  reasons. 

The  quantization  corresponds  to  the  ASR  pulsewidth  (1/16  nmi  or 
0.833  psec)  thus  simplifying  the  radar-Area  MTI  interface.  A  =  3q 
corresponds  to  a  quantized  range  window  of  ±  2q  *  ±  1/8  nmi,  which 
is  also  easy  to  implement.  Since  q  <  A,  the  design  has  taken  on  the 
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form  of  the  Target  Store  Area  MTI.  This  is  not  offered  as  the  only 
possible  design,  but  rather  as  a  reasonable  first-cut  design. 

With  the  above  design,  an  aircraft  moving  in  a  radial 
direction  must  have  a  ground  speed  of  about  120  knots  or  more  in 
order  to  attain  a  90%  probability  of  leaving  the  bin  in  range 
and  thus  escaping  cancellation  by  the  Area  MTI.  This  can  be  seen 
from  Figure  5(c): 

Pj  0.9  for  v  A/2T  +  q/2T  =  120  knots 


Angle  Parameters 


Selection  of  quantization  and  bin  size  in  angle  is  complicated 
by  the  fact  that  angle  measurement  errors  at  each  scan  effectively 
contribute  to  the  apparent  bird  velocity  fluctuations.  At  long  range 
(i.e.,  approaching  15  miles)  these  measurement  errors  dominate  the 
effective  velocity  standard  deviation  Thus,  the  angle  parameters 
must  be  range  dependent  to  ensure  PQ  =  0.1,  which  is  the  design  goal 

of  Section  4. 


An  RMS  angle  measurement  error  ofi  =  0.2°  is  a  reasonable 

value  for  the  ASR,  This  results  in  an  effective  bird  velocity 
standard  deviation  in  the  tangential  or  angle  direction  as  shown  in 
Figure  7.  Also  shown  is  the  apparent  RMS  velocity  of  an  aircraft, 
for  which  the  only  contributing  term  is  the  angle  measurement  error 


Since  the  highest  probability  of  leaving  the  cell  in  angle 
occurs  when  the  wind  heading  is  tangential  (i.e.,  perpendicular  to 
radial)  a  worst-case  30  knot  mean  tangential  bird  velocity  is  assumed. 


The  simplest  design,  for  angle  parameters  considered  first, 
is  the  Target  Store  Area  MTI  with  infinitely  fine  quantization 
(corresponding  to  the  jump  probability  of  Figure  5(b)).  With  this 
design,  A  must  be  large  enough  so  that  A/2T  is  about  J.28  standard 
deviations  away  from  the  mean  velocity  of  30  knots,  since  a  Gaussian 
variable  has  a  10%  chance  of  being  more  than  1.28  c-units  larger  than 
its  mean.  This  is  illustrated  in  Figure  8.  The  lightly  shaded  area 
is  the  probability  that  birds  leave  the  bin  tangentially,  P^;  this 

must  be  set  to  0,1.  Thus  given  c^,  the  bird  velocity  standard 
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deviation,  the  bln  size  must  satisfy: 


*  30  knots  +  1.28  o  , 
2T  vb 


(10) 


The  resulting  A  is  given  in  Figure  9  as  a  function  of 
range.  A  varies  from  0.15  nmi  at  close  range  to  0.27  nmi  at  15  nmi< 

With  this  bin  size  A,  how  fast  must  a  tangential  aircraft 
travel  in  order  to  achieve  a  90%  probability  of  leaving  the  box? 
Figure  9  indicates  the  answer  is: 


V 

a 


_A 

2T 


+  1.28 


(ID 


With  these  values  of  V  and  A,  the  dark  shaded  area  of 

& 

Figure  8  can  be  held  to  10%.  The  resulting  aircraft  ground  speed 

is  shown  in  Figure  9;  V  varies  between  70  knots  at  close  range  and 

a 

210  knots  at  15  nmi. 


Up  to  this  point  is  has  been  assumed  that  infinitely  fine 
quantization  has  been  used.  In  fact,  angle,  quantization  of  11  bits 
(LSB  ■  0.175°)  turns  out  to  give  false  alarm  and  detection  performance 
essentially  identical  to  infinitely  fine  quantization.  Use  of  10  bits 
or  less  essentially  inhibits  detection  of  aircraft  moving  less  than 
250  knots  between  10  and  15  nmi.  Thus  11  bits  of  angle  quantization 
are  necessary. 


Effects  of  High  Angel  Densit} 


Equation  (1)  and  (2)  of  Section  3  give  the  conditional 
display  probabilities  for  aircraft  and  birds,  given  that  no  other 
birds  were  detected  in  the  target’s  new  bin  on  the  previous  scan. 

If  one  or  more  were  detected  last  scan,  display  would  be  inhibited 
this  scan.  The  unconditional  display  probabilities  are  the  right 
hand  side  of  Equations  (1)  and  (2),  multiplied  by  the  probability  that 
no  birds  were  detected  in  the  given  bin  last  scan. 


It  is  assumed  in  this  model  that  birds  bunch  randomly;  i.e., 
have  positions  completely  independent  of  each  other.  If  the  two 
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dimensional  bird  density  Is  p  birds/nmi2  before  single-scan  processing, 
then  it  is  p  birds/nmi2  after  single  scan  processing.  We  can 

now  conceptually  divide  a  bin  with  area  A  into  N  cells,  each  with 
area  A/N.  If  N  is  large,  A/N  is  incrementally  small,  so  that,  the 
probability  that  one  or  more  birds  are  in  a  cell  is  p  pfa  A/N.  The 

probability  that  there  are  no  birds  in  a  cell  is  1  -  p  p,  A/N.  The 

D 

probability  that  there  are  no  birds  in  a  bin  (of  area  A)  is  just 
N 

(1  ■  P  P|j  A/N)  which,  for  infinitely  large  N,  approaches 

[no  birds  in  bin]  =  e~P  A  pb  (12) 


Combining  this  with  equations  (1),  (2)  we  have: 

(!Wa  '  '"a  (1  -  P.)  +  Pa2  V  *  Pb  (13> 

<Pdis>b  ”  (Pb  (l  -  V  +  pb2  Pjb>  e'P  A  Pb  (U> 


Figure  10  shows  the  effect  of  varying  clutter  density 
(p  A)  on  display  probabilities.  Each  of  the  curves  represents  the 
locus  of  possible  pairs  of  aircraft  and  bird  display  probabilities 
for  constant  p  A  and  varying  first  threshold  (i.e.,  varying  p  and 

p^).  Note  that  for  p  A  0.5,  the  Area  MTI  always  does  worse  than 

no  Area  MTI  at  all  (i.e,  the  solid  lines  are  below  the  dotted  line 

meaning  the  Area  MTI  provides  a  lower  (P  ,,  )  for  the  same  (P„  ). 

dis  a  v  dis  b 

as  no  Area  MTI).  For  p  A  <_  0.1,  the  Area  MTI  is  useful,  but  only 
if  no  single-scan  processing  is  used  (i.e.,  p  »  p  «  1). 

3  D 

The  important  conclusions  to  be  drawn  from  this  model  are 
that  heavy  angel  clutter  renders  the  Area  MTI  useless  and  that  Area 
MTI  helps  in  light  clutter  only  if  no  single-scan  processing  is  used. 


jerator  Detection  with  Area  MTI 


The  Area  MTI  mathematical  model  at  this  point  (reflected 
by  equations  (12)  and  (14))  describes  the  probability  that  a  given 
bird  or  a  given  aircraft  will  appear  on  a  display  fed  by  information 
from  an  Area  MTI.  Additional  bird/aircraft  discrimination  must  be 
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supplied  by  an  operator.  Unfortunately,  even  if  an  aircraft  i£ 
displayed,  it  may  be  rejected  by  the  operator  as  a  valid  target 
if  it  is  in  a  group  of  false  alarms  due  to  birds.  This  is  called 
masking .  At  present  we  have  only  a  rough  qualitative  concept  of 
the  effects  of  bird  masking  on  operator  detection;  psychometric 
experimentation  would  be  needed  to  quantify  these  concepts.  However, 
a  (somewhat  speculative)  model  is  proposed  here  to  indicate  the 
general  trend.  This  model  is  due  to  J.  Miller  (Reference  3). 


The  basis  of  the  model  is  a  masking  function,  F(n)  having 
the  properties  0  £  F(n)  ^  1»  E(0)  *  1,  and  F(n)-*0  as  n  -*■  00 .  In  the 
masking  model,  F(n)  is  the  probability  that  the  operator  detects 
an  aircraft,  given  that  the  aircraft  is  displayed.  In  other  words, 
it  is  the  probability  that  the  aircraft  is  not  masked.  Here  n  is 
the  "effective  average  number  of  masking  birds".  For  example, 
suppose  the  operator  looks  over  an  area  Aq  on  the  PPI  (suppose 


each  bin  consisting  of  area  A).  The 

,  The  average  number 

vhich  "survive"  Area  MTI  processing  is  n  «  p  A  ^  ^jis^b*  T*lis  is 


A  consists  of  many  bins  , 
op  * 

average  number  of  birds  in  this  area  is  p  A 


the  average  number  of  birds  which  are  effective  In  masking  an  aircraft. 
The  unconditional  probability  that  an  aircraft  is  detected  by  the 
operator  is: 


Pd  *  (p^Ja  F<n)  (15) 

dis  a 

The  function  F(n)  is  unknown,  but  Reference  3  suggests  two  possible 
functions  satisfying  the  requirements:  F(n)  =  exp  (-C^n)  and 

F(n)  =  exp  (-C^n2).  We  will  consider  the  first  of  these  two. 

Thus, 

Pd  ‘  <Pdis>a  e*p  (-Cl »  VWb> 

or 

Fd  ■  <Pdls>.  ex!>  (‘C  »  A  <PdlsV  “6> 

where  C  »  C.  A  /A  atiu  sf , .  )  and  (P . .  ),  are  found  from  equations 
1  op  dis  a  dis  b 

(13)  and  (14).  Figure  11  shows  operator  detection  probability  Pd  as 
a  function  of  single  scan  aircraft  detection  probability  P  for 

varying  clutter  levels  p  A.  The  masking  constant  C  is  10  in  Figure  11. 
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Note  that  for  light  angel  clutter  (i.e.,  p  A  <  0.1),  opcimum  Pd 
occurs  for  p  =  1,  i.e.,  no  single-scan  processing.  For  heavy 

cl 

angel  clutter  (p  A  ^  0.5),  optimum  Pd  is  attained  for  p  smaller 
than  unity.  a 

Figure  11  shows  the  effect  of  operator  masking  with 
Area  MTI.  Without  Area  MTI,  similar  masking  would  occur,  yielding 

Pd  =  p  exp  (-C  p  A  p  4)  (17) 

a  a 

Curves  similar  to  Figure  11  were  found  for  the  "no  Area  MIX"  case, 
and  the  optimum  Pd  attainable  was  cross-plotted  in  Figure  12.  Also 
plotted  were  optimum  Pd  with  and  without  Area  MTI  for  light  masking, 

C  =  1.  As  was  noted  before,  the  masking  constant  C  is  unknown. 
However,  C  =  10  seems  more  reasonable  since  it  implies  that  with  0.1 
birds/bin,  the  masking  function  is  exp  (-  (0.1) (10))  =  37%,  meaning 
that  an  aircraft  has  a  63%  chance  of  being  lost  due  to  operator 
confusion.  If  C  =  1,  than  there  is  a  9%  chance  of  being  lost  due  to 
operator  confusion  with  0.1  birds/bin.  The  former  of  these  two  seems 
more  reasonable. 

Figure  12  supports  the  following  major  conclusion.  Area  MTI 
improves  slightly  on  single-scan  processing  i_f  the  angel  clutter  is 
light  (p  A  <  0.4)  and  if  masking  of  targets  due  to  ange.l  breakthrough 
is  heavy.  If  masking  is  light  and/or  clutter  is  heavy.  Area  MTI  is  of 
little  use. 
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FIGURE  12.  OPTIMUM  DETECTION  PROBABILITY  WITH  AND  WITHOUT  AREA  Mil 


APPENDIX  C-2 


ENCLOSURE  A 
BIRD  MODEL 


In  Reference  4,  an  exhaustive  survey  of  existing  angel 
literature  formed  the  basis  for  a  mathematical  model  of  angel 
returns  due  to  birds.  One  part  of  the  model  was  the  distribution 
of  angel  speed  (i.e.,  magnitude  of  the  velocity  vector),  which  was 
found  to  be  chi-square  with  six  degrees  of  freedom  in  the  absence 
of  wind: 

v 

2  — 

f(v)  =  --~-g—  e  V0,  vo  =  9  knots  (Al) 

o 

This  distribution  does  not  lend  itself  to  calculation  of 
the  distributions  of  radial  and  tangential  velocity  components. 
However,  a  slight  modification  of  f(v)  to  g(v): 

g(v)  =  Jr  -  (A2) 

permits  one  to  find  the  radial  and  tangential  components  of  angel 
velocity. 


We  assume  that  the  angel’s  heading,  <J>,  is  a  random  variable 
between  0  and  2ir  (which  is  reasonable  for  local,  non-migratory  flight). 
Thus,  defining  a  as  $-8,  where  0  is  the  angel  bearing,  we  see  that  a 
too  is  uniform  from  0  to  2*.  The  radial  and  tangential  components  of 
velocity  are  then 

v^  »  v  cos  a  (A3) 

v  *  v  sin  a  (A4) 

respectively.  Simple  transformation  of  variables  shows  that  v^  and 

v  are  both  Gaussian  with  zero  mean  and  variance  o2.  In  view  of  the 

simplicity  of  this  result,  the  alternate  bird  speed  model  (A2)  was 
adopted  with  a  ~  30  knots.  This  results  in  an  average  air  speed  of 
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=  37.6  knots  from  (A2)  which  is  consistent  with  observation 
of  angels.  (Reference  1  and  2). 

The  above  model  describes  bird  velocities  with  no  wind  effects 
However,  the  Gaussian  model  is  excellently  suited  to  include  the 
effect  of  a  mean  wind  by  adding  the  radial  (or  tangential)  wind 
component  in  as  a  non-2ero  mean.  In  addition,  these  wind  effects  must 
be  considered  in  analyzing  an  Area  KTI  which  would  be  seriously 
affected  by  additional  bird  velocity.  Thus,  if  the  wind  direction  and 
magnitude  are  B  and  W  respectively,  then: 

vr  -  v  cos  a  +  W  cos  (8-6)  (A5) 

v  <=  v  sin  a  +  W  sin  (6-0)  (A6) 

i.e.,  radial  and  tangential  components  of  the  net  angel  velocity  are 
Gaussian  with  a  mean  component  in  the  direction  of  the  wind. 
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ENCLOSUP.E  B 


EFFECT  OF  AREA  Mil  WHEN  P 


«#  p 

ja - 


It  is  shown  here  that  any  aircraft  which  treats  aircraft 
and  birds  equally  (i.e.,  P^a  *  P^)  actually  degrades  performance 

relative  to  no  Area  MTI  at  all.  We  first  prove  the  following 
theorem: 

Theorem:  Given  h(x)  and  g(x),  with  h(x)  £  x,  g(x)  £  x,  h(0)  =  g(0)  =>  0, 

h  (x)  >_  0  and  S  (x)  £  0,  then  h(g(x))  £  g(h(x)). 

•  • 

Proof:  Since  h(x)  is  concave  upward  (h  (x)  £  0),  we  must  have 

h(A  x1  +  (1-A)  x2)  £  A  h(X])  +  (1-A)  h(x2>,  0<A  <  1  (Bl) 


In  particular  let  x^  =  0  and  define  t  *  A  x^,  whence: 

h(t)  £  ~  h(x)  0  £  t  £  x  (B2) 

Similarly.  g(x)  is  concave  downward,  yielding 

g(t)  £  “  g(x)  0  £  t  £  x  (B3) 

Since  g(x)  £  x,  we  can  substitute  g(x)  for  t  in  (B2) : 

h(g(x))  <  g(x)  ^  (B4) 

Since  h(x)  £  x,  we  can  substitute  h(x)  for  t  in  (B3): 

g(h(x))  >  h(x)  (B5) 

From  (B4)  and  (B5)  we  have 

h(g(x))  £  g(h(x))  (B6) 
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We  now  let  h  be  the  function  relating  single-scan  detection 
probabilities,  =  h(pa).  This  function  is  the  well-known  Receiver 

Operating  Curve  (ROC)  which  satisfies  the  hypotheses  of  the  theorem 
for  h  (Reference  5 

We  let  g(x)  =*  x  -(1— P)  x2,  where  P  =  P  =  P^.  Clearly, 
g(x)  satisfies  the  hypotheses  of  the  theorem. 

Given  any  p  ,  then  we  have 

cl 

pb  =  h(pa) 

<Pdis>a  =  8(Pa) 

<Pdi.>b  =  g(h(pa»  >h<S<Pa» 


(Pdis>b  i  h«Pdis>a> 


(B7) 


for  the  Area  MTI,  while  with  no  Area  MTI,  we  have  (P„  )  =  p  =  h(p  ) 

dis  b  *b 


(IWa  ■  pa 


(Pdis^b  “  h(<PdlsV  <BB> 

Thus  if  P  a  P  and  a  certain  value  of  (P_.  )  is  required,  the  Area 
j **  jo  dis  d 

MTI  will  have  a  higher  (?dis)b  than  with  no  Area  MTI. 
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ENCLOSURE  C 

CALCULATION  OF  JUMP  PROBABILITY 


Appendix  A  describes  a  model  for  the  distribution  of  radial 
bird  velocity,  resulting  in  a  Gaussian  distribution  with  mean  V  knots 
and  standard  deviation  a  knots,  where  0  £  V  30  and  a  =  30. 

The  conditional  jump  probability  (given  radial  velocity) 
is  Pj(v)  and  is  shown  in  Figure  Cl  (reproduced  from  Figure  5(c). 

Superposed  is  the  radial  velocity  distribution  f(v). 


1  2  3  4  5  6  7 


Figure  Cl 


Since  the  unconditional  jump  probability  is 


.  00 

P.  13  /  P  (v)  f(v)  dv 
J  —00  ^ 


(Cl) 


we  have  only  to  integrate  over  the  seven  regions  indicated  in  Figure  (Cl) 
yielding  the  general  result: 
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SIMPLIFIED  THEORY  OF  OPERATION  FOR 
PATTERN  RECOGNIZER  TECHNIQUES 


Operation 

Basically  a  pattern  recognizer  uses  measured  data  to 
characterize  the  object  which  the  data  represents.  Normally,  the 
characterization  is  accomplished  by  selecting  among  several  possible 
alternatives  which  most  closely  fit  the  measured  and/or  extracted 
data. 

To  arrive  at  this  decision,  all  pattern  recogniation  devices 
follow  a  three  step  procedure,  illustrated  in  general  terms  by  Figure  1. 
The  sample  space  is  formed  from  the  set  of  ail  measurements  of  a  given 
target  which  will  be  used  by  'be  recognizer.  For  this  report,  all 
37  alternate-dwell  amplitude  samples  from  the  range  cell  of  the  data 
collect  matrix  containing  the  target  formed  the  sample  space.  If 
these  37  samples  are  considered  as  independent,  then  the  resultant 
sample  space  is  considered  to  be  37  dimensional  with  one  dimension 
for  each  amplitude  sample-  Each  complete  measurement  will  then  consist 
of  a  full  37  samples  which  are  viewed  as  producing  a  l^-componenC 
vector*  in  this  space.  Each  target  return  produces  a  vector  in  this 
sample  space.  Presumably,  similar  targets  will  produce  similar  vector's, 
possibly  allowing  3  partitioning  of  this  space  into  regions  containing 
only  vectors  for  each  type  cf  targe' 

Feature  Selection 

Generally,  ether  'his  spare  !s  cf  "co  high  a  dimensionality, 
or  a  •om?1ete  separation  of  target  types  is  net  possible  with  the 
measurements  alone,  s  rhat  some  sort  of  data  reduction  is  necessary. 
This  reduces  the  dimensionality  and  increases  the  separation  of  target 
types,  if  the  appropriate  reductions  are  used.  The  results  of  these 
reductions  are  called  features  and  form  another  spare,  called  feature 
space  Now,  ea'h  reduced  result  is  considered  as  'ndependent  and 


An  n-vector  is  an  n-dimensional  vector,  l.e-,  has  n  components. 
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occupying  a  separate  dimension  Thus,  the  measurements  of  the  target 
are  converted  into  features,  forming  a  single  feature  vector  for 
each  target  observed.  In  one  of  the  approaches  considered  in  this 
report,  20  alternate  dwell  samples  containing  the  target  were  selected 
to  be  used  as  features  (the  amplitude  features  set).  In  addition, 

16  other  statistical  parameters  including  maximum  amplitude,  average 
amplitude,  run  length,  etc.,  were  found  for  these  20  samples  fthe 
statistical  features  set).  These  16  additional  quantities,  along  with 
the  selected  20  target  samples,  were  then  treated  as  forming  a  36-vector 
in  feature  space  (the  combined  features  set).  For  this  example,  the 
dimensionality  of  the  feature  space  is  essentially  the  same  as  the 
sample  space  mainly  due  to  the  limited  nature  of  the  s^udy.  In  order 
to  gain  the  maximum  benefit  within  the  time  allowed,  features  were 
liberally  selected  so  as  not  to  overlook  any  possibly  significant 
parameter.  Once  the  significant  features  are  found,  the  other  features 
can  be  eliminated,  reducing  the  dimensionality  of  the  feature  space 
required  for  the  decision  making  process  and  the  complexity  of  the 
system.  However,  the  limited  scope  of  the  task  precluded  major  effort 
toward  quantitatively  identifying  the  most  significant  features,  or 
to  reducing  the  number  of  features  to  an  optimum  level 

Decision  Algorithms 

Once  these  appropriate  features  have  been  found,  *,he  recognizer 
uses  a  decision  algorithm  which  operates  on  this  feature  °pace  and 
.dentlfies  the  target  type  by  the  feature  vprtor  fo-nd  fc  each  unknown 
*arget.  Basical1 >,  the  decision  algorithm  consists  c*  a  known 
partitioning  of  the  feature  space  into  subregions.  Previous  experience 
with  the  algorithm  is  used  to  associate  each  subregion  with  a  target 
type,  and  so,  when  a  feature  vector  for  an  unknown  target  falls  within 
a  subregion,  it  is  identified  with  the  associated  target  type  For  the 
angel  clutter  reduction  study,  thp  feature  space  was  partitioned  into 
two  subregions,  one  applicable  to  angels  and  the  other  to  aircraft. 


Ideally,  if  angels  and  aircraft  were  to  separate  into  two  distinct 
groups,  they  could  be  separated  by  a  hyperplane  drawn  between  the 
two  distinct  groups,  resulting  in  a  simple  partitioning  cf  the 
feature  space.  The  term  hyperplane,  as  used  here,  is  the  generalization 
of  a  three-dimensional  plane  to  higher  dimensional  spaces.  Since  the 
dimension  of  the  feature  spaces  used  in  this  study  were  generally  larger 
than  three  (the  dimensionality  of  the  example  considered  above  was  36, 
i.e.,  20  samples  and  16  statistical  features),  hyperplane  techniques 
for  separating  target  types  were  necessary.  Obviously,  illustrations 
cannot  convey  the  appearance  of  such  high  dimensional  spaces,  so 
examples  will  be  illustrated  in  two  dimensional  spaces  where  the  easy 
generalization  to  higher  dimensions  can  be  envisioned  (a  hyperplane 
in  two  dimensions  is  a  straight  line) ,  Figure  2  illustrates  the 
situation  which  would  occur  for  two  simply  separated  target  types 
with  two  features.  Here,  sr  in  higher  dimensions,  the  decision  is 
made  by  determining  on  which  side  of  this  hyperplane  the  known  target's 
feature  vector  lies. 

Unfortunately,  it  is  not  always  possible  to  separate  target 
types  so  easily,  as  is  illustrated  in  Figure  3.  This  situation  repre¬ 
sents  the  next  level  of  complication,  in  that  the  target  types  are 
8 till  topologically  separable,  although  now  a  more  complicated  hyper¬ 
surface  is  required  to  partition  this  space.  Once  the  hypersurface  is 
known,  however,  the  decision  is  still  made  by  noting  on  which  side 
of  the  hypersurface  the  unknown  target's  feature  vector  lies.  The 
hypersurface  itself  may  have  a  rather  complicated  structure,  in  general, 
although  for  most  practical  purposes,  it  may  be  approximated  by  the 
intersection  of  several  hyperplanes.  This  is,  in  fact,  one  of  the 
procedures  rsed  in  this  report. 

There  is,  of  course,  one  more  case  to  consider j  fbat  of 
two  target  types  which  are  not  topologically  separable  in  feature 
space.  This  is  the  case  most  often  realized  in  practice  and  is 
characv irized  by  the  groupings  of  the  two  target  types  overlapping 
to  some  degree.  Figure  i  indicates  the  case  for  a  moderate  degree 
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FIGURE  3  HYPERSURFACE  SEPARATING  TWO  TARGET  CLASSES 
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of  overlap  for  two  target  types.  Here  it  is  no  longer  possible  to 
pass  a  hypersurface  between  the  two  target  types,  so  that  no  completely 
successful  decision  making  technique  can  be  devised.  Hypersurfaces 
can  be  found,  however,  which  will  either  minimize  the  number  of 
incorrect  identifications  or  use  a  three-level  logic  to  split  targets 
into  either  type  #1,  type  If 2,  or  "don't  know".  This  gives  at  least 
some  information  concerning  the  target.  Further  identification 
techniques  can  then  be  applied  to  the  "don't  know"  targets,  to  complete 
the  process.  Occasionally  this  dilemma  may  be  resolved  i, 'L  another 
measurement  is  made,  or  another  feature  extracted,  adding  another 
dimension  to  the  problem.  If  the  additional  feature  is  chosen  properly, 
the  target  types  may  be  separated  in  this  higher  dimensional  space. 

One  is  never  guaranteed  of  finding  the  required  feature.  (It  must 
exist,  however,  otherwise  the  separate  target  categories  would  not 
exist  to  begin  with) . 

Recognizer  Development 

Developing  a  pattern  recognizer  follows  two  phases,  a 
learning  phase  and  a  testing  phase.  The  learning  phase  is  normally 
performed  off  line  with  a  training  data  set  and  some  learning  algorithm 
which  will  eventually  generate  the  desired  hypersurface.  For  this 
report  the  Widrow-Hoff  technique  was  used  for  the  learning  phase, 
as  the  software  necessary  to  implement  this  technique  had  already 
been  developed  at  Sendix,  During  the  testing  phase,  the  hypersurface 
developed  in  the  learning  phase  is  used  by  a  recognizer  (as  simulated 
in  a  computer)  to  sort  out  a  new  set  of  targets  to  verify  the  proper 
choice  of  hypersurface.  In  the  testing  phase  for  the  simple  case  of 
two  regions  separated  by  a  single  hyperplane,  the  decision  procedure 
consists  of  determining  on  which  side  of  the  hyperpiane  the  test 
target's  feature  vector  lies.  Analytically  this  corresponds  to  a 
comparison  of  the  quantity: 


B.  Widrow  and  M.  E.  Hoff,  "Adaptive  Switching  Circuits",  Technical 
Report  #15553-1,  Contract  N0NR-225(24)  Standford  Electronics  Labs, 
Standford,  California,  30  June  1960 
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S  -  E  fiw.  (1) 

i-1 

to  a  threshold  WQ,  where 

(f^)  *  the  feature  vector 

(w^)  «  the  hyperplane  (weight)  vector 

Physically,  S  represents  the  distance  along  the  hyperplane 
vector  (given  by  w^)  from  the  origin  of  the  feature  space  to  the  tip 
of  the  feature  vector  (given  by  f^).  This  equation  comes  from  a 
generalization  of  the  familiar  scalar  product  used  in  vector  analysis. 

W^,  on  the  other  hand,  is  the  distance  of  the  hyperplane  from  the  origin. 
If  S  is  greater  than  Wq,  i.e., 

S  >  W  (2) 

o 

then  the  feature  vector  is  on  the  far  side  of  the  hyperplane;  if  S  is 

less  than  W  ,  i.e. , 
o’  * 

S  <  W 
o 

the  feature  vector  is  on  the  near  side  of  the  hyperplane,  providing 
a  decision  criteria.  More  complex  hypersurfaces  proceed  in  a  similar 
manner,  except  that  now,  several  thresholds  are  used,  one  for  each 
hyperplane  used  to  approximate  the  hypersurface.  Depending  on  the 
technique  used,  all  the  thresholds  need  not  be  implemented  for  each 
target. 

There  is  a  useful  alternative  interpretation  of  equation  (1) . 
This  interpretation  views  the  components  of  the  hyperplane  vector  (w^> 
as  weights.  Equation  (1)  then  correlates  the  weights  with  the  features, 
weighting  features  favorable  to  one  target  type  positive,  and  the  other 
target  type  negative.  Interestingly,  the  magnitude  of  the  weight  gives 
the  relative  effectiveness  of  the  corresponding  feature  in  contributing 
to  the  correct  target-type  identification.  Thus  these  weights  can  be 
viewed  as  a  measure  of  the  relative  merit  of  the  various  features 
used  in  the  learning  phase.  A  feature  corresponding  to  a  low  weight 
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therefore  contains  little  information  concerning  the  identification 
of  the  type  of  target,  and  would  be  of  little  value  in  practical 
applications.  Heavily  weighted  features,  on  the  other  hand,  contain 
valid  information  concerning  target  type  and  may  be  considered  as 
distinguishing  characteristics  of  the  target.  Thus,  in  addition 
to  providing  a  systematic  approach  to  utilization  of  several  sources 
of  information  to  make  target  identifications,  the  procedure  also 
yields  an  evaluation  of  the  relative  effectiveness  of  the  information 
sources. 

The  object  of  the  Bendix  contract  was  to  use  the  angel  data 
collected  by  APL,  extract  appropriate  features,  and  use  techniques 
available  at  Bendix  to  select  the  appropriate  hypersurface  or  hyperplane, 
to  produce  the  maximum  correct-identification  rate  for  angel  and  aircraft 
targets.  This  process  of  selecting  a  hyperplane  (or  hypersurface) 
is  called  "training  a  hyperplane"  (or  hypersurface) .  The  computer  is 
first  given  an  initial  hyperplane,  then  based  on  successive  iterations 
through  the  training  set  of  data,  it  moves  the  hyperplane  around  in  the 
feature  space  until  the  maximum  number  of  correct  decisions  is  made. 

If  a  perfect  separation  of  target  types  is  not  obtained  by  this  hyper¬ 
plane,  successive  hyperplanes  can  be  generated,  using  "don't  know" 
targets  from  the  previous  hyperplane,  eventually  generating  a  fairly 
complicated  hypersurface.  Once  this  hypersurface  is  found,  it  is 
then  used  in  the  pattern  recognizer  to  make  all  future  decisions. 

Changes  in  the  distribution  of  the  target  types  are  not  accounted 
for,  as  presumably,  the  training  set  was  representative  of  all  situa¬ 
tions  the  pattern  recognizers  will  encounter. 

Of  the  techniques  discussed,  normally  one  of  them  will  be 
applicable  to  almost  any  type  of  decision  problem.  Obviously,  the 
advantage  of  these  techniques  is  that  they  can  make  use  of  Information 
from  several  channels  at  the  same  time,  resulting  in  increased  decision 
making  capabilities  over  techniques  using  only  one  channel.  On  the 
other  hand,  the  obvious  problems  presented  by  the  pattern  recognition 
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technique,  is  that  a  correct  choice  of  features  must  be  made  in 
order  to  allow  the  simplest  separation  of  target  types.  In  addition, 
the  selection  of  the  appropriate  hyperplane  or  hypersurface,  requires 
fairly  sophisticated  techniques,  particularly,  when  such  large 
dimensional  spaces  are  used. 

In  general,  the  approach  can  be  effective  for  most  problems, 
provided  appropriate  features  have  been  selected,  the  hypersui^ace 
has  been  properly  trained,  and  the  training  data  is  representative. 

The  approach,  except  for  feature  selection,  is  systematic,  and  makes 
use  of  all  the  information  available  to  make  decisions,  instead  of 
relying  on  only  one  or  two  sources.  The  results  appear  quite  promising, 
and  should  be  explored  more  thoroughly,  especially  with  an  eye  toward 
optimizing  the  features  used,  and  reducing  the  complexity  of  the 
decision-making  hardware  or  software. 
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Automatic  Aircraft/Angel  Discrimination  for  Angel  Clutter  Rejection 

I  Introduction 

This  report  describes  a  series  of  experiments  performed  at  the  Bendix 
Communications  Division  which  applied  pattern  recognition  methodology  to  the 
problem  of  angel  l-adar  clutter  reduction.  Target  amplitude  versus  azimuth  data 
were  entered  into  the  Bendix  pattern  recognition  simulation  facility  to  determine: 

(1)  What  level  of  automatic  aircraft/angel  discrimination 
reliability  performance  may  be  ultimately  feasible  and 

(2)  What  level  of  performance  can  be  expected  to  be  realized 
within  the  practical  limitations  of  a  modest  processor 
installed  at  a  radar. 

Results  of  these  efforts  are  contained  in  the  next  section  of  this  report  and  are 
related,  as  requested,  to  similar  studies  reported  in  JIIU/APL  reports  MSO-F-148 
and  MSO-F-156. 

II  Summary  of  Results 

The  discrimination  experiments  conducted  extended  the  techniques  previously 
applied  in  three  areas: 

(1)  Hyperplane  decision  rules  were  trained  to  make  the  necessary 
discrimination  by  interative  algorithms, 

(2)  Unmodified  amplitude  data  were  used  as  a  portion  of  the  input 
parameter  set  to  the  decision  rule. 

A  number  of  other  data  statistics  were  combined  with  the 
amplitude  data  to  make  the  angel/clutter  discrimination. 
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These  extensions  enabled  a  hyperplane  recognizer  to  be  trained  to  perform 
perfect  decisions  (error  free)  on  the  training  sample  of  130  aircraft  and  84  angels. 
When  tested  on  the  entire  normally  occurring  set  of  259  aircraft  and  168  angels, 
this  recognizer  rejected  91%  of  the  angel  clutter  while  maintaining  a  96%  aircraft 
acceptance  probability.  By  varying  a  threshold  in  this  recognizer,  aircraft  detection 
probability  can  be  traded-off  against  angel  rejection  probability.  The  combinations 
of  detection  and  rejection  probabilities  obtained  by  this  threshold  variation  are 
plotted  as  the  highest  curve  in  figure  1.  Figure  1  includes  curves  obtained  by  API. 
studies  for  comparison  purposes.  Bendix  hyperplane  decision  rules  operated  on 
every  other  dwell  amplitude  over  an  interval  of  40  dwells  for  a  total  of  20  dwell 
amplitudes.  Thus  the  Bendix  curves  are  best  compared  to  the  APL  curves  generated 
from  36  dwell  data. 

The  significant  performance  improvement  afforded  by  hypcrplane  discrimina¬ 
tion  is  obvious  from  the  curves,  which  indicate  that  an  operational  system  could  lie 
deployed  that  rejects  between  70  and  80%  of  the  angels  without  any  aircraft  loss. 

The  lower  hyperplane  performance  curve  plots  the  detection  and  rejection 
probability  combinations  for  all  the  aircraft/angel  data  available.  This  curve  should 
be  regarded  as  representing  the  performance  that  can  be  expected  under  the  very 
worst  operating  conditions.  Thirty-three  percent  of  the  angels  in  this  data  set  are 
labeled  as  an  "anomalous  target",  treated  as  a  bird  but  consistently  detected  as  an 
aircraft.  This  "anomalous  target"  is  responsible  for  pulling  the  lower  curve  to  the 
left  of  the  upper  curve,  i.e.  for  the  large  loss  in  angel  rejection  probability.  Also, 
61%  of  the  aircraft  in  this  data  set  are  flying  tangential  to  the  radar  and  fading  in 
and  out  of  the  MTI  notch.  The  hyperplane  threshold  setting  which  produced  the 
96%  detection  rate  quoted  above  for  normal  aircraft  only  degraded  *o  90%  for  such 
tangential  aircraft.  Thus  the  difference  in  the  two  Bendix  curves  in  figure  1  is 
almost  entirely  due  to  the  "anomalous  target". 
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Since  it  is  highly  unlikedly  that  aircraft  will  fly  tangential  to  th-n  radar  61% 
of  the  time,  or  more  importantly  that  33%  of  the  time  anomalous  targets  will  be 
present,  Bendix  hyperplane  discrimination  performance  can  be  expected  to  be  at 
or  near  the  upper  curve  plotted  in  figure  i. 

Fixed  target  clutter  data  was  also  made  available  to  Bendix.  Two  hyper- 
plane  decision  rules  were  trained  to  separate  both  aircraft  and  birds  from  such 
fixed  target  clutter  using  half  of  the  data.  When  tested  with  the  entire  data  base, 
these  decision  rules  rejected  over  90%  of  the  clutter  while  losing  only  2%  of  the 
aircraft/birds.  This  clutter  rejection  performance  is  plotted  in  figure  2.  A 
hyperplane  discriminator  of  this  type  could  be  used  in  conjunction  with  the  angel 
rejecter  to  greatly  improve  PPI  legibility. 

It  should  be  stressed  that  all  the  techniques  investigated  and/or  reported 
by  Bendix  are  amenable  to  hardware  implementation  within  the  cost  constraints 
of  individual  radar  site  deployment. 

Ill  Feature  Extraction 

Bendix  was  provided  with  every  other  dwell  amplitude  of  A  PL's  data  base. 
From  this,  further  selection  was  made. 

1,  To  the  Data  Base 

In  APL's  track  3-2,  fixed  targets  (or  ground  clutter)  appear  in  the  data. 
This  data  is  always  excluded  from  bird/angel  recognition  problems  and  is  treated 
as  a  fixed  target  in  a  separate  fixed  target  classifier. 

Tracks  16-1  and  16-3  marked  "an  anomalous  target"  were  not  used  in 
training  and  treated  as  angels  in  generalization. 

The  bulk  of  APL's  aircraft  data  is  without  MTI,  so  tracks  21-1  and  22-1 
(with  MTI  and  MTI  shut-off  in  mid  run)  were  not  used  for  training. 
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2.  Dwell  Reduction 

A  threshold  of  five  was  compared  to  each  dwell  amplitude.  Twenty 
consecutive  amplitudes  after  the  first  crossing  were  used  as  input  data.  (Twenty 
dwells  are  equivalent  to  forty  dwells  in  the  original  data. )  Only  one  aircraft 
target,  the  64th  run  of  track  6-2,  was  lost  by  this  threshold.  Inspection  revealed 
that  this  run  does  eventually  exceed  the  input  threshold  level,  but  too  late  to  take 
twenty  dwells  in  the  available  data. 

This  thresholding  has  a  strobing  effect  on  the  data  which  facilitates  the 
training  and  operation  of  a  pattern  recognizer  which  operates  on  waveform  amplitude 
data.  Strobing  avoids  simulation  of  data  propagating  through  the  recognizer  during 
generalization. 

3.  Extracted  Feature  Statistics 

SLxteen  amplitude  statistics  were  generated  from  the  twenty  dwells  selected 
by  the  threshold.  These  parameters  were  Combined  with  the  twenty  dwell  amplitudes 
for  a  total  of  thirty  six  features.  The  21st  feature  is  the  unbiased  estimate  of  the 
standard  deviation  of  the  first  twenty. 

The  22nd  feature  is  the  maximum  amplitude  and  the  23rd  is  the  average 
value  of  the  first  20. 

Features  24  through  33  were  counts  derived  from  amplitude  crossings  often 
different  double  threshold  settings.  Thresholds  were  set  at  10  &  5,  15/5,  15/10, 
20/5,  20/10,  20/15,  25/5,  25/10,  25/15  and  25/20.  A  count  was  made  of  decreasing 
wave  form  amplitudes  as  they  passed  through  each  threshold.  Passage  below  the 
upper  threshold  was  counted  as  a  +  1  while  passage  below  the  lower  was  counted 
as  a  -  1. 
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Features  34,  35  and  36  were  generated  by  counting  the  number  of 
consecutive  dwell  amplitudes  which  remained  above  Hie  thresholds  8,  12  and 
16  respectively.  Only  the  first  such  "run  length"  was  calculated  for  each  scan. 

IV  A  Hyperplane  Training  Example 

In  BCD's  pattern  recognition  routines,  a  set  of  features  can  be  selected 
from  the  20  amplitude  dwells  and  the  16  extracted  statistics.  A  feature  vector 
is  one  of  the  N  dimensional  vectors  generated  from  N  features.  In  the  example 
diagrammed  in  figure  3,  average  vector  1  (AVI)  is  the  average  of  all  aircraft 
feature  vectors  without  MTI.  Area  normalizatipr  was  applied  to  each  36  dimen¬ 
sional  feature  vector.  Average  Vector  2  (AV2)  is  the  average  feature  vector  for 
area  normalized  angel  data  (excluding  runs  16-3  and  16-1). 

The  difference  vector  is  AVI  -  AV2,  and  was  amplitude  normalized  for 
plotting  purposes. 

A  fourth  curve  set,  "Hyperplane  20",  is  a  representation  of  weighting 
values  as  a  recognizer  was  trained  on  these  data.  Since  the  difference  vector  was 
used  as  a  starting  point  in  this  training,  one  of  these  overlayed  curves  closely 
resembles  the  difference  vector  plot.  In  the  first  training  iteration  -  all  data  - 
area  normalized  -  using  the  difference  vector  as  a  hyperplane  the  recognizer  made 
67  errors,  calling  7  of  the  278  aircraft  angels  and  60  angels  of  168,  aircraft. 

As  training  proceeded,  the  total  number  of  errors  were  reduced  and  became 
more  balanced  between  the  two  classes.  With  each  iteration  a  hyperplane  was  plotted 
and  on  its  movement,  speed  of  convergence  etc.  further  experiments  were  planned. 

The  hyperplane  decision  rules  employed  by  Bendix  simply  compute  the  inner 
product  between  the  input  feature  vector  and  the  weight  vector  composed  of  the  N 
weight  parameters.  If  this  inner  product  crosses  a  similarly  trained  threshold,  the 
input  feature  vector  is  decided  to  come  from  an  aircraft.  If  it  does  not  cross  the 
threshold,  tho  target  is  rejected  as  an  angel. 
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V  Discrimination  Experiments 

A  Discrimination  based  on  extracted  statistics  alone 

Angel  rejection  expe rimer!  s  were  conducted  using  the  16  extracted 
statistic  features  alone  as  input  to  the  hyperplane  decision  rule.  This  was  done 
to  determine  whether  dwell  amplitude  information  was  needed  and  to  relate  hyper¬ 
plane  decision  performance  to  APL  results  achieved  with  feature  parameters  of 
the  same  type.  It  was  concluded  that  extracted  statistic  information  alone  is 
insufficient  for  reliable  angel  clutter  rejection  and  that  dwell  amplitude  information 
must  be  included. 

The  following  describes  the  experiments  which  lead  to  that  conclusion. 

When  a  single  hyperplane  decision  rule  cannot  handle  a  discrimination  problem, 
other  hyperplane  decision  rules  can  be  added  to  it  to  assist  it  through  the  following 
"sequential"  structure.  Such  "sequential"  hyperplane  structures,  since  they  are 
nonlinear  processors,  can  handle  discrimination  problems  which  are  linearly 
inseparable.  The  angel  clutter  rejection  problem  from  Bendix  experiments  appears 
to  be  linearly  Inseparable,  since  a  sequential  structure  was  necessary  to  achieve  zero 
errors  in  a  training  sample.  To  create  a  "sequential"  structure,  the  threshold  of 
the  hypcrplane  is  increased  and  decreased  until  no  errors  occur.  The  two  extreme 
threshold  settings  which  eliminated  the  errors  are  termed  a  double  threshold.  All 
those  feature  vectors  which  cause  hyperplane  innerproduct  values  to  lie  between 
the  double  threshold  are  difficult  cases  which  that  hyperplane  cannot  discriminate. 
Therefore,  a  second  hyperplane  is  trained  to  discriminate  just  those  feature  vectors 
which  proved  difficult  for  the  first  one.  Because  the  second  hyperplane  does  not 
have  to  consider  the  feature  vectors  that  were  successfully  handled  by  the  first 
hyperplane,  the  "difficult"  discrimination  problem  is  "easier"  for  it.  Should  the 
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second  hyperplane  decision  rule  also  not  be  able  to  achieve  zero  error  discrimination 
without  a  double  threshold,  the  procedure  is  concatenated  by  establishing  a  double 
threshold  for  it  and  training  a  third  to  discriminate  its  difficult  cases.  During 
operation,  if  the  first  hyperplane  can  make  a  decision,  i.e.  if  both  or  neither  of 
the  double  thresholds  were  crossed,  it  does.  If  not,  then  and  only  then  the  sub¬ 
sequent  hyperplanes  are  consulted  to  resolve  the  more  difficult  cases. 

Figure  4  documents  the  training  of  the  first  hyperplane  using  the  16 
extracted  statistic  features.  The  two  plots  to  the  right  of  the  figure  depict  the 
gradual  changes  in  the  hyperplane  weights  during  each  successive  training  iteration. 
Since  the  dimensionality  of  the  feature  vector  is  16,  the  weight  vector  has  a  dimen¬ 
sionality  of  16  and  therefore  has  16  components.  Wien  the  weight  corresponding 
to  a  feature  value  increases  in  absolute  magnitude,  this  indicates  that  the  decision 
rule  trainer  is  finding  out  that  this  feature  is  gaining  significance  in  the  overall 
discrimination  problem.  On  the  other  hand,  when  the  magnitude  of  a  weight 
decreases,  it  can  be  concluded  that  the  corresponding  feature  value  is  not  as 
significant  as  the  other  features  when  used  ,ln  combination  with  those  features. 

Note  that  the  14th  weight  continually  decreases  in  magnitude,  indicating  that  the 
14th,  viz.  the  run  length  over  a  threshold  setting  of  8,  could  probably  be  eliminated 
without  effect  on  the  performance  of  this  decision  rule.  This  hyperplane  was  termed 
"No.  33".  Hyperplane  numbers  are  assigned  in  the  Bendix  pattern  recognition 
simulation  facility  as  names,  and  do  not  reflect  either  chronologic  events  or  any 
other  information.  Indeed,  they  rarely  are  consecutive. 

Hyperplane  33  in  this  case  could  not  achieve  zero  error  operation. 
Therefore,  Hyperplane  36  was  trained  on  its  difficult  cases,  which  is  documented 
in  figure  5.  Hyperplane  36  made  errors  also.  A  third  level  In  the  sequential  logic 
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was  therefore  attempted.  However,  it  was  noted  in  the  average  vector  plots  at 
the  beginning  of  that  training,  which  are  plotted  in  figure  6,  that  the  aircraft 
and  angel  average  vectors  were  almost  identical  both  in  shape  and  amplitude. 

This  indicates  that  the  information  contained  in  these  16  features  is  inadequate 
for  discriminating  difficult  aircraft  and  angels.  This  conclusion  was  further 
verified  when  the  attempted  training  of  a  third  hyperplane  ended  without  signifi¬ 
cantly  increasing  the  number  of  feature  vectors  successfully  discriminated. 

B  Training  and  generalization  of  a  single  hyperplane  decision  rule 
with  36  features 

The  20  consecutive  dwell  amplitudes  after  the  first  strobe  threshold 
crossing  (i.e.  amplitude  greater  than  5)  were  then  added  to  the  id  extracted 
statistic  features  to  form  a  feature  vector  with  a  dimensionality  oi  36.  Half  of 
the  aircraft  feature  vectors  (130)  and  half  of  the  angel  vectors  (84)  were  picked 
at  random  for  hyperplane  training.  Hyperplane  No.  22  was  trained  on  this  data. 

Its  weights  are  plotted  in  figure  7.  Note  that  the  16  extracted  statistic  features 
at  the  end  of  the  feature  vector  still  appear  to  be  more  significant  than  the  dwell 
amplitudes.  The  training  of  this  hyperplane  is  documented  in  figure  Q  *ud  0. 

Using  the  threshold  produced  by  the  training  algorithm,  17  aircraft  and 
5  angels  were  misclassified  by  Hyperplane  22.  Curve  A  of  figure  10  demonstrates 
how  aircraft  detection  probability  and  angel  rejection  probability  can  be  traded-off 
through  threshold  variation.  One  of  the  points  on  this  curve  represents  the  17  ar.d 
5  errors  made  by  the  unmodified  threshold. 

Without  any  modification  of  Hyperplane  No.  22,  it  was  tested  against 
aircraft  and  168  angels.  Curve  B  of  figure  10  plots  the  probability  trade-off 
curve  for  this  larger  data  base.-  Note  how  closely  the  knee  of  this  curve  approaches 
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the  performance  achieved  on  the  training  sample  alone.  This  excellent 
"generalization”  of  Hyperplane  No.  22  guided  the  construction  of  a  sequential 
structure  based  on  it  and  described  in  the  next  section.  Point  C  in  figure  10 
is  the  generalization  of  that  sequential  hyperplane  decision  rule.  It  is  plotted 
here  simply  for  reference  and  to  demonstrate  that  a  sequential  structure  can 
improve  generalization. 

Generalization  can  sometimes  be  facilitated  by  training  a  hyperplane 
decision  rule  on  "hand  selected"  data  which  are  very  representative  of  target 
classes  to  be  discriminated.  A  training  sample  o ?  87  aircraft  and  85  angels 
was  so  selected.  Hyperplane  No.  38  plotted  in  figure  11  was  trained  on  that 
data  base.  Note  again  that  the  extracted  statistic  features  are  still  relatively 
more  significant  than  the  dwell  amplitudes.  The  training  of  this  hyperplane  is 
documented  in  figure  12.  The  generalization  probability  trade-off  curve  of 
Hyperplane  No.  38  is  plotted  in  figure  13.  Note  through  comparison  with  curve 
B  of  figure  10  that  the  generalization  of  this  hyperplane  is  inferior  to  that  of  No. 

22.  The  knee  in  this  curve  is  lower  and  less  pronounced.  In  this  case,  a  selected 
training  sample  actually  hindered  generalization. 

Based  on  these  results,  a  sequential  structure  was  trained  using  No.  22 
as  the  first  hyperplane  decision  rule. 

C  Training  and  generalization  of  a  four  hyperplane  decision  rule  with 
36  features 

Using  the  same  training  sample  of  randomly  selected  feature  vectors, 
three  additional  hyperplanes  were  trained  to  achieve  a  sequential  structure.  These 
three  hyperplanes  were  numbered  23,  24  and  26  respectively.  Their  training  is 
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documented  in  figure  14,  15  and  16.  Note  that  the  average  feature  vectors  remain 
distinctive  in  the  20  dwell  amplitude  values.  This  indicates  why  these  higher 
dimensional  feature  vectors  are  more  effective  in  providing  sufficient  information 
for  the  reliable  discrimination  of  "difficult"  aircraft  and  angels. 

Figure  17  plots  the  weight  vectors  of  all  four  hyperplane  decision  rules. 

Note  that  in  the  first  two  hyperplanes  in  the  sequential  structure,  viz.  Nos.  22  and 
23,  the  16  extracted  statistic  features  are  still  relatively  more  significant  than  the 
dwell  amplitude  information.  However,  in  the  latter  two  hyperplanes.  Nos,  24  and 
26,  all  36  features  contribute  more  or  less  equally  to  the  final  discrimination. 

This  result  confirms  the  conclusion  that  the  dwell  amplitude  information  is  essential 
for  discrimination  of  "difficult"  aircraft  and  angels.  Reliable  angel  clutter  reduction 
can  only  be  achieved  through  inclusion  of  this  type  of  input  data. 

The  first  three  hyperplanes.  Nos.  22,  23  and  24,  all  have  double  thresholds. 
If  an  innerproduct  value  for  any  of  these  three  decision  rules  lies  between  the  double 
thresholds,  that  hyperplane  does  not  make  a  decision,  but  refers  the  feature  vector 
to  the  next  hype  zplane.  However,  the  last  hyperplanc,  No.  26,  does  not  have  a 
double  threshold.  It  is  forced  to  make  a  decision  on  all  feature  vectors  which  are 
presented  to  it.  Which  hyperpiane  ultimately  makes  a  decision  in  a  sequential 
structure  like  this,  can  be  used  to  denote  the  inherent  credence  in  the  decision. 

If  the  first  hyperplane  makes  the  decision,  the  discrimination  can  be  regarded  as 
firm.  On  the  other  hand,  if  the  third  or  fourth  hyperplaneo  make  the  discrimination, 
less  credence  might  be  placed  in  the  result.  No  attempt  was  made  during  these 
efforts  to  incorporate  this  type  of  credence  information  into  the  angel  rejection 
problem. 
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Of  prime  importance  is  tlk  result  that  this  4  hyperplane  recognizer 
achieved  error  free  operation  on  its  randomly  selected  training  sample.  This 
perfect  performance  is  plotted  as  Point  A  in  figure  18.  When  generalizing  to 
the  large  sample  of  259  aircraft  and  168  angels,  Curve  C  resulted.  Curve  C 
constitutes  a  major  result  of  this  effort  and  is  discussrd  in  the  Summary  of 
Results  section  of  this  report.  The  probability  trade-*  if  to  generate  this  curve 
was  performed  by  shifting  the  double  threshold  of  Hyperplane  No.  22,  the  first 
one  in  the  sequential  chain.  This  threshold  variation  could  easily  be  controlled 
by  a  PPI  operator  to  optimize  the  angel  rejection  characteristics  in  the  face  of 
varying  target  and  clutter  conditions. 

Curve  C  establishes  that  significant  angel  rejection  is  possible  with 
zero  aircraft  loss.  However,  improved  performance  over  that,  given  by  Curve  C 
could  be  obtained  in  deployed  a/scema  by  implementing  a  sequential  hyperpl&ne 
decision  rule  that  had  the  benefit  of  a  large  training  sample.  Training  on  the  full 

feature  vector  set  was  not  attempted  during  this  effort. 

The  sequential  decision  rule  was  tested  with  the  "anomalous  angel"  and 

"tangential"  aircraft  feature  vectors.  Performance  against  this  data  is  plotted 
as  Point  D  in  figure  18.  Ncte  that  the  "anomalous  angel  target"  was  called  an 
aircraft  by  the  hyperplane  decision  rule  over  75%  of  the  time.  Combining  this 
difficult  data  with  the  normal  generalization  sample  yields  the  probability  trade¬ 
off  curve  labeled  E  in  figure  18.  Curves  E  and  C  a.  n  replotted  in  figure  1  and 
discussed  in  the  Summary  of  Results  section  of  this  report. 

D  Experiments  with  20  dwell  amplitude  features  alone 
Once  reliable  angel  rejection  had  been  demonstrated  using  36  dimensional 
features,  attempts  were  made  to  achieve  similar  levels  of  performance  with  feature 
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vectors  of  reduced  dlmensIonaHty.  Three  reasons  motivate  this  approach: 

(1)  Fewer  features  require  less  feature  extraction  signal  processing. 

(2)  Fewer  features  permit  more  efficient  hyperplane  decision  rule 
training. 

(3)  Training  with  36  dimensional  vectors  has  indicated  that  some 
of  the  features  receive  insignificant  weights  during  training. 

Hyperplane  No.  5  plotted  in  figure  19  was  trained  on  a  selection  of  167 
aircraft  and  168  angels  using  only  the  20  dwell  amplitudes  as  features.  This 
training  is  documented  in  figure  20.  Note  that  the  general  shape  character  of 
the  difference  vector  has  been  carritd  over  into  the  final  hyperplane  weights. 

Figure  21  plots  the  probability  trade-off  curves  for  this  hyperpiane.  Curve  A 
foi  this  hyperplane  plots  the  probabilities  obtained  when  tested  with  the  training 
sample.  Comparison  with  the  similar  curve  for  Hyperplane  No.  22,  Curve  A 
of  figure  10,  shows  that  No.  5  has  superior  performance,  viz.  a  higher  and  more 
pronounced  knee.  This  comparison  is  all  the  more  significant  when  realizing  that 
No.  5  has  a  much  larger  training  sample  and  uses  only  the  20  amplitude  information. 
This  indicates  that  possibly  no  statistic  feature  extraction  may  be  necessary  to 
achieve  reliable  angel  rejection.  On  the  other  hand,  comparison  of  the  B  Curves 
for  these  two  hyperplanes,  which  plot  the  probabilities  obtained  when  tested  agtinst 
all  the  normally  occurring  data,  shows  that  the  36  feature  hyperplane  No.  22, 
has  superior  generalization.  The  final  conclusion  as  to  whether  or  not  statistic 
feature  extraction  is  required  can  only  be  answered  by  the  development  of  a 
sequential  structure  based  on  20  features  similar  to  that  developed  for  the  36 
featured  vectors.  This  report  in  a  later  section  strongly  recommends  to  pursue 
that  development,  since  the  resolution  of  this  question  makes  significant  impact 
on  the  amount  of  hardware  required  for  on-site  implementation. 
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Curve  C  of  figure  21  documents  generalization  for  all  aircraft  and  angel 
data  available. 

Figure  22  plots  the  weights  of  Hyperplane  No.  37  whichwas  trained  on  the 
randomly  selected  training  sample.  Figure  23  plots  the  probabilities  obtained 
when  this  hyperpiar.e  is  tested  with  all  normally  occurring  aircraft  and  angel  data. 
Clearly  this  performance  is  inferior  to  others  previously  obtained. 

E  Fixed  target  (clutter)  rejection 

Hyperplane  decision  rules  were  trained  on  a  randomly  selected  training 
sample  of  fixed  target  (clutter)  and  aircraft/angel  feature  vectors.  Hyperplane 
Nos.  29  and  32  resulted  from  that  training  and  are  plotted  in  figures  24  and  25. 

No.  32  was  obtained  from  the  same  training  session  as  29,  but  received  more 
iterations  through  the  training  set.  This  is  reflected  in  the  weight  plot  shapes. 

No.  29  is  smoother  and  No.  32  had  more  time  to  make  finer  adjustments  to  its 
weights.  The  training  of  these  hyperplanes  is  documented  in  figure  26.  Note 
how  the  weight  plot  in  the  right-hand  portion  of  the  figure  gradually  changed 
from  that  of  No.  29  to  that  of  No.  32. 

The  excellent  generalization  properties  of  these  two  hyperplanes  are 
plotted  and  discussed  in  figure  2  presented  in  the  Summary  of  Results  section. 
Hyperplane  No.  29  had  a  97.7%  aircraft/angel  detection  rate  and  a  92.7$  fixed 
target  rejection  rate.  Comparable  percentages  for  No.  32  are  98.1$  and  90.  6%. 

Figure  27  diagrams  the  manner  in  which  this  fixed  target  clutter  rejection 
capability  could  be  used  in  conjunction  with  an  angel  rejection  function.  The 
serial  logic  flow  conserves  signal  processing  resource  in  that  only  target  data 
known  to  be  an  angel  or  an  aircraft  are  presented  to  the  angel  rejection  decision 
rule.  This  serial  logic  flow  is  permissible  because  of  the  zero-aircraft-loss 
operation  of  the  strobe  threshold  and  the  fixed  clutter  rejection  rule. 
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VI  Recommendations 

A  General  Recommendations 

The  M,ork  reported  herein  establishes  that  very  high  levels  of  angel 
clutter  reduction  are  possible  with  zero-aircraft-loss.  A  recommended  imple¬ 
mentation  is  described  in  the  next  section  of  this  report  which  is  economically 
acceptable  for  on-site  deployment. 

Therefore,  major  recommendations  can  be  made  at  thiB  time: 

(1)  Continue  hyperplane  decision  rule  experimentation  to  develop 
a  reliable  rule  trained  on  a  large  representative  sample  of  aircraft/angel 
data  using  only  the  20  dwell  amplitudes  and  as  few  of  the  special  processors 
as  possible. 

(2)  Conduct  a  systems  studv  based  on  the  results  of  this  expermen- 
tatlon  and  finalize  hardware  design. 

(3)  Construct  a  prototyc*  angel  reduction  system  and  install  at  an 
active  T.adtr  site. 

(4)  Monitor  and  evaluate  site  operations. 

(5)  Finalize  system  deployment  specification. 

(6)  Produce  and  deploy  angel  reduction  processor  systems  to  appropriate 

sites. 

Because  of  its  combined  in-depth  experience  in  all  phases  of  these  required 
tasks,  Bendix  Communications  Division  feels  especially  qualified  to  perform  them. 
The  integration  of  decision  rule  training,  systems  study,  hardware  design,  radar 
signal  processing,  site  Installation  and  equipment  production  capabilities  within  one 
contractor's  expertise  would  greatly  facilitate  action  on  these  recommendations. 
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B  Preliminary  Implementation  Recommendation 

The  implementation  shown  in  figure  28  conforms  exactly  to  that  simulated 
on  the  computer.  Thus  it  is  recommended  only  to  the  extent  that  all  of  the  various 
special  processors  are  needed.  Further  study  will  probably  show  that  only  a  few 
are  truly  useful.  However,  to  gain  some  idea  of  the  worst  case  cost  of  an  eventual 
system  this  maximum  system  was  blocked  out  and  estimates  based  on  it.  These 
budgetary,  quantity  estimates  are  under  $12,000  for  a  single  channel  system  and 
under  $23,000  for  a  dual  system. 

This  system  configuration  Is  presented  here  solely  to  substantiate  the  fact 
that  the  signal  processing  developed  during  die  reported  work  can  be  implemented 
within  the  economic  constraints  of  large  scale  site  deployment. 

In  developing  a  proposed  hardware  realization  both  minicomputers  and 
special  purpose  processors  were  considered.  The  special  purpose  circuits  won 
out  because  of  the  large  number  of  multiplications  required  in  near  real  time  and 
because  the  fixed  and  parallel  nature  of  the  operations  required  take  little  advantage 
of  the  general  purpose,  sequential  capabilities  of  a  computer. 

The  major  cost  centers  considered  in  arriving  at  quantity  budgetary 
estimates  are  listed  below: 

Material 

256  Bit  Shift  Registers  A/D  Converter 
400  Bit  Shift  Registers  Square  Root  ROM 
2048  Bit  Shift  Registers  D/A  Converters 
8  x  20  RAM's  Miscellaneous  Modules 

Multipliers 

Cabinet 
PC  Boards 
Power  Supplies 

Miscellaneous  Hardware  &  Cable  Set 


$12, 000 


Dual  Channel  System 


$23, 000 
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The  operation  of  the  system  is  as  follows: 

1.  MTI  Video  is  converted  to  5  bit  digital  words  at  a  sample  rate  of 

0. 4  usee  (approximately  one  half  the  pulsewidth).  All  processing  is  then  carried 
on  digitally.  When  video  Is  mentioned  it  should  be  understood  as  digital  video. 

2.  Whenever  the  video  exceeds  a  threshold  a  strobe  pulse  is  generated 
which  places  a  sampling  gate  around  the  range  of  the  strobe.  (Typically  the 
gate  is  3  range  cells  (1.2  usee)  wide  centered  on  the  strobe.)  This  sampling 
gate  is  repeated  every  other  PRF  period  at  the  same  range  for  40  PRF  periods. 

3.  The  sampling  gate  allows  video  to  flow  to  the  peak  detector  which 
determines  the  maximum  value  in  the  three  contiguous  range  cells.  The  output 
of  the  peak  detector  is,  therefore,  sets  of  20  samples  taken  every  other  PRF 
period  at  each  range  where  a  target  appears.  In  this  system  256  range  cells  are 
considered.  This  gives  a  total  coverage  of  1. 2  usee  x  256  *  307. 2  nsec  or  24. 8  miles. 

4.  The  samples  are  next  processed  by  the  special  purpose  circuits.  These 
handle  the  data  in  real  time  keeping  running  sums  for  each  range  increment.  For 
example,  the  average  value  circuit  (figure  29)  multiplies  each  sample  by  -|q  then 
adds  it  to  the  value  already  recirculating  in  the  shift  register  memory.  In  this 

i 

way  — -  or  S  is  available  only  one  multiplication  and  addition  after 

the  last  sample.  Multiply  and  add  times  are  typically  300  nanoseconds  each. 

Thus  die  answer  is  ready  in  less  than  die  1. 2  usee  range  increment. 

The  other  special  processor  circuits  operate  in  the  same  general 
manner  with  varying  degrees  of  complexity. 

5.  The  20  amplitude  samples  are  fed  to  the  four  recognizers  in  parallel 
where  they  too  are  processed  in  real  time.  Figure  30  shows  the  amplitude  section 
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of  a  typical  recognizer.  The  weight  memories  an  sir  Jail,  solid  state,  random 
access  'Anita  which  operate  in  300  nanoseconds.  Thus  the  time  required  for 
access,  multiply  and  isum  is  still  under  1.2  usee  keeping  the  process  real  time. 

6.  Since  die  outputs  of  the  special  processor  for  any  particular  target 
are  not  available  until  after  the  last  sampe  (the  20th) ,  they  are  weighted  and  added 
to  the  amplitude  sum  during  the  time  beyond  the  clutter  range  in  the  40th  PRF 
period. 

7.  With  the  recognition  sums  complete,  the  output  of  each  recognizer 
is  fed  to  the  logic  circuit  which  implements  the  simple  decision  tree  described 
previously  in  the  simulation  sections.  The  output  of  the  logic  circuit  for  aircraft 
is  a  pulse  at  the  target  range. 

8.  The  pulses  from  the  logic  circuit  are  sent  to  the  azimuth  width 
generator  which  produces  pulses  at  the  target  range  for  a  number  (K)  of  succeeding 
PRF  periods.  This  output  is  called  the  decision  video  and  when  applied  to  a  PPI 
will  show  only  targets  that  the  recognizer  considers  to  be  aircraft.  The  apparent 
target  azimuth  widths  can  be  varied  by  varying  K. 

9.  Provision  is  made  for  generating  gated  MTI  video  or  MTI  video  so  that 
by  using  the  PPI  mixing  circuits  various  video  and  background  combinations  can 
be  obtained. 

10.  Since  the  decision  making  process  takes  40  PRF  periods  all  videos 
and  azimuth  change  pulses  must  be  delayed  correspondingly.  Videos  which  are 
to  be  integrated  (i.e.  all  primary  radar  videos)  should  be  delayed  less  since  the 
integration  process  also  produces  a  delay. 
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C  Further  hyperplane  decision  rule  experimentation 

Work  can  proceed  Immediately  to  investigate  whether  comparable 
performance  can  be  demonstrated  with  dwell  amplitude  information  alone 
without  the  use  of  extracted  statistic  features.  If  this  can  be  accomplished, 
significant  reductions  in  signal  processing  hardware  requirements  would  be 
possible.  This  would  favorably  impact  both  system  cost  and  reliability/ 
maintainability  considerations. 

If  satisfactory  performance  cannot  be  achieved  with  dwell  amplitude 
information  alone,  surely  some  of  the  16  extracted  statistic  features  can  be 
dispensed  without  noticeable  performance  degradation.  Prime  candidates  for 
such  deletion  from  the  feature  vector  are  some  of  the  double  threshold  features 
and  possibly  the  amplitude  standard  deviation  estimate.  Any  such  reduction  in 
input  data  requirement  would  again  favorably  impact  site  implementation. 

When  using  dwell  amplitude  information,  it  has  been  demonstrated 
during  this  program  that  area  normalization  of  dwell  amplitude  feature  vectors 
can  significantly  improve  decision  rule  performance.  Amplitude  waveform 
samples  can  be  efficiently  area  normalized  in  hardware  implementations.  Thus 
experiments  should  be  conducted  which  seek  to  establish  the  benefits  accruing 
through  area  normalization. 
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DESIGN  DATA  FOR  RECEIVER  GAIN  CONTROL  CIRCUITRY 

In  order  to  accommodate  different  density  angel 

situations,  different  STC  attenuation  profiles  should  be  provided. 

Based  upon  the  information  presented  in  Section  4.1.1,  the  best 

-4 

profile  is  proportional  to  R  ;  however,  the  best  attenuation 
level  is  dependent  upon  angel  altitude.  To  provide  for  these 
variations,  the  STC  generator  design  provides  for  the  three 
profiles.  It  should  be  emphasized  that  this  circuit  approach  can 
be  used  to  generate  any  desired  profile  shape.  The  attenuation 
profile  is  maintained  in  to  2  nmi.  Within  2  nmi,  constant  atten¬ 
uation  is  used. 

Figure  E-l  is  a  schematic  diagram  of  the  STC  generator. 
The  remote  lines  select  the  appropriate  Programmed  Read-Only 
Memory  (PROM);  they  also  select  the  clock  frequency.  Since  the 
total  STC  range  is  different  for  each  curve  and  since  there  are 
64  range  bins  in  each  curve,  the  clock  frequency  is  a  function 
of  the  particular  curve  selected. 

The  control  signals  can  be  easily  traced  by  referring 
to  the  Timing  Diagram,  Figure  E-2.  The  ASR  STC  trigger  is  buffered 
by  the  input  amplifier.  The  STC  trigger  is  used  to  reset  the 
address  counters  and  the  divide  by  14,  12  and  10  counters.  One- 
shots  it 6  and  #7  generate  one  clock  pulse  during  the  radar  dead¬ 
time  and  set  the  address  counter  to  binary  "1".  The  output  word 
of  the  selected  PROM  is  then  decoded  by  the  D/A  converter.  The 
analog  voltage  at  the  converter  output  is  then  at  the  correct  DC 
level  for  the  start  of  the  sweep.  Thus  the  IF  attenuation  will 
be  constant  out  to  two  nautical  miles. 

One-shots  #3  and  #4  along  with  flip-flop  2  are  used 
to  gate  the  clock  into  the  address  counter  at  a  range  of  two 
nautical  miles.  As  the  counter  counts  from  1  to  63  the  stored 
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words  in  the  PROM  which  correspond  to  consecutive  amplitude  points 
on  the  curve  are  brought  out  on  data  lines  and  converted  to  an 
analog  voltage  by  the  D/A  converter.  The  operational  amplifier 
following  the  converter  provides  for  adjustment  of  the  STC  sweep 
amplitude. 

When  the  address  reaches  "63"  the  output  of  the  eight- 
input  NAND  gate  goes  to  ground  which  triggers  one-shot  It  5  and  its 
output  pulse  then  clears  FF2  thus  inhibiting  the  clock  pulses. 

The  contents  of  the  address  counters  remain  at  "63"  and  the  output 
of  the  D/A  converter  will  be  at  its  maximum  value.  The  next 
STC  trigger  will  generate  a  counter  reset  and  the  procedure  is 
repeated. 

The  operation  used  to  terminate  the  STC  sweep  at  a 
selected  range  operates  as  follows:  The  control  switch  is  thrown 
from  the  "normal  sweep"  to  the  "terminate  sweep"  position.  If 
the  sweep  is  to  be  terminated  at  15  nmi,  for  example,  the 
potentiometer  is  adjusted  for  a  delay  of  15  nmi  x  12.36  ua  (delay 
from  Mod  Trig  to  STC  Trig  *  197  ub) . 

At  a  range  of  15  nmi  the  trailing  edge  of  the  pulse 
at  one-shot  //2-Q  output  triggers  one-shot  It  1  and  a  1  us  pulse  is 
generated  which  prevents  FF1 .  The  ground  at  the  Q  output  of  FF1 
is  inverted  and  logic  "ones"  will  be  present  on  address  lines. 

This  address  places  amplitude  word  number  63  on  the  PROM  output 
lines,  thus  causing  the  output  of  the  D/A  converter  to  reach  its 
maximum  amplitude.  The  address  counters  are  still  counting  clock 
pulses  until  a  count  of  "63"  is  reached  and  the  eight-input  NAND 
gate  generates  a  clock  inhibit  via  one-shot  it5  and  FF2, 

The  next  STC  trigger  clears  FF1  and  resets  all  counters 
and  the  procedure  is  then  repeated. 
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AZIMUTH  CORRELATOR  DESIGN  DATA 

A  schematic  diagram  of  the  Azimuth  Correlator  is 
shown  in  Figure  E-3  and  E-4.  The  top  row  of  elements  in  Figure 
E-3  is  20  1000-bit  shift  registers.  The  outputs  are  gated  to 
a  set  of  full  adders  for  determining  the  number  of  ones  at  the 
output  of  the  shift  registers.  The  gates  are  for  range  control 
of  n.  This  sum  is  then  compared  with  the  threshold  (m)  in  the 
two  4-bit  comparators.  If  the  threshold  is  exceeded,  a  correlated 
synthetic  video  output  pulse  is  Indicated.  Figure  E-4  is  a 
schematic  of  the  control  necessary  to  change  the  values  of  n  and 
m  at  a  specified  range.  The  two  m  values  are  indicated  as  m^  and 
and  are  represented  by  5  bits  each.  When  the  contrcl  line  is 
high  one  set  of  5  is  gated  to  the  threshold.  When  the  control 
line  is  J  w,  the  other  set  of  5  are  used.  The  operation  for  n 
is  similar  except  a  decoding  matrix  is  necessary  to  provide  for 
the  proper  selection  of  the  gates.  With  4  lines  for  ^  and  n2, 
n  can  be  varied  from  5  to  20.  The  timing  diagram  is  included  in 
Figure  E-5. 
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FIGURE  E-3  AZIMUTH  CORRELATOR  DETECTOR  SCHEMATIC  DIAGRAM  (  SHEET  1  OF  2  ) 
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FIGURE  E-4  AZIMUTH  CORRELATOR  DETECTOR  SCHEMATIC  DIAGRAM 
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Affara  s-3 

ADAPTIVE  QUANTIZER  DESIGN  DATA 

An  adaptive  quantizer  is  a  video  signal  processor 
which  thresholds  each  range  cell  using  a  threshold  which  is 
derived  from  a  sample  of  the  environment  surrounding  the  range 
cell  of  interest.  The  purpose  of  the  adaptive  threshold  is  to 
maintain  a  constant  or  approximately  constant  false  alarm  rate 
in  distributed  clutter  and  receiver  noise. 

An  in-depth  study  of  the  several  forms  of  adaptive 
quantizers  suitable  for  unc  with  the  ASR  radar  was  conducted 
under  Task  I  of  this  contract*.  The  recoamended  approach,  which 
is  now  being  used  in  the  Enhanced  ARTS  Radar  Processing  Subsystem, 
is  a  rank  quantizer  (Figure  E-6) ,  which  compares  the  target 
with  a  number  of  range  cells  ahead  of  and  behind  the  target  cell. 

A  detection  is  declared  when  the  target  cell  exceeds  k  out  of  the 
total  number,  n,  of  range  calls.  The  false  alarm  rate  is  then 
independent  of  the  distribution  of  the  input  interference  (clutter, 
noise,  etc.)  and  the  processor  is  said  to  be  "distribution  free" 
or  "non-parametrlc",  with  a  false  quantization  rate  given  by: 

p  .  n  +  *  ~  fc 
n  n  +  1 

The  false  quantization  rate  can  also  be  changed  by  multiplying 
the  target  cell  voltage  by  a  constant,  K^,  for  normal  video. 

(An  additive  constant  K^,  is  used  for  log  video.)  If  is  greater 

*APL  Report  MSO-F-183,  ARTS  Enhancement  Support  Program  Multisenoor 
Sy tea  Study,  31  January  1973 
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FIGURE  E-6  RANK  QUANTIZER 
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than  unity,  the  false  quantization  rate  decreases.  Unfortunately, 
if  i*  1  (or  i*  0) ,  the  processor  looses  its  non-parametricity; 
the  false  quantization  rate  becomes  more  dependent  upon  the  inter¬ 
ference  distribution  as  the  target  cell  modifying  constants  depart 
from  these  values. 

Figures  E-7  and  E-8  are  schematics  of  an  Adaptive  Rank 
Quantizer.  In  Figure  E-7  the  radar  video  is  integrated,  sampled, 
and  A/D  converted  with  two  alternating  sets  of  hardware.  This 
is  to  allow  adequate  time  for  the  A/D  converters  to  complete  the 
operation.  The  outputs  are  multiplexed  as  a  continuous  stream  of 
6  parallel  bits.  In  Figure  E-8  the  6  parallel  bits  are  stored  in 
16  serial  hex  flip-flops.  The  outputs  of  these  flip-flops  are 
all  compared  with  the  central  flip-flop  in  digital  comparators. 

The  outputs  of  the  comparators  are  summed  in  full  adders  to  deter¬ 
mine  the  number  of  cells  exceeded  by  the  center  cell.  This  number 
is  then  compared  with  a  threshold  in  another  digital  comparator  to 
generate  the  output.  The  output  consists  of  a  single  pulse  if 
the  threshold  is  exceeded.  The  timing  diagram  is  contained  in 
Figure  E-5  of  Appendix  E-2, 
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FIGURE  E-8  ADAPTIVE  QUANTIZER  SCHEMATIC  DIAGRAM  (SHEET  2  OF  2  ) 
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SCAN  HISTORY  DISPLAY  DESIGN  DESCRIPTION 


The  Scan  History  Display  (SHD)  is  a  specialised  PP1 
built  to  display  a  maximum  of  256  target  positions  as  they 
appeared  during  the  most  recent  eight  radar  scans.  The  target 

positions  are  entered  into  the  SHD  from  the  AN/SYS-1  Target 
Information  Processing  System  developed  by  APL  for  the  U.S. 

Navy.  As  video  processor  target  detections  are  correlated  in 
the  SYS-1  software,  hit  centroid  positions  are  generated  and 
transferred  to  storage  in  the  SHD  core  memory.  For  the  radars 
normally  used  with  the  SYS-1  system,  the  target  centroids  reported 
to  the  SHD  typically  are  generated  from  six  ungated  hits  but  may 
be  generated  from  only  one  ungated  hit.  Once  the  target  positions 
are  reported  to  the  SHD  they  are  stored  and  displayed  for  eight 
scans.  The  eight  scan  display  of  target  centroids  allows  a 
greatly  enhanced  manual  detection  and  entry  capability  in 
adverse  clutter  environments.  True  target  tracks  are  easily 
discriminated  against  a  clutter  background  by  the  distinctive 
trail  they  produce. 

The  Scan  History  Display  circuitry  is  fully  documented 
in  an  AN/SYS-1  report*  which  can  be  made  available  upon  request. 
Consequently  only  a  general  outline  of  its  design  will  be  provided 
here. 

System  Block  Diagram 

The  SHD  block  diagram  is  shown  in  Figure  E-  9  indicating 
major  functional  blocks.  There  are  three  computer  input  functions 
to  the  SHD,  an  x-coordlnate  word,  a  y-coordinate  word,  and  a  zero 
bearing  mark  (Figure  E-10).  These  words  are  buffered  through  line 


*APL  Fleet  Systems  Report  SMS-FS-551,  RDS/SHD  Scan  History  Display, 
August  1972. 


£-15 


STORAGE 


E-9  BLOCK  DIAGRAM 


APPENDIX  E-4 


receivers  and  held  in  storage  buffers  until  entered  into  the 
SHD  memory  core.  Until  the  SHD  has  completely  processed  the  input 
storage  buffer  contents,  which  takes  from  18  to  30  microseconds,  a 
buffer-ready  signal  inhibits  further  computer  data  transmission. 

An  update  target  position  entered  into  the  storage  buffer  is  given 
a  pseudo  target  number  pertinent  only  to  the  current  scan.  The 
same  target  entered  in  two  successive  scans  will  probably  have 
different  target  numbers.  The  target  number  is  actually  a  partial 
address  of  where  the  x-y  coordinates  are  stored  in  the  SHD  display 
memory  core.  The  target  number  generator  is  just  a  counter 
incremented  for  each  input  buffer  update  and  is  reset  at  the  zero 
bearing  mark.  The  target  number  specifically  designates  a  sixteen 
word  scan  block  that  the  input  buffer  contents  will  be  stored. 

The  SHD  memory  is  a  4K  word  by  16  bit  magnetic  core 
logically  arranged  into  256  blocks  of  16  words.  Each  16-word 
block  is  called  a  scan  block.  Since  each  target  requires  two 
words,  eight  scans  for  256  targets  can  be  stored.  The  first  word 
in  a  scan  block  corresponds  to  a  position  in  the  current  scan,  and 
the  last  word  in  the  scan  block  corresponds  to  a  position  in  the 
eighth  elapsed  scan.  The  two-word  sections  in  a  scan  block  are 
numbered  scan  #1  through  #8.  Each  target  update  is  entered  into 
the  scan  #1  position  of  a  scan  block,  and  all  previously  stored 
coordinates  in  the  scan  number  eight  location  are  dropped.  Since 
each  scan  block  consists  of  exactly  sixteen  words,  the  pseudo 
target  number  is  actually  the  eight  most  significant  bits  of  the 
12  bit  memory  code  address  word.  The  lower  four  bits  control  a 
shift  and  write  updating  process.  Each  computer  input  during  a 
scan  is  stored  in  consecutive  3can  blocks.  The  target  number 
counter  is  reset  to  zero  at  the  radar  bearing  zero  crossing.  The 
counter  reset  action  is  actually  a  process  in  which  the  8  bit  counter 
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is  clocked  from  Its  state  at  the  time  of  zero  bearing  crossing, 
through  roll  over  to  zero.  At  each  interim  count,  the  shift  and 
update  process  is  exercised  entering  a  zero  into  the  scan  #1 
position  of  each  scan  block.  During  this  process,  lasting  a  worst 
case  of  46  milliseconds,  the  input  is  inhibited  from  accepting 
new  data. 

Between  target  updates,  stored  target  coordinates  are 
displayed.  There  exists  a  free  running  address  counter  selecting 
consecutive  target  locations  for  display.  The  counter  is  sectioned 
i~>to  three  parts  or  functions.  The  upper  eight  bits  address  con¬ 
secutive  scan  blocks  resetting  itself  every  256  counts.  Thus 
targets  are  displayed  in  the  same  order  as  they  were  received. 

The  next  lower  three  bits,  (21-23)  address  the  scan  number  and  are 
decremented  after  a  variable  number  of  roll-overs  of  the  8-bit 
scan  block  counter.  These  three  bits  specify  which  scan  number  Is 
displayed  for  each  of  the  256  target  scan  blocks.  Inhibiting 
this  counter  from  decrementing  provides  the  mechanism  for  continual 
display  of  the  same  scan  number.  The  time  duration  for  display 
of  scan  number  one  of  each  scan  block  is  an  operator  entry  and 
sets  the  time  period  that  the  three  bit  down  counter  is  inhibited 
from  changing  from  a  state  of  zero.  Controlling  the  counting 
frequency  of  this  counter  is  another  operator  entry  and  is  the 
refresh  rate  function.  The  last  bit  of  the  display  address  is  the 
x-y  bit  state.  A  logical  zero  indicates  an  x-coordinate  is  read 
and  a  logical  one  indicates  a  y-coordlnate  is  read. 

After  each  x-y  read  cycle,  the  digital  coordinates  are 
D/A  converted  and  displayed,  The  x-y  coordinate  fields  occupy  the 
lower  twelve  bits  in  each  word.  The  upper  two  bits  in  the  x-coordlnate 
word,  serves  special  functions.  No  target  position  is  displayed 
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unless  bit  number  one  in  the  x-coordinate  word  Is  set.  If  a 
track  sense  switch  is  set,  only  those  targets  with  bit  number 
two  in  the  x-coordlnate  word  set  to  one  will  be  displayed.  The 
range  scale  is  an  operator  entry  and  varies  the  gain  of  the  output 
deflection  amplifiers. 
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